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ABSTRACT 


This  work  presents  an  experimental  and  theoretical  study  of  the 
interactions  between  acoustic  surface  waves  on  a piezoelectric  insu- 
lator and  charge  carriers  in  a neighboring  semiconductor  separated  by 
a small  airgap;  and  the  application  of  these  phenomena  to  signal  pro- 
cessing and  optical  imaging.  Although  the  two  media  are  mechanically 
isolated,  the  electric  field  associated  with  the  surface  waves  pene- 
trate into  the  semiconductor  and  interact  with  its  charge  carriers. 
This  interaction  is  nonlinear  and,  in  the  case  of  two  oppositely 
propagating  surface  waves,  real  time  convolution  is  performed. 

A new  theory  of  the  convolution  interaction,  taking  into  account 
the  possible  presence  of  a depletion,  inversion  or  accumulation  la/er 
at  the  surface  of  the  semiconductor,  resulting  from  a built-' i or 
applied  dc  field,  is  presented.  Experimental  measurements  for  both 
n-  and  p-type  Silicon  on  Lithium  Niobate  systems  are  in  very  good 
agreement  with  this  theory,  and  they  lead  to  a novel  technique  of 
measurement  of  the  surface  state  distribution  on  semiconductors. 

The  Silicon  on  Lithium  Niobate  structure  used  gives  rise  to  one 
of  the  largest  convolution  efficiencies  reported  to  date.  Several 
signal  processing  applications  are  presented.  They  show  the  versa- 
tility of  the  device  which  can  be  employed  as  a programmable  filter 
or  a parametric  processor.  In  addition,  the  first  autocorrelation 
is  demonstrated  by  cascading  two  convolvers. 


•FRSCEDUO  PACK  BLANK-NOT  ? 


The  acoustic  convolver  was  adapted  for  the  scanning  of  optical 
images.  The  convolution  efficiency  can  be  modulated  by  the  local 
variations  of  the  illumination  of  the  semiconducto.'  and  the  scanning 
is  performed  by  the  acoustic  surface  waves.  We  present  a simple 
method  of  reducing  the  dark  current  to  arbitrary  level  and  hence 
increase  the  optical  dynamic  range.  The  problem  of  light  sensitivity 
is  treated  qualitatively  for  the  various  semiconductor  regimes. 

A major  advance  is  the  use  of  the  device  to  obtain  transforms  of 
optical  images.  We  have  obtained  Fourier  and  Fresnel  trarsforms  of  an 
optical  image.  A spectrum  analyzer  or  a matched  delay  line  filter  is 
used  for  the  reconstruction  of  the  original  image.  By  this  means, 
arbitrary  scan  rates  were  achieved,  large  increases  in  signal-to- 
noise  ratio  realized  and  the  bandwidth  of  the  transmitted  signal  was 
radically  reduced.  Several  examples  of  these  types  of  optical  signal 
processing  are  presented,  and  analyzed. 
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CHAPTER  I 
INTRODUCTION 


In  this  work  we  are  concerned  with  acoustic  surface  waves  on  a 
piezoelectric  substrate  and  their  interactions  with  charge  carriers 
in  a neighboring  separated  semiconductor.  Although  the  piezoelectric 
and  the  semiconductor  media  are  mechanically  Isolated  due  to  the  pre- 
sence of  an  airgap  between  the  two  surfaces,  the  coupling  is  achieved 
through  the  electric  fields  which  accompany  the  surface  waves  and 
which  extend  outside  the  piezoelectric  into  the  semiconductor. 

We  are  particularly  interested  in  the  study  of  these  inter- 
actions for  the  purpose  of  VHF  signal  processing  and  optical  imaging. 
Extenclve  work  has  been  done  in  the  past  to  develop  acoustic  wave 
amplifiers,  which  are  based  on  linear  processes.  More  recently. 


nonlinear  coupling  mechanisms  have  been  investigated,  and,  based  on 


tliese  principles,  real  time  convolution  devices,^"^  acoustic  para- 


metric amplifiers,^  and  optical  imaging  devices®^^^^®  have  been 
developed.  The  semiconductor  medium  is  either  separated  from  the 


piezoelectric  (as  in  the  works  just  cited)  or  directly  deposited  on 


it,  thus  forming  a monolithic  structure. There  is  also  a third 


family  of  devices  which  employ  piezoelectric  semiconductors  like 

1 r 1 ■«  rr  ^ 


GaAs 


15,16,17 


The  structure  which  Is  treated  throughout  this  work  is  composed 
of  two  crystalline  materials.  Firstly,  a piezoelectric  crystal,  In 


1 


1 


the  shape  o.  a long  bar,  with  a very  well  polished  surface.  High 

10 

efficiency  interdigital  transducers  are  deposited  on  the  piezo- 
electric material;  and  when  an  rf  electric  signal  is  fed  into  them, 
acoustic  stresses  are  generated  at  the  surface  of  the  acoustic 
medium  and  acoustic  surface  waves  are  launched.  Several  types  of 

surface  waves  may  exist,  and  with  the  materials  employed  here,  we 

1 A 

are  only  concerned  with  the  so-called  Rayleigh  waves.  Although 
the  vibrational  field  is  limited  to  the  piezoelectric  material,  the 
concomitant  electric  field  extends  some  distance  outside  it.  There- 
fore, if  a semiconductor  - the  second  crystalline  material  - is  placed 
close  enough  to  the  surface  of  the  delay  line,  the  electric  fields 
associated  with  the  acoustic  waves  can  penetrate  into  the  semi- 
conductor and  act  on  the  charge  carriers.  The  carriers  will  react 
to  the  electric  field.  As  a reaction  to  the  longitudinal  component 
of  the  fields,  a bunching  of  carriers  will  occur  along  the  surface, 
resulting  in  wave  attenuation.  In  addition,  the  normal  component  of 
the  electric  field  will  tend  to  deplete  or  accumulate  the  surface  of 
the  semiconductor.  While  the  potential  dropped  across  an  accumulation 
layer  is  very  small,  the  potential  generated  across  a depletion  layer 
is  large  and  it  is  proportional  to  the  square  of  the  total  electric 
field  at  the  surface.  In  other  words,  the  device  measures  the  capa- 
city of  the  depletion  layer,  which  is  a nonlinear  function  of  the 
field.  The  nonlinear  behavior  of  this  potential  is  the  basis  for 
parametric  coupling  or  mixing  of  acoustic  waves.  In  the  case  of 
oppositely  propagating  surface  waves  at  the  same  frequency,  a uniform 
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potential  is  generated  at  the  double  frequency  and,  as  a function  of 
time,  it  represents  the  convolution  of  the  modulations  of  the  input 
signals . 

Such  a device,  capable  of  taking  the  convolution  of  two  functions, 

ig 

has  very  important  applications  in  the  signal  processing  field.  ^ 

Indeed,  if  one  of  the  inputs  signals  is  regarded  as  a reference,  the 
convolution  device  may  be  considered  as  a linear  filter  whose  impulse 
response  is  the  reference  signal  and  can  ^;herefore  be  programmed  at 
will.  A matched  filter  is  particularly  simple  to  make;  one  only  has 
to  use,  as  the  reference,  the  time  inverted  replica  of  the  signal  to 
which  the  filter  has  to  be  matched.  Other  applications  include  time 
inversion  schemes  and  hence  the  possibility  of  making  a correlator, 
and  real  time  Fourier  transform. 

The  basic  nonlinear  effect,  as  we  presented  it  earlier.  Is  the 
generation  of  a depletion  layer  at  the  surface  of  the  semiconductor, 
with  an  acoustic  read-out  of  its  value.  The  thickness  of  this  deple- 
tion layer  is  a function  of  the  carrier  concentration.  Therefore,  if 
extra  charges  are  optically  produced,  the  depletion  layer  is  changed  and 
so  is  the  nonlinear  output  voltage.  By  opposition  to  a classical 
photodetector,  this  device  does  not  measure  the  photocurrent  passing 
through  it,  but  the  change  in  depletion  layer  caused  by  illumination. 

The  read-out  is  therefore  non-destructive.  Thic  principle  was  first 

g 

suggested  by  Quate  and  with  the  original  system,  the  scanning  was 
done  by  a short  acoustic  wave  pulse.  This  pulse  moved  along  the 
delay  line  at  the  acoustic  wave  velocity  and  it  sensed  changes  in 
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the  carrier  density  caused  by  local  variations  in  the  intensity  of  the 
light.  This  scanning  has  the  merit  of  being  simple  but  it  raises 
several  problems,  :.amely,  a high  inherent  dark  signal  and  a limited 
sensitivity.  In  addition,  there  is  a trade-off  between  dynamic  range 
and  resolution,  and  finally,  the  scan  velocity  is  fixed. 

The  two  phenomena  of  convolution  generation  and  optical  imaging 
are  investigated  in  this  work.  Both  n-  and  p-type  silicon  have  been 
used  with  LiNbO^  or  Bij^^^^eO^Q  piezoelectric  substrates.  Lithium 
niobate  was  chosen  for  its  strong  piezoelectricity;  Bij^^^eO^^  for  its 
slow  acoustic  surface  wave  velocity.  The  only  semiconductor  employed 
was  silicon  because  of  the  extended  knowledge  of  its  surface  proper- 
ties and  because  of  the  good  controllability  over  its  parameters. 

This  study  is  divided  into  three  parts.  In  Chapter  II,  we  pre- 
sent a new  theory  of  the  acoustic  surface  wave  convolver.  Its  opera- 
tion is  now  well  understood  and  confirmed  experimentally  in  all  regimes, 
i.e.,  depletion  at  the  surface,  at  flat  band,  with  an  accumulated  sur- 
face or  with  an  inversion  layer  at  the  surface.  This  understanding 
was  necessary  for  the  design  of  high  efficiency  devices  and  this  study 
was  also  a vital  step  in  the  development  of  sensitive  optical  imaging 
devices.  Furthermore,  previously  published  theories  only  cover  the 
flat  band  case;  we  have  shown  that  this  is  not  the  optimum  regime  of 
operation  and  moreover,  very  few  semiconductors  are  naturally  operated 
in  this  regime,  due  to  the  presence  of  localized  states  at  the  surface 
of  the  semiconductors.  Throughout  this  study,  a new  method  of  measure- 
inent  of  the  surface  state"  distribution  on  silicon  has  been  proposed 
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and  made.  It  employs  the  acoustic  field  to  measure  the  surface  proper- 
ties and  offers  the  unique  feature  of  using  oxidized  silicon  without 
requiring  deposition  of  metal  films. 

A more  device-oriented  experimental  study  of  the  convolver  is 
presented  in  Chapter  III.  Particular  emphasis  is  put  on  the  terminal 
characteristics  and  performances  of  the  device  as  a signal  processor. 
Also  comparisons  are  made  with  other  competing  acoustic  systems . We 
have  experimented  with  various  signal  processing  schemes  such  as 
convolution,  code  recognition,  time  inversion  and,  for  the  first  time, 
real  time  correlation  was  demonstrated  by  cascading  two  devices. 
Finally,  several  technological  difficulties  are  discussed:  they  are 

related  to  the  problem  of  maintaining  a uniform  and  small  alrgap 
(a  few  thousands  angstroms)  between  the  piezoelectric  and  the  semi- 
conductor media.  In  addition,  the  semiconductor  surface  is  unpro- 
tected and  has  a tendency  to  become  contaminated  if  no  precautions 
are  taken. 

In  Chapter  IV  we  show  how  the  semiconductor  convolver  was  used 
as  an  optical  imaging  device.  A new  technique  has  been  developed  to 
reduce  the  dark  signal  to  arbitrarily  low  levels  and  hence  greatly 
enhance  the  dynamic  range  and  the  sensitivity  of  the  system.  A 
quantitative  study  of  the  light  sensitivity  of  the  device  is  given 
and  compared  to  simple  theoretical  predictions . 

With  the  acoustic  semiconductor  convolver,  we  have  been  able  to 
obtain  directly,  in  real  time,  in  the  device,  Fresnel  and  Fourier 
spatial  transforms  of  an  optical  image.  The  original  Images  have 
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been  reconstructed  by  carrying  out  the  inverse  transforms.  In  the 
case  of  the  Fresnel  transform,  vre  do  this  by  inserting  the  output 
signal,  which  is  in  the  form  of  an  rf  chirp,  into  a dispersive  delay 
line.  In  the  case  of  the  Fourl  '.r  transform,  we  take  the  inve.'se 
transform  by  using  a spectrum  analyzer  or  another  acoustic  convolver. 

The  employment  of  such  signal  processing  techniques  leads  to  an  improve- 
ment in  signal-to-noise  ratio  corresponding  to  several  times  the  number 
of  spots  in  a line.  This  is  basically  because  we  are  using  the  pulse 
compression  concepts  of  radar  in  optical  imaging.  It  implies  that  we 
read  off  the  signal  from  one  single  spot  in  the  line  time  rather  than 
in  the  scan  time  of  one  spot.  This  is  an  entirely  new  technique  for 
scanning  an  optical  image.  This  same  technique  also  provides  arbitrary 
line  time  and  output  bandwidth.  Several  imaging  systems  using  these 
signal  processing  schemes  have  been  demonstrated  and  an  experimental 
and  theoretical  study  of  their  characteristics  and  limitations  is 
presented  in  Chapter  IV. 
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CHAPTER  II 


GENERAL  CONVOLUTION  THEORY  AND  EXPERIMENTAL  RESULTS 


In  this  chapter  we  shall  give  a description  of  the  structure 
of  the  semiconductor  airgap  piezoelectric  convolver.  A simplified 
theory  relying  on  the  depletion  layer  approximation  will  then  be 
described:  this  will  introduce  the  concept  of  the  varactor-like 

behavior  of  the  system  and  will  show  from  where  the  nonlinear 
field-potential  relation  arises  and  how  the  system  can  be  used  to 
take  the  convolution  product  of  electrical  signals. 

The  structure  to  which  all  the  work  reported  here  applies,  is 
illustrated  in  Fig.  2.1-1.  An  acoustic  surface  wave  delay  line  is 
made  of  a piezoelectric  material  with  a well  polished  top  surface. 
Surface  waves,  i.e.,  waves  which  decay  rapidly  away  from  the  sur- 
face, are  generated  by  applying  electrical  signals  to  interdigital 
transducers  deposited  on  the  surface  of  the  material.  Because  of 
the  piezoelectric  nature  of  the  substrate,  an  acoustoelectric  field 
is  associated  with  each  wave.  A semiconductor,  such  as  Silicon,  is 
placed  in  close  proximity  to  the  surface  of  the  acoustic  delay  line. 
The  airgap  present  between  the  semiconductor  and  the  piezoelectric 
surfaces  is  required  to  keep  the  semiconductor  from  mass  loading  the 
piezoelectric.  When  the  airgap  is  small,  the  electric  field  pene- 
trates into  the  semiconductor  and  act  with  its  carriers. 

This  surface  wave-carrier  coupling  depends  on  the  geometry  of 
the  device;  it  is  also  sensitive  to  the  surface  properties  of  the 


semiconductor  . It  can  therefore  be  controlled  by  varying  the  state 
of  the  semiconductor  surface.  This  can  simply  be  done  by  applying  a 
biasing  voltage  to  the  MIS  structure.  To  be  able  to  choose  the  opti- 
mum geometry  and  the  optimum  operating  conditions  for  the  convolver 
and  the  optical  imaging  device  derived  from  it,  a detailed  study  of 
the  interaction  process  is  needed,  and  should  be  carried  out  for  all 
surface  conditions . 

5 5 20  21 

Several  theories  have  been  published.  ''  * Although  some 

of  those  theories  are  mathematically  rigorous,  it  is  difficult  to 

understand  from  them  the  basic  physical  processes  underlying  the 

convolution  phenomenon;  moreover,  except  for  the  work  by  Takada, 

21 

et  al.,  they  only  apply  to  the  flat  band  case.  To  avoid  these 
difficulties,  we  have  developed  a somewhat  different  theory  making 
certain  assumptions,  valid  in  all  the  cases  of  practical  interest, 
which  lead  to  considerable  simplification  in  the  theory  and  a better 
physical  insight  into  the  nonlinear  phenomenon  involved.  When  the 
energy  bands  are  bent  in  the  semiconductor,  under  the  influence  of 
an  external  biasing  voltage,  simple  physical  models  are  used  with 
the  same  goal.  This  study  also  provides  an  understanding  of  the 
effect  of  the  surface  states  present  on  the  semiconductor.  The 
results  of  the  analysis  suggest  a novel  technique  for  measuring 
their  density  distribution  across  the  energy  gap  of  the  semi- 
conductor . 
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2.1  PRELIMINARY  STUDY 


2.1.1  The  Basic  Convolution  Mechanism 

We  refer  to  the  schematic  shovm  in  Fig.  2.1-1,  which  represents 
the  convolver  structure.  When  the  airgap  is  small,  the  electric 
field  associated  with  a surface  wave  can  interact  with  the  carriers 
in  the  semiconductor.  To  obtain  a more  precise  idea  of  how  this 
interaction  takes  place,  we  consider  the  simplest  regime  of  operation, 

the  "depletion"  regime  and  we  make  use  of  the  depletion  layer  approxi- 

22 

mation.  Any  electric  field  E normal  to  the  surface  of  the  semi- 
conductor produces  a depletion  layer  of  thickness  I where 
E = -qN^i/e  . is  the  donor  density  of  the  assumed  n-type  semi- 

coiuiuctor,  q the  electronic  charge  and  e the  permittivity  of  the 
semiconductor.  In  turn,  this  implies  that  there  will  be  a potential 

p 

developed  across  the  depiction  layer,  of  value  0 = -qN^i‘^/2e  or 

= -eE^/2qN^  . (2.1-1) 

In  this  case,  -e/2qN^  is  the  nonlinear  coupling  coefficient.  Hence 
a normal  component  of  electric  field  E at  the  surface  creates  a 
potential  which  is  proportional  to  the  square  of  the  field.  This  is 
the  mechanism  that  takes  place  in  a varactor;  but  whereas  a varactor 
is  generally  voltage  driven,  the  semiconductor  convolver  can  be  con- 
sidered as  a displacement  current  or  an  E-field  driven  varactor. 

When  the  device  is  operated  as  a convolver,  two  waves  are  launched 
at  frequencies  and  cd^  respectively,  in  opposite  directions  by 
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means  of  interdigital  transducers.  If  we  take  the  center  of  the 
semiconductor  as  the  origin  2 = 0 and  neglect  the  propagation  loss 
for  the  time  being,  these  waves  have  fields  that  vary  as 

f^(t-z/v)  exp  [jojj(t  - z/v)]  and  f^(t  + 7/v)  exp  + z/v)] 

respectively;  v is  the  surface  wave  velocity  and  f^^  and  f^  are 
the  modulations  of  the  input  signals.  Because  of  the  nonlinearity 
associated  with  the  semiconductor,  a potential  will  be  developed 
across  the  depletion  layer  which  is  proportional  to  the  product  of 
these  two  signals,  which  has  a variation  of  the  form 

exp  J[(cDj +0D2)t  - (cDj^ -m2)z/v]  . 

In  the  degenerate  case,  the  two  input  signals  are  at  the  same  fre- 
quency CO  = and  the  potential  0 is  uniform  along  the 

length  of  the  depletion  layer.  It  can  be  detected  between  a metal 
plane  electrode  on  the  top  surface  of  the  semiconductor  and  a plane 
electrode  on  the  lower  surface  of  the  delay  line.  In  this  case,  the 
open  circuit  output  signal  takes  the  form 

L/2 

% =^exp(2jcot)  / f^(t  - z/v)  f2(t  + z/v)d2  (2.1-2) 

4/2 

where  is  a constant  evaluated  in  Sections  2.2.2  and  2.3.I,  and 

L is  the  length  of  the  semiconduct  . If  we  suppose  that  the  modu- 
lated signals  exist  for  a time  shorter  than  the  acoustic  delay  along 
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the  semiconductor,  we  can  take  the  limits  of  integration  to  infinity; 
and  with  the  change  of  variable  t - z/v  = t , the  output  signal  takes 
the  final  form 


V 

op 


,>■/ V exp(2joit) 


/ 


-00 


f^(T)f^(2t  -T)dT 


(2.1-5; 


The  output  is  seen  to  be  the  true  convolution  of  the  modulations  of 
the  two  inputs,  compressed  by  a factor  of  2,  and  it  is  the  modulation 
of  a carrier  with  a frequency  twice  the  input  frequency. 

Since  the  semiconductor  is  uniform,  it  is  best  suited  to  detect 
spatially  invariant  potentials.  For  this  reason,  the  semiconductor 
convolver  is  generally  operated  in  a degenerate  mode,  i.e.,  with 
identical  input  frequencies.  Non-degenerate  operation  is  nevertlieless 
possible  and  of  interest. 

More  generally,  we  can  assume  that  the  coupling  coefficient 
is  nonuniform  along  the  semiconductor  t^/=  if/ {z)  . The  output 

voltage  V should  then  be  expressed  as  follows: 


^op  = exp(2jcot)  / - z/v)  f^{t  + z/v)dz  (2.1-4) 

-L/2 

Equation  (2.1-4)  provides  a different  interpretation  of  the  mathe- 
matical operation  which  takes  place  in  the  device.  Instead  of  con- 
sidering the  output  as  the  convolution  product  of  the  two  electrical 
inputs,  we  can  view  it  as  an  integral  transform  of  the  silicon-piezo- 
electric coupling  distribution  ,j/(z)  . In  this  case,  the  maximum 
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duration  of  the  input  signals  is  not  necessarily  limited  to  the 
acoustic  delay  time  along  the  semiconductor. 

We  shall  make  extensive  use  of  these  two  remarks  in  Chapter  IV 
wlien  we  describe  the  optical  imaging  device.  In  this  application, 
iV(z)  will  be  a function  of  the  illumination,  i.e.,  of  the  optical 
input  signal  . 

2.1.2  Bf.sic  Explanation  of  the  Convolution  Variation  with  Bias 
Voltage 


;o 


In  the  preceding  section,  we  showed,  in  a simple  example,  how 
the  semiconductor  can  supply  a nonlinear  mechanism  and  how  this  can 
be  used  to  perform  certain  mathematical  functions,  namely  the  con- 
volution of  two  signals. 

We  are  now  going  to  define  and  study  qualitatively  the  efficiency 
of  the  device  and  how  it  depends  on  the  material  parameters  and  the 
operating  conditions.  To  do  this  we  must  determine  the  various  mecha- 
nisms leading  to  the  generation  of  the  output  convolution  signal . The 
flow  graph  in  Chart  1 will  be  used. 

We  denote  by  and  the  input  electrical  powers.  Due 

to  the  transducer  loss,  the  corresponding  acoustic  powers  in  the  acous- 
tic waves  right  after  their  generation  are  and  P^^^<  P^^^ 

ao  e ao  e 

As  the  surface  waves  propagate,  they  are  influenced  by  the  presence  of 

the  semiconductor  and  this  will  result  in  some  propagation  loss.  The 

displacement  fields  in  the  acoustic  waves  and  penetrate 

d Si 

into  the  semiconductor  and  generate  displacement  fields  and  D 

s s 


(2) 
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CHART  I.  Flow  graph  of  the  various  mechanisms 
present  in  a semiconductor  convolver 
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at  the  surface  of  the  semiconductor.  Then,  through  the  nonlinear 
coupling  mechanism  described  in  Section  2.1.1,  these  fields  give 
rise  to  the  convolution  voltage  which  is  then  transferred  to 

a load  and  the  output  power  is  detected. 

We  can  define  several  figures  of  merit.  The  propagation  loss 
due  to  the  semiconductor  is  denoted  by  Lg  and  will  be  expressed 
in  decibels.  The  terminal  convolution  efficiency  is  defined  as 

= 10  log  . p(2)j  ^ (2.1-5) 

The  powers  are  expressed  in  nW  and  in  dBm.  Finally  the  internal 

convolution  efficiency  is  defined  as^^ 


U = 


V w 
OP 


JT^ 


(2.1-6) 


where  w is  the  acoustic  beam  width.  cY/  does  not  take  into  account 
the  transducer  loss  , the  propagation  loss  Lg  and  the  loss 

introduced  by  the  output  circuit  R . also  normalized  to  the 

power  densities  in  the  acoustic  waves.  It  will  be  seen  that  the 
following  relation  holds: 


20  - A + F.J.  -f  Lg  + + R (2.1-7) 

where 

A = 10  log  (Rj^w^)  H-  50  , 

Rj^  is  the  load  resistance.  All  the  previous  coefficients  can  be 


measured  experimentally  except  ,((  which  is  computed  with  the  help 
of  Eq.  (2.1-7). 

Three  of  the  mechanisms  shown  in  Chart  I are  due  to  the  semi- 
conductor. They  are  the  propagation  loss,  the  linear  coupling  and 
the  nonlinear  coupling.  They  will  depend  very  strongly  on  the  car- 
rier density  and  carrier  mobility  at  the  semiconductor  surface.  We 
shall  describe  these  mechanisms  in  more  detail  later. 

As  the  surface  waves  propagate  under  the  semiconductor,  the 
longitudinal  component  of  their  E-fields  will  interact  with  the 
carriers  causing  a current  J -o  flow  which  depends  on  the  semi- 
conductor resistivity.  There  will  be  a power  J . E^/2  per  unit 
volume,  absorbed  from  the  acoustic  wave  which  will  therefore  be 
attenuated.  This  attenuation  is  directly  related  to  the  resistivity 
of  the  active  part  of  the  semiconductor,  i.e.,  the  region  which  is 
actually  affected  by  the  acoustic  field.  For  very  low  resistivity, 
metal-like  semiconductors,  the  loss  will  be  small  and  it  will  in- 
crease with  the  resistivity.  On  the  other  end  of  the  spectrum,  a 
very  high  resistivity,  dielectric-like  semiconductor  will  also  not 
attenuate  the  waves  very  strongly  since  there  are  very  few  carriers. 
The  attenuatiori  is  maximum  for  an  intermediate  value  of  resistivity. 
In  addition,  the  larger  the  airgap,  the  smaller  the  coupling  to  the 
semiconductor  and  hence  the  smaller  the  attenuation. 

Each  acoustic  wave  can  be  considered  as  a voltage  source  d 
coupled  to  the  semiconductor  through  the  airgap;  the  voltage  0 
can  be  considered  as  being  placed  across  the  gap  capacitance  in 
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series  with  an  equivalent  semiconductor  imped dnce.  The  displacement 
field  at  the  surface  of  the  semiconductor,  which  corresponds  to  the 
current  flowing  into  it,  depends  on  the  value  of  tiiese  two  impedances. 
The  relation  between  the  field  at  the  semiconductor  and  the  potential 
at  the  piezoelectric  substrate  defines  the  linear  coupling  between  the 
piezoelectric  and  the  semiconductor.  The  coupling  strength  will  de- 
crease as  the  impedance  of  the  gap,  i.e.,  its  thickness,  or  that  of 
the  semiconductor  increases;  the  latter  occurs  when  the  carrier  den- 
sity is  reduced.  In  addition,  if  a thin  charge  sheet  is  present  at 
the  surface  of  the  semiconductor,  as  is  the  case  when  a dc  inversion 
layer  is  formed,  this  layer  will  tend  to  screen  out  the  free  carriers 
in  the  semiconductor  bulk  from  the  piezoelectric  field.  This  will 
result  in  a reduction  of  the  E-field  inside  the  semiconductor  and 
consequently  a reduction  of  the  linear  coupling. 

Finally,  the  convolution  output  voltage  is  a nonlinear  function 
of  the  surface  displacement  fields.  In  Section  2.1.1,  we  derived  an 
expression  for  the  nonlinear  coupling  coefficient.  If  we  write 
Eq.  (2.1-1)  in  the  form 

where  = eE^  is  the  displacement  field  at  the  surface,  the 
nonlinear  coefficient  defined  as  B = 0/D^  , is  found  to  be 

B =-l/2qN^e  . In  deriving  these  results,  we  reasoned  that  the  rf 

field  creates  an  rf  depletion  layer  across  which  an  rf  potential  is 
formed.  When  e dc  depletion  layer  exists  at  the  semiconductor  sur- 
face, the  rf  field  will  be  able  to  modulate  fully  the  depletion  layer 
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width  and  oKd  can  expect  the  above  results  to  be  exact.  This  will  be 
proved  rigorously  in  Section  2.2.2.  On  the  other  hand,  when  the  semi- 
conductor surface  is  at  flat  band  (uniform  concentration)  or  in  accumu- 
lation  (increased  majority  carrier  concentration  at  the  surface),  there 
is  no  dc  depletion  layer  and  only  half  of  each  cycle  of  the  acoustic 
field  will  deplete  the  semiconductor^^  and  create  a nonlinear  poten- 
tial. The  other  half-cycles  will  accumulate  on  the  surface  and,  be- 
cause of  the  1/N^  dependence  of  B , the  nonlinearity  associated 
with  them  will  be  negligible.  One  would  therefore  expect  a reduced 
nonlinear  coupling  in  these  circumstances. 

To  check  the  previous  qualitative  assertions,  we  measured  experi- 
mentally, the  convolution  efficiency  and  the  propagation  loss  Lg 

in  two  devices,  as  we  varied  the  dc  gate  voltage  applied  to  the 

top  electrode  of  the  semiconductor  (see  Fig.  2.1-1).  The  results  are 
plotted  in  Figs.  2. l-2(a)  and  (b)  for  p-  and  n-type  silicon, 
respectively. 

We  first  consider  a p-type  semiconductor  convolver  [Fig.  2. 1-2 (a)]. 
A positive  biasing  vc Itage  will  force  holes  toward  the  surface  and 
hence  will  accumulate  the  surface  of  the  semiconductor.  On  the  contrary 
a negative  bias  will  deplete  the  surface  of  its  holes  and,  for  a suf- 
ficiently high  bias,  an  inversion  layer  of  negative  polarity  will  be 
formed.  The  flat  band  condition  - uniform  carrier  distribution  through- 
out the  semiconductor  - corresponds  roughly  to  a zero  bias. 

^ ^ band  condition  is  nearly  satisfied.  The 

surface  resistivity  of  the  silicon  is  relatively  large  (lOO  Q-cm)  and 
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FIG.  2.1-2 — Experimental  measurement  of  convolution  efficiency  (F„) 
and  propagation  loss  (Lg)  vs  bias  voltage  (V^),  for 

(a)  an  n-type  Silicon  on  YZ  LiNbO,  convolver.  The  full 
line  curve  corresponds  to  the^dc  bias  measurements; 


the  dashed  line  curve  to  the  pulsed  bias  measurement, 


(b) 

FIG.  2.1-2— Experimental  measurement  of  convolution  efficiency  (F,,) 
and  propagation  loss  (Lg)  vs  bias  voltage  (V^ ) . for 

(b)  a p-type  Silicon  on  YZ  LiNbO^  convolver  (dc  bias  only) 
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the  propagation  loss  is  rather  high.  As  the  bias  is  increased  posi- 
tively, the  surface  is  driven  into  accumulation;  thus  the  surface 
resistivity  decreases  sharply  and  so  does  the  loss.  For  n''gative 
bias,  a dc  depletion  layer  is  first  formed,  which  acts  as  an  extra 
insulating  gap  and  reduces  the  coupling  to  the  semiconductor  and  hence 
the  loss.  As  the  bias  is  made  more  negative,  an  inversion  layer  is 
generated  and  the  surface  resistivity  decreases  monotonically  from 
almost  infinity  (insulating  depletion  layer)  to  zero.  The  loss 
exhibits  a maximum  as  we  predicted  earlier. 

As  for  the  convolution  efficiency,  it  decreases  in  the  accumula- 
tion regime  because  of  the  increased  carrier  density.  In  depletion, 
the  nonlinear  mechanism  becomes  more  efficient  because  the  fields  can 
now  fully  modulate  the  depletion  layer;  this  shows  up  as  an  increase 
in  convolution  efficiency.  Finally,  for  la.'ge  negative  bias  voltages, 
the  inversion  layer  charge  shields  the  semiconductor  and,  at  the  same 
time,  attenuates  the  waves,  both  of  these  effects  resulting  in  a 
sharply  reduced  convolution  efficiency. 

The  n-type  semiconductor  convolver  behaves  in  a similar  fashion. 

A negative  voltage  is  now  necessary  for  accumulation,  and  a positive 
one  for  depletion  and  inversion.  The  loss  at  flat  band  is  low  because 
of  the  lower  resistivity  of  the  silicon  used  here  (10  Q-cm).  In 
inversion,  the  propagation  loss  increases,  but  the  maximum  of  loss 
corresponds  to  a bias  voltage  outside  the  experimental  range  and  it 
is  not  observed  in  this  experiment. 
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A quantitative  theory  is  required  to  explain  the  exact  shape  of 
the  curves  plotted  in  Figs.  2.1-2(a)  and  (b).  It  is  developed  in  this 
chapter.  It  is  necessary  to  take  account  of  the  charges  trapped  at 
the  surface  of  the  semiconductor .^^^25-2?  According  to  the  dc  bias 
used,  these  traps  are  emptied  or  filled;  this  gives  rise  to  a surface 
state  screening  effect  so  that  higher  voltages  are  necessary  to  obtain 
inversion  than  would  be  the  case  if  no  surface  states  were  present. 

In  Fig.  2.1-2(b),  we  also  show  plots  of  the  convolution  effi- 
ciency and  of  the  propagation  loss  with  a pulsed  bias.  F.j,  and  Lg 
are  measured  right  after  the  bias  has  been  applied.  It  is  interesting 
to  note  that,  in  this  case,  when  a pulsed  bias  is  used,  no  charge  can 
be  exchanged  with  the  traps  because  of  their  long  tima  constants. ^5 
Very  large  fields  are  therefore  present  at  the  surface  and  if  there 
exists  a source  of  minority  carriers  at  the  surface  of  the  semicon- 
ductor (like  a p+  impurity  center),  a large  number  of  minority  car- 
riers will  immediately  be  injected.  These  carriers  will  be  concen- 
trated in  a thin  sheet  of  very  low  resistivity  at  the  surface,  under 
the  influence  of  the  applied  field.  The  propagation  loss  L will 

s 

therefore  be  very  small  and  the  convolution  decreases  because  of  the 
screening  introduced  by  this  minority  charge  sheet . This  seems  to  be 
the  mechanism  wl  ich  controls  the  n-type  silicon  convolver  whose  F 

T 

and  Lg  are  plotted  in  dashed  lines  in  Fig.  2. 1-2^3).  Minority 
carrier  sources  may  be  created  during  the  preparation  of  the  samples. 

In  particular,  the  samples  are  cut  with  a diamond  saw;  this  may  intro- 
duce p+  impurity  centers  at  the  extremities  of  the  semiconductor. 
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Let  us  note  however  that  other  silicon  samples  did  not  show  any  sign 
of  minority  charge  injection;  the  corresponding  remained  at 

approximately  a constant  level  when  a positive  bias  was  applied^  as 
we  would  expect  from  a deep  depletion  theory  without  any  inversion. 

2.2  COMPLETE  THEORY 

2.2.1  General  Discussion 

In  this  section,  we  are  concerned  with  a quantitative  study  of 
the  linear  and  nonlinear  coupling  between  Rayleigh  surface  wave  field 
and  carriers  in  a neighboring  semiconductor.  We  shall  evaluate  the 
convolution  efficiency and  the  propagation  loss  L_  for  all  sur- 

d 

face  conditions,  i.e.,  all  gate  voltage  V_  . The  structure  used 
is  that  of  Fig.  2.1-1. 

We  ,'ssume  that  the  properties  of  the  Rayleigh  wave  propagating 

along  the  electrically  free  surface  of  the  piezoelectric  substrate 

18 

is  known.  Because  of  the  weak  coupling  approximation,  i.e.,  the 
fact  that  the  electrical  power  is  a fraction  of  the  mechanical  power 
in  the  piezoelectric,  the  acoustic  wave  propagation  constant  is  only 
slightly  affected  by  the  electrical  loading  due  to  the  semiconductor. 
This  loading  can  be  considered  as  a small  perturbation  of  the  unper- 
turbed system  composed  of  the  bare  delay  line.  This  technique  was 

28 

developed  by  Kino-Reeder  and  we  use  their  results  to  soLve  the 
linear  coupling  problem  in  Section  2.2.5. 

To  evaluate  the  nonlinear  coupling  due  to  the  space  charge  layer 
in  the  semiconductor,  we  assume  thermal  equilibrium,  i.e.,  we  assume 
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that  the  C3i.riers  follow  the  acoustic  field  variations.  This  is  only 
a valid  assumption  for  majority  carrier  concentrations,  and  particular 
attention  must  be  paid  to  its  use  near  intrinsic  and  inversion  condi- 
tions. We  shall  discuss  this  further  in  Section  2.2.2. 

In  this  study  we  compute  the  open  circuit  convolution  output 
voltage  denoted  by  and  we  shall  refer  it  to  the  acoustic  powers 

on  the  delay  line.  This  means  that  we  do  not  take  into  account  the 
transducer  and  output  circuit  loss.  We  write  as 


V = ON.  . 
op  0 


Vq  is  the  open  circuit  output  voltage  computed  assuming  no  propaga- 
tion loss  under  the  semiconductor;  the  loss  Lg  = -20  log  CC  is  the 
propagation  loss  due  to  the  presence  of  the  semiconductor.  If  a 
normal  surface  displacement  field  is  created  at  the  surface  of 

the  semiconductor,  we  shall  see  that  a potential  0^  is  generated 
across  the  device,  which  has  the  following  form: 


0 = AD  + BD 

'^s  s s 


We  now  assur.e  two  oppositely  traveling  waves  are  present  with  fields 

and  at  the  surface  of  the  semiconductor,  with  acoustic 

s s 

potentials  at  the  surface  0^^^  and  0^^^  and  powers  and 

3d  3 3 

The  "nonlinear  convolution"  voltage  generated  at  the  sum  frequency  of 
the  input  frequency  is 


V = 

0 s s 
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It  can  be  rewritten  as 


•“  B is  the  nonlinear  coupling  coefficient;  it  is  derived  in 
Section  2.2.2  by  integration  of  Poisson's  equation. 

--  The  term 

,av(3)  = ( ,(i)j 

Is  the  linear  gap  coupling;  it  is  derived  in  Section  2.2.5  by  use  of 
an  equivalent  circuit  concept. 

--  The  term 

^u,v(2) 

— ■ ■■  , = 2Z 

Jp(l)p(2) 

* a a 

is  twice  the  acoustic  impedance  of  the  piezoelectric  used.  This  con- 
cept is  reviewed  in  Section  2.2.h. 

The  loss  term  Lg  is  a linear  term  and  hence  is  obtained  in 
Section  2.2.5. 

2 .2 .2  The  Nonlinear  Coupling 

We  define  the  potential  0 as  zero  in  the  bulk  of  the  semi- 
conductor; 0 is  therefore  measured  from  the  intrinsic  Fermi  level 
of  the  conduction  band  level  in  the  bulk.  We  consider  an  n-type 


Ik 
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semiconduccor  with  a donor  density  . At  room  temperature,  most 

of  the  donors  are  ionized,  and  the  bulk  electron  concentration  is 
; from  the  mass  action  law,^^  the  hole  concentration  in 
th.,  bulk  is  given  by  =,  Intrinsic 

concentration.  If  we  assume  non-degenerate  conditions  .ad  thermal 
equilibrium,  Boltzmann  statistics  apply  and  the  carrier  concentra- 
tions are  the  following  functions  of  the  potential  0 : 

n = exp  (q0/kT) 

P ■■  "i/Njj  exp  (-q0/1cT) 

We  then  solve  the  relation  between  and  0^  at  the  semicon- 

ductor surface  by  using  standard  MIS  diode  theory,  and  perturbing 
the  dc  solution  so  as  to  determine  0^  to  second  order  in  the  per- 
turbing displacement  field  . The  potential  as  a function  of  y 

n.ay  be  determined  by  writing  the  one-dimensional  Poisson's  equation 

d^0/dy^  = -p/e 

where  p is  the  space  charge  density  given  by 

p = q(N^  - n + p)  . 

Poisson's  equation  can  therefore  be  written  as 


exp  1 

(?) 

1 

2 

exp 

(- 1^)  - ii 

e 

- 

- 

^^d 

. 

' kT' 

When  the  surface  is  well  depleted j the  potential  is  strongly  negative 
over  the  extent  of  the  depletion  layer  and  the  charge  density  in  the 
depletion  layer  is  essentially  . In  this  case,  the  depletion 

layer  assumption  can  be  used  and  one  can  define  an  equivalent  deple- 

27 

tion  width  of  length  f over  which  the  semiconductor  is  fully 
depleted.  Equation  (2.2-2)  reduces  to 


d^0/dy^  = 

-qN^/e 

for  0 < y < £ 

0=0 

for  y > £ 

By  integration,  we  can  evaluate  the  surface  potential 

and  the  surface  displacement  field  D =-e^  (v  = 0) 

s dy  ' 


0 


s 


> 


0(y  = O) 


0 


s 


and  D = -qN^f 
s ^ d 


Elimination  of  £ yields  the  relation 


0 = -D^/2qN^e 


which  is  the  simple  result  used  previously  in  Section  2.1.1. 

In  general,  we  integrate  Eq.  (2.2-2)  once  by  multiplying  by 
d0/dy  } we  find  at  the  surface  of  the  semiconductor 


dJ  = N^kTe  F^(q0^/kT) 


(2.2-5) 


where  the  function  F is  defined^^  a 


F(v)  = 

where  v is  a dummy  variable. 


2 

n. 

exp  V - V - 1 + 

^d 


1/2 


exp  ( -V ) + V - 1 


, (2.2-U) 
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If  a u bias  is  applied  to  the  semiconductor,  we  write  the 

potential  and  the  displacement  field  as 


" ^sO  ^sl 

respectively,  where  the  subscripts  0 and  1 refer  to  dc  and  rf 
terms,  respectively.  0^  and  must  satisfy  Eq.  (2.2-3)  and  a 

general  solution  for  0^  can  be  written  as  a series  expansion  in 
the  successive  powers  of  D 

s 

^s  " ^sO  ^^°s  "^so^  ®^^s  ■*■  •••  (2.2-5) 

The  convolution  output  voltage  corresponds  to  the  second  order  non- 
linear term;  we  therefore  solve  Eq.  (2.2-5)  for  0^  to  second  order 
in  Dg  . It  follows  that 


A = 2D 


sO 


i(exp  Uq  - 1)  - b|exp(-u^)  - ij 


-1 


(2.2-6) 


B = 


i(exp  Uq  - 1)  - b(exp(-u^)  - l) 


-1 


is.  1 

kX  ' 

1 a exp 

b exp(-u^)) 

|a(exp(u^)  - l] 

1 -b| 

(exp(-u^)  - l) 

f 

(2.2-7) 


where 


“O  = ^ ; b = 2qn^eA^ 

(2.2-8) 

^sO  " • sgn(0^^) 
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For  an  n-type  semiconductor  and  flat  band  condition.  D =0 

’ sO  ' 

A = ^ “ -l/6qN^e  . When  the  surface  is  well  depleted 

A - j B = -l/2qN^e  where  d is  the  length  of  the  depletion 
layer  and  the  Debye  length.  This  is  again  the  result  derived 

in  Section  2.1.1.  When  the  surface  is  strongly  accumulated  A -> 0 
and  B ->  0 

So  far  we  have  assumed  thermal  equilibrium,  i.e.,  equilibrium, 
at  all  times,  between  the  carriers  and  the  total  electric  field.  The 
total  field  here  is  the  superposition  of  the  dc  biasing  field  and  the 
rf  acoustic  field.  To  stay  in  equilibrium,  the  carriers  must  be 
supplied  at  a rate  higher  than  the  rf  frequency.  For  a frequency 
of  100  MHz,  the  time  constant  associated  with  this  transport  of 

must  therefore  be  smaller  than  10  nsec . In  the  accumulation 
and  depletion  regimes,  the  principal  carriers  are  of  the  majority  type 
and  the  time  constant  involved  is  of  the  order  of  lo  psec;^^  so  we 
can  use  the  thermal  equilibrium  assumption.  On  the  contrary,  in  the 
case  of  a dc  inversion  layer  being  formed  at  the  surface  of  the  semi- 
conductor, the  supply  of  minority  carriers  from  the  bulk  is  very  slow, 
because  of  the  very  long  generation-recombination  times  involved. 

These  cariiers  will  not  be  able  to  follow  the  normal  component  of  the 
acoustic  E-field,  which  is  the  cause  of  the  nonlinear  potential.  The 
only  source  of  nonlinear  coupling  will  be  the  depletion  layer  located 
behind  the  inversion  layer  and  B will  remain  equal  to  -l/2qN^e 
This  is  equivalent  to  the  well  known  result  for  the  behavior  of  the 
capacity  of  MOS  diodes  with  frequency:  the  high  frequency  capacitance 
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-■  

remains  equal  to  the  depletion  capacitance  before  inversion, 

A complete  plot  of  -BqN^e  vs  q0^^/kT  is  given  in  Fig,  2,2-1,  The 
low  frequency  curves  have  been  drawn  assuming  thermal  equilibrium, 
but  as  we  have  just  discussed,  the  high  frequency  curve  should  be 
used  in  the  Inversion  regime. 

Let  us  note,  however,  that  minority  carriers  can  move  quite 
easily  along  the  surface  of  the  semiconductor.  The  surface  mobility 
is  typically  half  of  the  bulk  mobility Hence  we  can  assume  that 
the  minority  carriers  are  in  equilibrium  with  the  tangential  compo- 
nent of  the  electric  field.  This  coupling  between  tangential  field 
and  carriers  is  responsible  for  the  propagation  loss  and  the  so- 
called  longitudinal  convolution  effect  The  former  effect  will 

be  discussed  in  Section  2.2.5. 

Equation  (2.2-8)  allows  us  to  relate  the  bias  voltage  V to 

G 

the  dc  surface  potential  0^^  . Because  of  the  low  capacitance  of 

the  piezoelectric  material  (thickness  y^  = 1 to  5 mm  , and  for  YZ 

LiNbO^  the  permittivity  is  ~ 70  e^)  compared  to  the  airgap 

capacitance  (height  h=0,l  to  0.5  pm,  permittivity  e^)  , all  the 

voltage  is  dropped  across  the  piezoelectric.  If  D is  the  dc  dis- 

sO 

placement  field  in  the  piezoelectric  material,  we  can  write 

\ = °s0  ^p/^y  ^2-2-9) 

and  from  Eq,  (2,2-8),  0^^  is  given  by  the  relation 

Vq  = ^ . F(q0^^kT)  . sgn(0^Q)  . (2.2-10) 
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FIG.  2.2-l--Plot  of  the  nonlinear  coupling  coefficient  B vs  the  dc 
surface  potential. 
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2.2.5  Equivalent  Circuit 


i 


To  relate  the  first  order  term  in  d to  the  value  of  D , 
we  use  an  equivalent  electrical  circuit  in  which  "current"  is  displace- 
ment current  I = joD  and  "potential"  is  electric  potential  0 . The 

admittance  is  then 


Y = joJ)/0 

28 

We  shall  make  extensive  use  of  the  Kino-Reeder  perturbation  theory 
for  the  Rayleigh  wave  amplifier  to  derive  this  circuit. 

Flat  band  conditions. 

Following  Kino-Reeder,  we  compute  the  admittance  y”  at  the  sur- 
face of  the  semiconductor  in  two  different  ways;  1)  looking  at  the 
piezoelectric  through  the  airgap  thus  yielding  the  "piezoelectric 
coupling  equation",  2)  looking  into  the  semiconductor  assumed  semi- 
infinite,  this  yields  the  "electronic  equation". 

To  derive  the  piezoelectric  coupling  equation,  we  write  the  total 
rf  potential  as  the  sum  of  the  acoustic  potential  0^  and  the  space 
charge  potential  0„ 

+ 0S  • (2.2-11) 

For  y < h ; 0g  can  be  computed  by  replacing  the  volume  charge  in  the 
semiconductor  by  a fictitious  charge  p'  at  y = h . Solving 
Poisson's  equation  in  the  layered  system  so  formed  and  matching  the 
boundary  conditions  at  y = 0 and  y = h , it  is  shown  in  Eq . I8  of 
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reference  28  that 


^2fg(h)  = (p'/^^q)  M(ph) 


(2.2-12) 


where 


M(ph) 


*^anh(ph) 

^^0  tanh(3h)  j 


(2.2-13) 


3 is  the  propagation  constant  of  the  wave,  e = fe  p - p^  -j- 

p yy  zz  yz'^ 

the  effective  permittivity  of  the  piezoelectric  material  in  the  y 
direction . 

In  the  normal  mode  theory  of  Auld-Kino,^^  the  acoustic  potential 
0^  Is  the  term  associated  with  the  excitation  of  the  principal  acous- 
tic mode,  the  Rayleigh  surface  wave.  It  is  found  to  be  given  by  the 


relation 


^gpgWZa(3h) 

3'^  -3^ 


(2.2-14) 


where  3'  Is  the  perturbed  propagation  constant,  3 the  unperturbed 
one;  w is  the  width  of  the  acoustic  beam,  Z (3h)  the  acoustic 
impedance  at  the  surface  of  the  semiconductor.  With  the  help  of  the 
boundary  conditions  at  y = h , the  admittance  on  the  lower  side  of 
the  fictitious  charge  is 


Y = J3cd€q  - ja^'/0(h)  . (2.2-15) 
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The  piezoelectric  coupling  equation  follows  from  Eqs.  (2.2-11),  (2.2-12) 
(2.2-lU),  and  (2.2-15): 


where 


M(ph)  3wZ  (3h) 


(2.2-16) 


Yq  = . 

The  electronic  equation  is  obtained  by  writing  the  current,  con- 
tinuity and  Poisson's  equations  in  the  semiconductor.^®  The  general 
solution  is  the  superposition  of  two  waves:  the  Laplacian  wave  which 

corresponds  to  a zero  rf  charge  density  and  has  a variation  exp(-J3z±3y) 
and  the  diffusion  wave  which  propagates  as  exp(-j3z±7y)  . Assuming 
the  semiccnductor  is  semi-infinite,  it  can  be  shown  that  the  admit- 
tance at  the  surface  y = h is 


cm)  €3(7-3) 

Y = Y + — 

7U)  - jo)^3 


(2.2-17) 


where  Y^  - jPcoe  is  the  admittance  of  the  free  semiconductor, 

“c  ^ dielectric  relaxation  frequency  and  p the  mobility 

of  the  majority  carriers.  We  make  the  further  approximation  that 

1 /2 

(%%)  ' or,  written  in  another  way,  3 « l/\^  , where 

o 

0^  = v7(kT/q)p  is  the  diffusion  frequency.  This  restricts  the 

validity  range  of  the  present  theory  to  approximately  f < 50O  MHz  and 

111 

^d  - • Yhis  being  the  range  of  practical  interest,  the 
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assumption  is  sound  and  it  will  make  the  study  simpler  and  tae  equiva- 
lent circuit  easier  to  picture.  The  attenuation  constant  y or  the 
e-folding  distance  of  the  rf  charge  density,  reduces  to 

^ ^ > (2.2-18) 

where  is  the  Debye  length.  The  electronic  equation  becomes 


Y - Y 


(2.2-19) 


+ J(“c/%) 

With  the  use  of  Eqs.  (2.?-l7)  and  (2.2-lg),  „e  can  no»  draw  the  equiva- 
lent circuit  shown  In  Fig.  2.2-2  where  we  Ignore,  for  the  foment,  the 

two  cotnporion. t s Y,  ^nd  V tua 

inversion  ^depletion  ‘ acoustic  proper- 

ties  of  the  piezoelectric  substrate  are  modeled  by  an  inductance  £ 

across  which  the  potential  developed  is  0^  . The  gap  coupling  is 

represented  by  the  capacitance  and  the  bulk  of  the  semiconductor 

contains  a lossy  term  and  a capacitive  term  . The  respective 

expressions  are  given  below: 


1 pwZ  (ph) 

= ^ 

° a.3'^-g2 

=gap“  esg/ilteh) 

■S,  = 


(2.2-20) 
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FIG.  2.2-2--Equivalent  circuit  used  in  the  evaluation  of  linear 
coupling,  and  loss. 


It  will  be  noted  that  thanks  to  the  assumption  o)  « (o)  o^)  we 
can  ignore  the  effect  of  the  terms  and  which  should  appear  as 

capacitors  in  parallel  with  bulk  admittances.  Let  us  also  note  that 

for  the  semiconductors  of  interest  and  hence  the  resistive 

term  is  small, 

R;,  « l/(C^t»)  , 

thus  inferring  that  the  propagation  loss  due  to  the  presence  of  the 
semiconductor  is  small. 

The  "linear  gap  coupling"  term  D^/0^  is  now  straightforward  to 
evaluate 


= 


X 

1 j 


M(ph)  [ eM(3h) 


cu^eM(ph) 


-1 


(2.2-21) 


The  last  bracket  in  Eq.  (2.2-21)  is  equal  to  unity  to  within  10^  or 
less  for  10  to  lOO  fl-cm  silicon. 


Accumulation  conditions 

The  majority  carrier  concentration  is  larger  near  the  surface  than 
the  bulk  value  (N^)  . The  density  of  this  accumulation  layer  varies 

essentially  exponentially  with  distance  and  can  be  calculated  by  inte- 
gration of  Eq.  (2.2-2);  a useful  measure  of  the  extent  of  the  accumula- 
tion layer  is  the  distance, from. the  surface,of  the  center  of  mass  of  the 
excess  carrier  distribution.  This  distance,  called  , is  found^^  to 

be  equal  to  the  bulk  Debye  length  X^  for  weak  accumulation  and  is  a 
slowly  decreasing  function  of  the  dc  surface  potential.  Experimentally, 


"r-T  r 
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the  stronpett  accumulation  we  have  been  able  to  achieve  corresponds 

to  q0gQAl  = 5 , when  has  decreased  to  O.5  . The  total 

29 

excess  charge  can  also  be  computed,  ^ and  we  can  thus  define  an 
equivalent  uniform  surface  concentration  N ^ as  follows: 


accn 


= N , + m/21. 

accn  d ' c 


(2.2-22) 


where  ZiN  is  the  total  excess  charge  per  unit  area.  The  gap  coupling 
is  still  given  by  Eq.  (2.2-21)  where  X^  is  now  computed  for  a concen- 

Only  the  correction  term  is  affected  and  it  becomes 


tration  N 


accn 
even  smaller. 


Depletion  conditions. 

We  use  the  simple  depletion  layer  approximation.  This  introduces 

an  extra  capacitance  in  series  with  C where 

. Sep 

^depletion  “ ' (2.2-25) 

1 /2 

I = (2e^gQ/qNj)  ' is  the  depletion  length.  In  other  words,  the  air- 
gap  h has  now  increased  to  h + and  all  previous  results  are 

directly  applicable. 

Inversion . 

The  depletion  width  continually  increases  until  inversion  sets  in. 
This  occurs  when  becomes  approximately  twice  the  bulk  potential 

0g  . For  larger  depletion  width  will  level  off  and  any 

increase  of  surface  potential  will  correspond  to  an  inversion  of  the 
surface.  The  maximum  depletion  width  is  given  by  the  expression 

'n.,K  = ■ 
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For  10  to  100  n-cm  silicon,  is  10  to  11,  hence  the  mavimum 

depletion  width  is  approximately  6.5  . The  inversion  extends  into 

the  semiconductor.  To  keep  the  calculations  simple  but  nevertheless 
account  for  the  charge  spreading,  we  model  the  inversion  layer  as  an 
infinitesimally  thin  charge  sheet  located  at  its  true  center  of  mass, 
a distance  away  from  the  surface.  This  decreases  the  effective 

depletion  length  by  and  increases  the  airgap  by  L^^eQ/'e  . it 

is  necessary  to  take  into  account  the  spreading  of  the  inversion  charg 
because  we  are  concerned  with  weak  Inversion  and  with  calculating  the 
capacity  of  the  depletion  layer  as  accurately  as  possible.  The  effect 

of  the  finite  thickness  of  the  inversion  layer  L^.  is  not  negligible 
Since  for  weak  inverslon29 

Sion  layer  that  we  have  generated  experimentally,  L . has  a value 

c i 

which  is  still  of  the  order  of  > 1 ^ 

^ . The  total  inversion  charge  per 

unit  area  is  given  by  the  relation 


- VdS' 


(2.2-211 


This  relation  deflnea  the  parameter  g+  . And  aa  in  the  case  of 
accnmnlatlon,  g*  and  are  fnnctiona  of  and  they  have -bet 

plotted  by  Many,  et  al .=9  gaf erring  again  to  Kino-Reeder,^^  rhe  admi 
tance  of  a sheet  of  charge  is 


1 - ja./o^. 


(2.2-2 


where  and  are  the  dielectric  relaxation  and  diffusion  fr( 

quenolea,  reapeotively,  for  minority  oarrlers.  b la  the  width  of  tt 


equivalent  charge  aheet.  From  Eqs.  (2,2-21,)  and  (2.2-25).  we  obtain 
the  admittance  of  the  inversion  layer  as 


Y. 

inv 


(2.S-26) 


1 - Jo,/n^^ 

and  we  can  now  draw  the  full  equivalent  circuit.  It  is  ahown  in  Fig.  2.2-2. 
To  iirat  order,  la  .ainly  capacitive  and 

tive  introducing  aome  extra  loaa  and  hence  reducing  the  current  I flowing 
into  the  aemlconductor,  therefore  reducing  the  convolution  output.  Thla 
reduction  occura  vainly  in  the  inveralon  regime;  it  ia  therefore  of  Inter- 
eat  to  evaluate  the  ratio  of  I m the  preaence  of  an  Inveralon  (charac- 
terizod  by  g ) to  I at  maximum  depletion  (g"*^  = o)  ; 

12 

(2.2-27) 


5?  = 


2 

_i(g''=o) 

[0a/l>d(8*'O). 

The  ..reduction  factor..  * la  atr.ightforward  to  evaluate  ualng  the  equiva- 
lent circuit,  and  the  plota  for  and  g*  , it  la  found  that  B la 

a monotonlcally  decreaalng  function  of  g-^  and  for  gl  ~ lo  , Mch 
corresponds  to  a medium  inversion,  » is  of  the  order  of  lo'^.  Plots 

el  R are  given  in  Fig.  2.2-5  in  the  two  particular  cases  aummarlted  in 
Table  2.1. 


TABLE  2.1 


Piezoelectric 

YZ  LiNbO^  g ^ ^ 

Semiconductor 

Silicon  length  = 2 cm 

B = 7.10“m-5  c/cg  = 11.7 

U.J  = 0.025m^/Vs  e 0.07m^/Vs 

n-type 

p-type 

Airgap 

1500  S 

3000  X 
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REDUCTION  FACTORS 


FIG.  2 .2-3“"Reduction  factor  9 vs  inversion  strength  for  the 

two  cases  given  in  Table  2.1. 
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Propagat-ion  loss  due  to  the  semiconductor . 

The  dispersion  relation  for  the  semiconductor-airgap-piezoelectric 

system  is  obtained  by  writing  the  condition  for  total  impedance  in  the 

equivaler:,  circuit  to  be  zero.  In  the  expression  for  i.  it  is  con- 

a 

venient  to  assume  that  the  perturbed  acoustic  wave  has  a propagation 
constant  of  the  form 

jp  = J3  + (2.2-28) 

where  and  are,  respectively,  the  incremental  loss  and  phase 

delay  constants  for  the  system.  Equating  real  and  imaginary  parts  to 
zero  in  the  dispersion  r-lation  yields  the  propagation  loss  a_ 

(in  Nepers  per  meter). 

2.2.U  Acoustic  Impedance 

The  acoustic  potential  at  the  surface  of  the  piezoelectric  is 

related  to  the  acoustic  power  through  the  acoustic  impedance  as 

follows : 


0_(O)  = 2P  Z (0) 
3 a a J 


(2.2-29) 


is  given  by  the  perturbation  formula  Eq . (85)  in  reference  28: 


Z(O)  zz 


wa>(ep  + Co) 


Av(0) 


(2.2-50) 


|Av(0)/v|  IS  the  relative  surface  wave  -elocity  change  when  a shorting 
plana  is  introduced  at  y = 0 . Since  the  rf  electric  field  associated 

»lth  the  surface  wave  must  vary  as  exp(-py)  in  the  region  y>0  ,it 


- 42  - 


follows  that 


Z (3h)  - Z (0)  exp(-2ph)  . (2.2-51) 

a 3 

2.5  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
2.5.1  Flat  Band  Conditions 


We  first  examine  the  flat  band  results  since  they  will  provide 
simple  design  curves  to  select  the  airgap  thickness  and  the  conduc- 
tivity of  the  semiconductor. 

The  convolution  efficiency  can  be  normalized,  thus  defining  the 
.y// -value, as  - Vqw|2P^^^P^^^)  . For  flat  band  conditions, 

the  full  exp?;ession  is  derived  from  Eqs.  (2.2-1),  (2.2-7),  (2.2-21), 
(2.2-50),  and  (2.2-51).  We  obtain  for  , 


.U  = 


Av(0) 


m exp ( -2ph ) 
v'  (ph ) 


1 + + j 

eM(3h) 


0)€, 


U)^eM(ph) 


-2 


(2.5.1) 


and  for  the  attenuation  constant  a 


a 

r 


HqN^ 

Av 

V 

V 

exp(-2ph) 


+ e 

0 p 


1 + 


(V5 


Cl)  e 

+ M(3h)  — 
ctx;„ 


-1 


. (2.5-2) 


Recalling  the  assumption  o)  « (0)^0^)^/^  we  can  greatly  simplify  the 
expression  for  and  write  it  in  the  form 


2 

0)  Av 


a 


nqN^  V 6(6p  + eQ) 


M’^Oh) 


(2-5-5) 
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The  loss  expiTessed  in  dB  in  a length  L is 

Lg  = 8.6  a^L  . (2.5-1+) 

The  loss  is  seen  to  be  proportional  to  the  resistivity  of  the  semicon- 
ductor. As  mentioned  earlier  in  Section  2.2.5,  the  last  bracket  in 
the  expression  for  , //  is  a correction  term  and,  in  first  approximation, 
we  may  neglect  it.  So  with  the  help  of  Eq.  (2 .5-5),, 7/  can  be  written 
simply  as  follows: 

++ 

_ (2.5-5) 

5 

If  we  now  include  the  prop  Ration  loss  term,  we  obtain  the  open  voltage 
V 

op 

“-//exp  (-ar //2)  =^/7lo"^s/^°  . (2.5-6) 

When  the  approximations  leading  to  Eqs.  (2.5-5)  and  (2.5-5)  are  valid, 
we  obtain  the  following  important  result,  first  given  by  Otto:  For  a 

given  semiconductor  and  piezoelectric  material,  the  convolution  is 
optimum  when  the  propagation  loss  introduced  by  the  semiconductor  is 
1 Neper  or  8.6  dB . In  Fig.  2.5-l(a)  and  (b),  we  show  plots  of  the 
loss  Lg  and  the  ..//-values  t//  and  10  S^^  vs  conductivity  for 

airgaps  ranging  from  0.1  +jm  to  1.0  (ira.  The  expressions  given  by 
Eqs.  (2.5-I)  and  (2.5-2)  are  used.  Both  n-  and  p-type  silicons  are 
considered  and  the  operating  frequency  is  100  MHz.  The  silicon  length 
is  L = 2 cm  and  the  piezoelectric  material  is  YZ  LiNbO  . 
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SILICON-LiNbO; 

^=0.25m^/Vs 
f=  100  MHz 
P-TYPE 


PROP.  LOSS 
INCLUDED 


CONDUCTIVITY  (mho/m) 


Propagation  loss  and  //-values  vs  conductivity  for  various 
airgaps . 

(b)  p-type  silicon. 


The  approximation  made  earlier  limits  the  validity  range  of  these 
results.  A direct  comparison  with  the  full  theory  of  Otto^^  shows  that 
the  error  is  less  than  10^  if  the  frequency  f ^ the  conductivity  cr 
and  the  mobility  ^ are  such  that 

f < 120.^^  10^ 

for  a LiNbO^ -Silicon  system.  The  present  theory  is  therefore  limited 
to  approximately  f < 50O  MHz  and  tr>0.1  U/m  . In  practice,  inter- 
digital transducers  are  not  easily  made  to  operate  above  1 GHz  and 
silicon  wafers  of  conductivity  below  0.1  U/m  tend  to  be  unstable. 

This  theory  is  therefore  well  adapted  to  the  design  of  practical  con- 
volvers . 

2 -3 -2  Results  for  an  Untreated  Silicon  Surface  - Deep  Depletion 
Model 

The  case  of  untreated  silicon  (n-type  especially)  is  dealt  with 
separately  because  of  its  practical  importance:  it  has  been  used  in 

most  of  the  early  convolvers  and  it  provides  a very  large  nonlinear 
coupling . 

It  is  now  a well  established  fact  that  the  surface  of  a semi- 
conductor, like  silicon,  consists  of  an  oxide  film,  some  10  X thick, 
together  with  various  chemicals  adsorbed  from  the  etchants  and  the 
atmosphere.  This  creates  electronic  states  - the  "surface  states"  - 
whose  energy  levels  lie  within  or  around  the  silicon  bandgap.  When 
the  surface  has  been  given  no  special  treatment,  like  thermal  oxidation. 
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we  may  assuni^  that  the  surface  state  density  is  very  high  - in  the 

,^15  -2 

iO  cm  range  comparable  to  the  value  found  under  a Schottky 
barrier.  We  may  estimate  the  surface  potential  by  using  the  results 
obtained  for  Schottky  barriers  in  which  it  is  normally  assumed  that 
there  is  an  interfacial  oxide  layer  of  atomic  dimension  6 at  the  sur- 
face  of  the  semiconductor.  In  addition^  it  is  assumed  the  presence 
of  surface  states  at  the  interface  between  semiconductor  and  the  layer 
of  thickness  6 . The  Schottky  barrier  treatment  leads  to  the  evaluation 

of  the  surface  potential  at  a free  surface  and  the  surface  states.  It 
is  found  that  the  Fermi  level  is  pinned  down  approximately  O.3  V above 
the  valence  band  edge  and  the  surface  state  density  is  of  the  order  of 
10  cm  . For  an  n-type  silicon  of  10  Q-cm  resistivity,  this  means 
a dc  potential  given  by  the  relation 

q^so/'"'^  = -20  , 

i.e.,  a slight  surface  inversion  (the  parameter  g'*"  defined  in 
Section  2. 2. 3 is  less  than  O.Ol).  Although  the  surface  is  slightly 
inverted,  we  shall  compute  the  nonlinear  coupling  as  if  the  surface 
were  in  deep  depletion  (see  discussion  in  Section  2.2.2).  Figure  2.3-2 
theoretical  plots  of  the  ..//-value  in  the  depletion  and  accumula- 
tion regimes.  has  been  normalized  to  its  value  at  flat  band.  As 
the  depletion  layer  widens,  the  . //-value  only  moderately  increases 
under  the  combined  effect  of  the  increase  of  the  nonlinear  coupling 
and  of  the  decrease  of  the  linear  coupling. 
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RELATIVE  ,^-VALUE 


FIG. 


2.5-2-- 


. //-values  in  the  case  of  depleted  or 
semiconductor  surface. 


accumulated 
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Several  experimental  devices  were  constructed  and  tested  to  check 

i'l 


these  results.  Both  LiNbO^  and  were  used  in  conjunction 

with  various  silicon  wafers.  The  output  power  developed  across  a 
50  n load  was  measured  and  the  convolution  efficiency  computed:  this 

has  been  defined  in  Section  2.1.2  as  , the  ratio  in  dB  of  the  out- 
put untuned  power  ^ (in  50  n)  to  the  product  of  the  input  elec- 
trical powers  and  P^^^  , all  powers  in  nW . To  obtain  the 

open  circuit  output  voltage  and  hence  , the  o itput  electrical 

circuit  must  be  known.  Accurate  measurements  of  the  output  impedance 
of  the  convolver  were  carried  out  i.n  the  range  150~500  MHz  by  using 
an  HP  netTiork  analyzer.  From  these  data^  the  equivalent  circuit 
shown  in  Fig.  2.3-3  was  derived.  The  various  elements  are  easily 


' identified.  C.  is  the  capacitance  between  the  semiconductor  and  the 

H 1 


ground  electrode  over  the  acoustic  beam.  is  the  extra  capacitance 

to  ground,  is  the  resistance  of  the  output  gold  leads,  is  the 

inductance  of  the  output  leads  and  grounding  connections,  and  is 

the  extra  capacitance  to  ground  in  the  lead-connector  interface.  By 
theoretical  fitting  of  the  measured  impedance  curves  vs  frequency,  the 
values  of  4 C^,  and  were  obtained.  The  ratio  Cj^/(Cj^+C2) 

was  determined  from  a purely  geometrical  argument:  it  is  equal  to  the 

ratio  of  the  acoustic  beamwidth  to  the  semiconductor  width.  This  cir- 
cuit applies  very  well  to  convolvers  made  with  LiNbO  . When  Bi.^GeO-. 

5 12  20 

is  used,  the  capacitors  and  have  to  be  replaced  by  parallel 

RLC  resonant  circuits  to  account  for  some  internal  acoustic  resonances 
in  the  acoustic  medium.  Once  the  output  circuit  is  known,  its  "voltage 
transfer  function"  can  be  evaluated.  This  is  defined  as  the  ratio  of 
the  voltage  actually  dropped  across  a 50  0 load  to  the  open  circuit 
output  voltage  . If  | is  this  transfer  function,  the  power 

loss  R in  the  output  circuit  is  given  by  the  relation 


R = 20  log  I 

Some  of  the  experimental  results  iox  are  given  in  Table  2.2  and 
compared  to  the  theoretical  predictions  for  deep  depletion  conditions. 
As  can  be  seen,  the  agreement  is  very  good. 
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TABLE  2.2 


. //^-Values  Theory  and  Experiment 


Semiconductor:  n-type  Silicon;  untreated;  2 cm  long;  p = 0.1  m^/Vs 

Piezoelectric:  YZ  LiNbO^;  w = 1 .25  mm;  h = I500  X;  f = I08  MHz 

(111),  (in)  w = 1.25  mm;  h = 50OO  f = 52 


Piezoelectric 


Si  Resistivity  (Q-cm) 
Concentration  (m’^) 

2-port  Insertion  Loss^^^ 
^TR 

Convolution  Efficiency 

(dB)^^^ 

Output  Circuit  Effi- 
ciency 

at  flat  band  (theory) 

depletion/,.//' flat^^^ 

band 

depletion  (theory) 

(experiment) 

(V-m/W) 


LiNbO, 

3 

®h2‘='°20 

7 

7 X 10^° 

7 

7 X lo^° 

9 

5.5  X 10^° 

7 

7 X 10^° 

14 

22 

20 

18 

-56 

-58 

-66 

-80.5 

0.4 

0.7 

0.2 

0.18 

4.0  X lo"5 

4.0  X io"5 

5 -5  X 10  ^ 

1.0  X 10“^ 

1.50 

1.50 

1.4o 

2.10 

6.0  X 10”^ 
5.6  X 10'^ 

6.0  X 10'^ 
6.3  X 10'^ 

7.7  X lo’^ 

7.0  X 10'^ 

2.1  X 10’^ 
1.7  X 10'^ 

o D “ ^^ausuucer  and  propagation  losses. 

(in  50  n)  to  product  of 

. terminal  powers.  All  powers  in  iriW. 

(3)  atio  of  voltage  actually  developed  across  a 50  Q load  to 
open  circuit  output  voltage. 


2.5«5  Results  for  Oxidized  Silicon  — Off  Flat  Band  Regimes  — Surface 
States 

Although  convolvers  built  with  untreated  silicon  are  among  the  most 
efficient,  the  very  large  surface  state  density  reduces  the  light  sensi- 
tivity of  the  semiconductor  and  it  makes  it  almost  impossible  to  vary 
the  dc  surface  potential  by  applying  any  reasonable  bias.  If  the  sur- 
face of  the  semiconductor  is  carefully  cleaned  and  thermally  oxidized, 
then  annealed  in  an  inert  atmosphere  such  as  nitrogen,  the  surface  state 
concentration  can  be  reduced  by  several  orders  of  magnitude. 

Like  the  bulk  impurity  levels,  a surface  state  can  be  a donor  state 
or  an  acceptor  state.  A donor  state  is  either  neutral  or  positive,  when 
it  donates  an  electron.  An  acceptor  state  is  either  neutral  or  negative, 

when  it  accepts  an  electron.  The  probability  of  occupancy  of  a surface 

22 

state  by  an  electron  is  given  by  Feimi-Dirac  statistics:  a state  will 

be  occupied  if  its  energy  level  lies  below  the  Fermi  level.  It  follows 
that  a donor  state  will  be  positive  (no  electron)  if  its  energy  level 
lies  above  the  Fermi  level;  and  an  acceptor  state  will  be  negative 
(electron  trapped)  if  its  energy  level  lies  below  the  Fermi  level.  In 
all  cases,  as  the  Fermi  level  crosses  a state  level  upward,  the  surface 
becomes  more  negative.  In  the  following,  we  shall  only  be  interested 
in  the  net  surface  charge  and  shall  no  longer  specify  the  type  of  the 
states  . 

By  applying  a bias  voltage  V to  the  structure,  the  energy  bands 
can  be  bent  at  the  surface  relative  to  the  Fermi  level.  The  probability 
of  occupation  of  the  surface  states  will  thus  vary  as  a function  of  the 
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applied  bias;  so  will  the  charge  located  in  these  states,  denoted  by 
Q (V  ) . In  the  case  of  thermally  oxidized  silicon,  it  is  found 

that,  in  the  absence  of  any  bias  (V  = O),  the  electron  concentration 
is  increased  at  the  surface:  an  n-type  silicon  will  therefore  be 

accumulated  and  a p-type  depleted  or  inverted.  This  n-type  surface 
conductance  is  attributed  to  the  presence  of  a net  donor-type  surface 
state  concentration  resulting  in  a positive  charge  Q (V_  =0).  As 
a consequence,  the  theoretical  curves  will  be  shifted  towards  inver- 
sion or  accumulation  by  an  amount 


V = 


oxide 

A further  voltage  shift  will  occur  as  the  bias  is  increased. 

We  use  the  results  derived  in  Part  2.2  to  obtain  a theoretical 

plot  of  the  relative t/d-value  vs  bias  voltage  in  the  accumulation  and 

depletion  regimes.  As  shown  in  Figs.  2.J>-h{a)  and  2.5-5(a)>  we  obtain 

a very  good  fit  with  the  experimental  data  (dots)  if  we  allow  for  the 

presence  of  a charge  Qgg(O)  ^t  the  Si-SiO^  interface.  The  p-type 

semiconductor  starts  out  slightly  depleted  corresponding  to  a Qgg(0)/q 
9 “2 

of  8 X lO^cm  (90  volt  shift).  The  n-type  semiconductor  starts  out 

9 -2 

accumulated  with  a Q (0)/q  of  9 X lO^cm  (100  volt  shift). 

s s 

Figures  2.3“^(a)  and  2.5“5(a)  only  show  relative  values  oii4(  ; abso- 
lute values  are  compared  in  Table  2.3,  now  that  v?e  know  what  the  flat 
band  voltages  are.  The  agreement  is  good  here  too. 

As  the  bias  forces  the  surface  of  the  semiconductor  to  invert, 
the  Fermi  level  begins  to  scan  through  the  surface  states  energy  dis- 
tribution near  one  edge  of  the  bandgap;  this  will  introduce  a further 


shift  of  the  characteristics.  In  Figs.  2.5-4(a)  and  2.3-5(a)  the  dashed 
curves  are  drawn  in  the  absence  of  any  more  surface  states;  to  account 
for  them  these  curves  must  be  shifted  by  an  amount  AV  shown  in 
Figs.  2.3"^(b)  and  2.3“5(b)  as  functions  of  the  surface  potential 
q0so/^^  • The  surface  state  distribution  per  unit  energy  N is 

then  obtained  by  graphical  differentiation 

1C., 

^ss  ^ ~ ‘^^ss/^^^sO  ^ ^ ‘^(^V)/d(q0^Q/kX)  states/cm  /eV  . 

As  seen  in  Figs.  2.3“^(b)  and  2 , this  distribution  exhibits 

maxima  near  the  band  edges.  The  total  number  of  surface  states  measured 

10  -2 

in  this  way  is  Q = 2.5  X 10  cm  near  the  conduction  edge  and 
10  -2 

Q = 5-2  X 10  cm  near  the  valence  band  edge.  Since  it  is  com- 
s s 

paral'ly  harder  to  accumulate  strongly  than  to  invert  a semiconductor, 
the  surface  state  distribution  could  only  be  studied  at  one  extremity 
of  the.  bandgap  with  a given  semiconductor.  The  surface  state  concen- 
tration is  rather  small  because  of  the  careful  cleaning  and  oxidation 
procedures  followed  (see  Appendix  B). 

The  measure  of  the  t--// -value  of  the  semiconductor  convolver  there- 
fore provides  a new  method  of  measurement  of  the  surface  state  density. 
This  density  is  plotted  vs  energy  in  Fig.  2.3-6  and  compared  to  the 
results  obtained  by  other  authors These  authors  used  MOS 
diodes  and  measured  either  their  capacitance  or  their  conductance  as 
a function  of  the  applied  gate  voltage  or  by  varying  the  temperature. 

The  distribution  we  obtain  is  more  sharply  peaked  around  two  definite 
energy  levels,  0.2  eV  and  O.83  eV.  The  surface  state  density  was  not 
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FIG.  2.3-4(a) 


.^-value  and  propagation  loss  vs  bias  voltage  (n-tvne 
near  the  valence  ban'rcLr 


0.25  0.24  0.23  0.22  0.21  0.20  0.19 
ENERGY  FROM  VALENCE  BAND  EDGE  (eV) 


FIG.  2.3.4(b)-Voltage  shift  and  surface  state  distribution  (n-tyn 
silicon). 
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FIG.  2. 3. 5(a)--. .//-value  and  propagation  loss  vs  bias  voltage  (p-type 
silicon).  The  dots  are  the  experimental  points.  The 
full  line  curve  corresponds  to  the  complete  theory; 
the  dashed  line  curve  to  the  theory  in  the  absence  of 
surface  states  near  the  conduction  band  edge. 
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FIG.  2.3-5(b)--Voltage  shift  and  surface  state  distribution  (p-typ 
silicon). 


TABLE  2.3 

•^-Values  Theory  and  Experiment 


Semiconductor:  Silicon;  1000  X Si0„  layer;  2 cm  long  . 5 «■  A 

^ . .07 

Piezoelectric:  YZ-LiNb03;  w = I.25  mn;  f = I08  MHz 


Silicon  type 

.X  (*) 

concentration  (m  ■^) 
resistivity  (n-cm) 

N 

7 X 10^° 
~ 12 

P 

3.5  X 10^° 
~ 70 

Airgap  (X) 

1500 

3000 

2-porc  Insertion  Loss^^^  (Ho  bins) 
(■») 

16 

50 

Convolution  Efficiency (No  bias) 
Ft  (dB) 

-66 

-7it 

Output  Circuit  Efficiency^^^  5 

0.25 

0.20 

e/d<-value  (Experiment,  No  bias) 

2.8  X lo"^ 

8.8  X 10‘5 

^ at  flat  band/M  at  OV  Bias^^^ 

1.7 

0.6 

e^-value  at  flat  band  (Expt) 

it  .8  X 10*^ 

5.3  X 10"5 

n^-value  at  flat  band  (Theory) 

(V-mA) 

it. 5 X lo"^ 

5 X lo’^ 

n/df-value  at  optimum  voltage 

9. it  X lo'5 

8.8  X 10*^ 

(*)  from  MOS  C-V  Measureraents 

(1)»  (2),  (3)  - See  Table  2.2 

(It)  from  Figs.  2.3-it(a)  and  2.'-'^va' 
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DENSITY  OF  SURFACE  STATES 


ENERGY  (eV) 


FIG.  2.3-6-Surface  state  density  ,£  SI-SIO  structure  as  obtained  from 

(aj  present  work, 

(b ) After  Nicollian  and  Goetzberger^^ 

(c)  After  Castro  and  Deal .52 


measured  very  close  to  the  band  edges  E and  E because  it  is 

C V 

almost  impossible  to  have  the  Fermi  level  get  close  to  these  band  edges 

by  the  use  of  an  applied  gate  voltage  only.  One  needs  to  use  the 

22 

temperature  method.  The  measurement  of  the  voltage  shift  between 
ideal  and  real  curves  leads  to  the  evaluation  of  the  surface  state 
density.  This  measurement  is  most  inaccurate  in  very  weak  inversion, 
i.e.,  when  the  two  curves  begin  to  split.  In  addition,  in  this  same 
range,  the  surface  potential  is  a rapidly  varying  function  of  the  bias 
voltage  and  hence  is  determined  more  inaccurately.  These  two  sources 
of  errors  may  account  for  the  discrepancies  between  curve  (a)  and 
curve  (b)  or  (c)  in  the  energy  range  0.25  to  0.1+  eV  and  O.65  to  0.8  eV . 
In  addition,  it  should  be  pointed  out  that  the  method  used  by  Nicollian 
and  Goetzberger,  in  particular,  and  the  indirect  way  in  which  they 
derived  the  surface  state  density,  is  very  much  argued  upon  at  this 

57 

time;  no  definite  answer  has  yet  been  cast  on  this  matter. 

There  are  two  new  features  about  the  method  we  propose  here.  First 
we  measure  a second-order  characteristic,  namely  the  derivative  of  the 
capacitance  w.r.t.  the  surface  potential.  In  addition,  this  method  is 
non-destructive,  since  no  metal  is  deposited  on  the  silicon  oxide, 
therefore  no  extra  impurity  centers  are  added  to  the  surface.  We  must, 
however,  note  that  a careful  measurement  should  be  made  under  vacuum  so 
that  the  possible  contaminants  from  the  surrounding  air  do  not  affect 
the  long  term  properties  of  the  surface  of  the  semiconductor. 

With  this  knowledge  of  the  voltage  shift  introduced  by  the  surface 
states,  we  can  now  plot  the  theoretical  propagation  loss  due  to  the  semi- 
conductor in  all  regimes  [see  Figs.  2. 5-2+ (a)  and  2.5-5(a)].  For  the 
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inversion  regime,  we  used  the  plot  shown  in  Fig.  2.5~7-  It  shows  the 
theoretical  loss  vs  the  inversion  strength  g^  . We  note  that  the 
loss  is  maximum  when  1 < g^  < 10  , which  correspond  to  the  minority 

carrier  resonance.  The  depletion  layer  approximation  is  not  very 
accurate  in  accounting  for  the  propagation  loss.  The  theoretical  loss 
is  somewhat  higher  than  experimentally  observed  in  inversion.  The 
inversion  layer  uniformity  is  Indeed  a very  sensitive  function  of  the 
uniformity  of  the  airgap  and  we  noticed  that  only  portions  of  the  sur- 
face were  inverted;  hence  the  transmission  loss  was  lower  than  expected. 
Besides,  the  model  we  used  for  the  inversion  layer  is  quite  simple  and 
may  not  be  very  accurate  for  light  inversion  when  the  inversion  layer 
extends  several  Debye  lengths  into  the  bulk.  But  the  general  shape 
and  orders  of  magnitude  can  be  correctly  predicted  and  in  the  case  of 
the  cY'/' -value,  the  agreement  is  excellent. 

We  have  only  discussed,  so  far,  the  effect  of  a dc  biasing  voltage. 
If  a pulsed  bias  voltage  is  now  used  and  if  the  measurement  is  made 
right  after  the  bias  is  applied,  no  charge  will  be  exchanged  with  the 
surface  traps.  And  for  all  purposes  we  may  ignore  them.  No  addi- 
tional voltage  shift  will  take  place  and  large  surface  fields  may  be 
created.  As  a second  consequence  of  the  use  of  a pulsed  bias,  no 
minority  carrier  will  be  generated  at  the  surface  from  the  bulk  of 
the  semiconductor  because  of  the  long  generation-recombination  times 
involved  in  this  process.  One  would  therefore  expect  the  surface  to 
be  driven  in  deep  depletion.  The  .///-value  would  then  remain  at  the 
constant  depletion  level  (see  Fig.  P.J-S)  and  the  loss  would  be  very 
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N-TYPE 


FIG.  J"*  inversion  strength  for  the  two  cases 


low  because  of  the  large  depletion  layer.  Although  some  convolvers 
behaved  in  this  manner,  others  did  not.  The  response  of  one  of  these 
is  represented  by  the  triangular  dot  curves  of  Fig.  2.'^-h(a).  The 
mechanism  believed  to  be  responsible  for  this  behavior  is  the  presence, 
at  the  surface  of  the  semiconductor,  of  sources  of  minority  carriers 
like  p"*"  impurity  centers.  Under  the  high  surface  fields,  these 
centers  can  inject  a large  quantity  of  minority  carriers  which  will 
stay  confined  at  the  surface  because  of  the  biasing  field.  This  is  the 
mechanism  that  takes  place  in  a field-effect-transistor  where  one  of 
the  end  p contacts  (the  "source")  injects  minority  carriers  and 
thus  form  the  inverted  channel  of  the  transistor.  Since  these  car- 
riers are  of  majority  type  for  the  generation  centers,  the  supply  will 
be  very  fast.  Hence,  with  a pulsed  bias  an  inversion  layer  is  formed, 
and  one  would  expect  the  ^^-value  to  decrease  as  predicted  by  the  pre- 
vious general  theory  in  the  absence  of  surface  states- [see  the  dashed 
curve  in  Fig.  2.3-k(a)].  The  results  shown  in  Fig.  2.5-J+(a)  confirm 
this  assertion  well.  In  addition,  the  maximum  of  the  propagation  loss 
is  now  visible  because  the  inversion  now  occurs  for  lower  voltages, 
but  the  value  of  the  maximum  is  smaller  than  expected.  Once  again, 
alrgap  nonuniformities  may  well  be  one  of  the  reasons  for  this  dis- 
crepancy. 

Since  the  dc  bias  response  differs  from  the  pulsed  bias  response, 
we  have  been  interested  in  studying  the  transient  phenomena  involved, 
just  as  is  done  in  a normal  MIS  diode  measurement.  For  that  purpose, 
the  bias  voltage  was  turned  on  at  the  time  t = 0 and  the  output 
convolution  signal  was  displayed  against  time  by  inserting  into  the 
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acoustic  convolver,  two  identical  pulse  trains:  the  pulse  length  was 

U |isec  and  the  repetition  rate  approximately  J>0  kHz.  It  was  found 
that  for  bias  voltages  lower  than  5OO  V,  the  output  signal  is  almost 
constant  in  time  and  it  is  therefore  difficult  to  measure  any  decay 
time.  As  the  bias  is  increased,  the  convolution  output  variation  with 
time  becomes  essentially  exponential.  The  time  taken  for  the  steady 
state  to  be  reached  is  maximum  for  a bias  voltage  = UOO  to  60O  V , 
and  it  is  found  to  be  I5  ms.  This  corresponds  to  a time  constant  of 
approximately  2 ms.  These  values  of  corresponds  respectively  to 

a surface  potential  q^^^/kT  of  -I5  and  -2h;  in  other  words,  they 
correspond  to  a semiconductor  surface  in  deep  depletion  at  t = 0^ 

This  time  constant  is  associated  with  the  trapping  rate  in  the  surface 
states  as  we  saw  earlier  and,  on  this  basis,  the  largest  time  constants 
should  be  expected  when  the  background  total  concentration  of  mobile 
charges  is  minimum,  hence  to  a surface  in  deep  depletion.  The 
presence  of  impurities  consisting  of  shallow  acceptor  centers  will 
perturb  the  trapping  process  wnen  they  start  injecting  minority  car- 
riers at  the  surface;  as  a consequence,  the  trapping  time  becomes 
shorter.  No  matter  whether  we  consider  the  surface  states  as  being 
of  donor  or  acceptor  type,  the  concentration  of  mobile  minority 
carriers  (holes  here)  will  decrease  as  trapping  takes  place. 

To  show  this  experimentally,  we  measured  the  propagation  loss 

versus  time  after  a large  biasing  voltage  had  been  turned  on 

(V„  = 1500  V).  For  t < 0 , the  loss  is  approximately  2 dB  [see 

Fig.  2.1-2(a),  V = 0 V].  At  t = 0 , because  of  the  large  hole 

Cj  + 
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injection  due  to  the  presence  of  shallow  acceptor  centers,  tha  surface 
mobile  carrier  concentration  is  very  high  and  the  propagation  loss 
decreases  to  almost  0 dB  [see  the  dashed  curve  in  Fig.  2.1-2(a)  with 
Vg  = 1500  V] . As  time  evolves,  trapping  takes  place,  the  hole  con- 
centration at  the  surface  decreases  and  the  propagation  loss  is  seen 
to  increase  monotonically  until  it  reaches  the  steady  state  value, 
i.e.,  the  dc  bias  value,  of  30  dB  (see  full  line  curve  in  Fig.  2.1-2(a) 
with  = 1300  V).  If  we  now  look  at  the  convolution  output,  at 
t = 0^  , it  decreases  sharply  by  I5  to  20  dB  due  to  the  large  screen- 

ing effect  introduced  by  the  hole  surface  concentration.  With  time, 
this  concentration  decreases,  the  screening  becomes  weaker  and  the 
convolution  efficiency  increases;  but  the  increase  in  propagation  loss 
will  soon  be  preponderant  and  the  convolution  is  seen  to  go  through  a 
relative  maximum  before  it  decreases  again  to  finally  reach  the  dc  bias 
level,  some  25  dB  below  the  zero  bias  level  [see  Fig.  2.1-2(a)]. 

The  results  we  have  just  summarized,  have  shown  that  very  long 
time  constants  are  associated  with  surface  states  on  silicon  surfaces. 
These  experiments  led  Kayakawa  and  Kino^^  to  consider  a new  type  of 
storage  process  using  the  same  convolver  structure.  In  this  process, 
the  information  is  represented  by  the  amount  of  charge  stored  in  the 
surface  states;  it  can  be  written  in  or  read  out  acoustically  and  long 
storage  times  of  2 ms  or  more  have  been  demonstrated.^^  A careful 
measurement  of  the  storage  time  vs  applied  bias  voltage  was  done  by 
Hayakawa;  this  is  in  fact  a direct  measuremr.vt  of  the  surface  state 
time  constant,  and  it  was  found  to  be  in  good  agreement  with  previous 
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theory  and  measurement  by  Nicollian  and  Goetzberger . 

A detailed  quantitative  study  of  the  storage  is  beyond  the  scope  of 
this  work;  but  it  has  recently  been  the  subject  of  considerable  inter- 

5 5 3^  3^ 

est  ^ ' because  of  the  possibility  of  building  acoustoelectric 

correlators  with  memory. 

2.k  CONCLUDING  REMARKS 

It  has  been  shown  that  for  real  semiconductors,  the  effect  of  sur- 
face states  is  dominant  and  cannot  be  overlooked.  In  particular,  if 
untreated  n-type  silicon  is  used,  deep  depletion  must  be  assumed  at 
the  surface.  Several  experiments  were  conducted  and  confirmed  this 
assumption.  Furthermore,  if  a convolver  is  designed  to  be  optimized 
for  flat  band  conditions,  its  convolution  output  can  still  be  enhanced 
if  the  surface  of  its  semiconductor  is  depleted.  The  untreated  n-type 
silicon  surface  is  naturally  almost  intrinsic  and  it  is  therefore  a 

good  choice  when  only  electrical  signal  processing  applications  are 
the  aim. 

Oxidized  silicons  were  studied  in  detail.  The  propagation  loss 
and  convolution  efficiency  were  compared  to  theoretical  predictions  on 
the  basis  of  our  model.  The  same  conclusions  concerning  the  optimum 
operating  conditions  hold.  In  addition,  this  study  provided  a new 
measurement  technique  c£  the  surface  state  distribution  in  silicon. 
Detailed  information  about  the  surface  state  is  supplied  by  the  study 
of  the  variation  of  the  convolution  efficiency  and  of  the  propagation 
loss  with  respect  to  the  applied  bias  voltage.  In  addition,  the  use 
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of  pulsed  bias  and  the  study  of  the  transient  responses  provide  a direct 
measurement  of  the  time  constant  of  the  surface  states.  We  shall  also 
see  in  Chapter  IV,  that  illumination  of  the  semiconductor  affects  the 
convolver  characteristics  and  some  extra  information,  like  the  life- 
time of  minority  carriers,  can  be  retrieved. 

The  semiconductor-airgap-piezoelectric  convolver  is  therefore  a 
very  powerful  tool  for  the  Investigation  of  semiconductor  surfaces. 

The  effect  of  bias  has  been  treated  in  some  detail  here,  but  the  effect 
of  illumination  could  be  studied  more  deeply  and  so  could  the  effect  of 
temperature . 
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CHAPTER  III 


APPLICATION  OF  THE  SEMICONDUCTOR  CONVOLVER  TO 
ACOUSTIC  PARAMETRIC  SIGNAL  PROCESSING 

In  this  chapter^  we  show  how  the  separated  medium  semiconductor- 
piezoelectric  configuration,  studied  previously,  can  be  used  as  an 

efficient  signal  processor,  and  how  it  compares  with  other  acoustic 
devices . 

Parametric  signal  processing  using  acoustics  has  been  the  subject 
of  considerable  Interest  for  the  last  few  years,  because  of  the  new 
techniques  it  brought  to  this  field.  Among  them  are  real  time  con- 
volution,^ ^ and  Fourier  transform,^®  pulse  compression,^^  and  time 
inversion.  Under  the  label  "parametric"  we  include  all  the  tech- 
niques by  which,  from  two  input  signals,  one  can  generate  a third  one 
at  the  sum  or  difference  frequency.  In  an  acoustic  wave  device,  these 
signals  may  be  traveling  waves  or  standing  waves  with  electric  fields 
associated  with  them  if  the  material  is  piezoelectric.  A parametric 
processor  must  therefore  include  a nonlinear  mechanism  allowing  for 
the  coupling  of  the  waves.  There  are  two  major  such  nonlinear  mecha- 
nisms. The  first  is  the  acoustic  nonlinearity  of  the  acoustic  medium 
itself  on  or  through  which  the  acoustic  waves  propagate.  Quate  and 

Thompson  first  studied  bulk  wave  Interactions,  Svaasand^^  and  Luukkala 
ll0 

and  Kino  studied  surface  wave  interactions  for  this  purpose,  in  the 
second  mechanism,  the  electric  field  associated  with  the  acoustic  wave:, 
interact  with  charge  carriers  in  a semiconductor  located  within  or  near 
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the  piezoelectric  medium;  and  wave-wave  coupling  is  achieved  through 
the  space  charge  nonlinearity.  There  exists  two  classes  of  devices 
built  on  this  principle:  one  uses  a piezoelectric  semiconductor. 

This  was  first  investigated  by  Wang^^  and  Turner,  et  al.,^^  with  CdS, 
then  by  Ludvik  ^ with  GaAs.  In  the  other  class,  described  in  this 
report,  the  semiconductor  is  placed  in  close  proximity  to  a piezo- 
electric ^references  this  work)  or  the  piezoelectric  is 

deposited  on  the  surface  of  a semiconductor  thus  forming  a monolithic 
structure. 

The  piezoelectric  nonlinearity  is  limited  and  the  devices  making 
sole  use  of  it  exhibit  low  efficiencies.  Various  methods  of  beam 

46  4y  48 

guiding  ’ or  of  beamwidth  compressing  have  been  used  with  suc- 
cess to  increase  the  power  density  in  the  wave  and  hence  the  non- 
linear efficiency. 

On  the  other  hand,  the  space  charge  nonlinearity  can  be  very 
strong,  and  efficient  devices  have  been  built,  but  the  corresponding 
technology  is  more  difficult.  Some  of  the  problems  such  as  inter- 
action uniformity  in  non-mono lithic  devices  and  material  preparation 
in  monolithic  devices  require  much  care.  One  unique  feature  of  the 
semiconductor  systems  is  their  light  sensitivity  (see  Chapter  IV). 

In  Section  5.I,  we  shall  study  the  semi conductor-airgap- 
piezoelectric  convolver  from  a device  point  of  view.  We  shall  be 
mainly  interested  in  its  electronic  characteristics  such  as  time- 
bandwidth  product,  terminal  efficiency  and  dynamic  range.  We  shall 
compare  the  different  types  of  convolvers  presented  earlier. 
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In  Section  5. 2, we  shall  describe  other  signal  processing  applica- 
tions, in  particular,  code  recognition,  time  inversion,  correlation 
and  Fourier  transformation.  Here,  the  parametric  processor  may  be 
considered  as  a programmable  filter:  one  of  the  Inputs  acts  as  a 

reference  or  a filter  for  the  other.  Since  this  reference  can  be 
varied,  the  filter  is  programmable. 

3.1  CONVOLUTION 

The  standard  semiconductor-piezoelectric  configuration  is  illus- 
trated in  Fig.  3.1-1.  A surface  wave  delay  line  is  polished  optically 
flat  and  Interdlgltal  transducers  are  deposited  at  either  end.  Two 
materials  have  been  used:  Y-cut, Z-propagatlng  LiNbO^  because  of  its 

high  piezoelectric  coupling  and  (lll)-cut,  (111 )-propagating 
(BGO)  because  of  its  slow  acoustic  velocity.  The  standard  transducers 
have  four  finger-pairs,  the  aperture  is  1.25  mm  and  the  linger  width 
and  spacing  is  tt  pm.  On  LiNbO^  (v  = 3^80  m/s)  the  operating  frequency 
range  is  I08  ±11  MHz;  on  BGO  (v  = 1620  m/s)  it  is  52.5  ± 5.5  MHz. 

A silicon  chip  is  pressed  down  on  spacer  rails  deposited  on  the  piezo- 
electric surface.  The  rail  thickness  varies  trom  1000  X to  50OO  8 
and  the  system  relies  on  the  flexibility  of  the  semiconductor  to  follow 
the  rails  to  maintain  a uniform  airgap  and  hence  a uniform  interaction 
strength.  On  LlNbO^,  the  rails  are  made  of  sputtered  SiO;  because  of 
the  chemical  activity  of  BGO,  three  layer  rails  had  to  be  used  on  BGO 
(1000  X of  ThOF;  700  of  ZnS ; I500  X of  ThOF).  Electrodes  on  the  top 
and  bottom  of  the  structure  allow  for  a capacitive  pick-up  of  the 
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FIG.  3 •l“l““Schematic  of  the  semiconductor  convolver 


parametric  ojtput.  The  silicon  chips  are  typically  20  to  50  mm  long, 

5 mm  wide  and  25O  pm  thick.  When  the  semiconductor  is  untreated,  it 
is  only  cleaned  in  the  usual  solvants  and  kept  dust-free.  If  a fur- 
ther oxidation  is  needed,  a special  cleaning  procedure  is  followed  to 
remove  all  contaminants;  after  thermal  oxidation,  the  surface  is  annealed 
in  nitrogen  to  lower  the  surface  state  density  down  to  the  10^*^  cm"^ 
range.  Details  of  this  procedure  are  given  in  Appendix  A. 

Characteristics  of  some  of  the  convolvers  built  on  this  model  are 
given  in  Table  3.I  and  we  shall  now  describe  them  in  detail. 

From  a device  point  of  view,  an  important  parameter  is  the  termi- 
nal efficiency  ; it  is  defined  as  the  ratio  in  dB,  of  the  output 

electrical  power  to  the  product  of  the  input  electrical  powers  when  all 
powers  are  expressed  in  itW . Recalling  the  results  of  Section  2.^  F 

* } ip 

is  of  the  following  form: 

Fx  = 20  log  j I - Transducer  Loss  (dB)  - Output  Circuit  Loss  (dB) 

\ 2 w / 

(5.1-1) 

where  R^^  is  the  load  resistance. 

In  designing  a convolver,  the  choice  of  the  materials,  airgap  and 
resistivity  is  made  from  the  results  of  Section  2.5.  In  our  case, 

LiNbO^  and  BGO  were  used  to  be  able  to  operate  in  the  VHF  range  with 
efficient  transducers;  we  decided  upon  silicon  as  the  semiconductor 
because  of  its  availability  and  because  of  the  good  control  over  its 
resistivity  up  to  100  CJ-cm.  Optimum  choices  of  the  airgap  height  a’-e 
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TABLE  5.1 

Device  Characteristics 


Piezoelectric;  LiNbO^  YZ  or  (111) (111);  w = 1.25 

Semiconductor:  Silicon,  2 cm  long 


Piezoelectric 

LiNbO, 

BGO 

Airgap  (X) 

(1) 

(2) 

(5) 

(i+) 

1500 

1500 

5000 

3000 

Silicon  type 

N 

N 

p 

N 

resistivity  (Q-cm) 

7 

10 

70 

10 

2-port  insertion  loss 
(dB) 


Convolution  Untuned 

Efficiency 
(dB)  Tuned 


^-Value  (V-m/W) 


-Kk.^  -44  -61  -74 

5.6x10'^  9.4  X 10'^  8.8  X 10'^  2.1x10"^ 


Interaction  Uniformity 
(dB) 


Frequency  (MHz) 

Input  Bandwidth  (MHz) 
Processing  Time  (|js) 

Remarks 


± 1 ♦ 5 
108 


± 1.5 

108 

25 

5.5 


± 1.5 

108 

20 

5.5 


Optimum  Bias  Used 


(1)  Same  device  as  in  column  1 of  Table  2.2 

(2)  Sa.ae  device  as  in  column  1 of  Table  2.5  (oxidized  silicon) 

(3)  Sa  le  device  as  in  column  2 of  Table  2.5  (oxidized  silicon) 

(4)  Same  device  as  in  column  4 of  Table  2.2 


I 


m 


then  1000  X to  I500  X with  10  them  silicon  and  30OO  X »lth  100  n-cm. 
If  the  sen.lcondnctor  Is  oxidized,  a bias  voltage  will,  l„  general,  be 

necessary  to  drive  the  sorface  Into  depletion  and  hence  optimize  the 
I'^-value . 

As  seen  from  Eq.  (3.1-1)  the  acoustic  beamwldth  should  be  as 

small  as  possible;  beam  guiding  and  compressing  techniques  have  been 

used  by  other  authors'^.'^.’*'?  for  that  purpose.  In  the  absence  of 

any  such  guiding,  the  minimum  width  Is  limited  by  diffraction  of  the 

acoustic  wave,  and  the  choice  of  a reasonable  Impedance  for  the  Input 
transducers . 

The  output  electrical  circuit  must  be  carefully  designed.  In 
particular,  the  width  of  the  output  electrode  must  be  chosen  as  close 
as  possible  to  the  beamwldth  In  order  to  reduce  any  stray  capacitance 
to  ground  which  would  act  as  a capacitive  bridge  and  reduce  the  maximum 
available  voltage,  simple  spacer  rails  do  not  allow  for  a complete 
eancellatlon  of  the  stray  capacitance.  Output  tuning  will  also 
enhance  the  power  transfer. 

we  have  been  able  to  build  devices  with  terminal  efficiencies  a. 
high  as  . -l(ll  dB  . This  Is  among  the  highest  device  efficiencies 
reported  to-date  If  we  exclude  the  multistrip  coupler  convolver  of 
Shreve.  9 Probably  more  Important  than  the  efficiency.  Is  the  dynamic 
range  of  the  system.  It  Is  limited  by  the  saturation  effects  and  the 
spurious  signal  level.  In  our  devices,  the  3 dB  compression  of  the 
bilinearity  corresponds  to  an  Input  powe-  of  20  dBm;  the  noise  in  the 
detection  system  Is  at  -95  dBm,  so  we  could  expect  a dyraml-  range  of 
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more  than  85  dB . This  is  actually  the  case  when  the  two  input  powers 
are  identical  but,  in  practice,  one  of  the  input  signals  is  taken  as  a 
reference  of  constant  power.  The  dynamic  range  will  then  not  be  limited 
by  the  noise  level  but  by  the  reflections  of  the  signal  of  higher  power 
off  the  end  transducers  and  the  edges  of  the  silicon.  With  the  present 
devices,  the  dynamic  range  is  found  to  be  limited  to  approximately  50  dB . 
Many  solutions  have  been  proposed  to  reduce  the  echoes  on  a surface  wave 
delay  line  and  it  seems  certain  that  with  their  help,  we  could  increase 
the  dynamic  range  to  some  60  to  70  dB . 

A second  parameter  of  importance  is  the  uniformity  of  the  airgap. 
This  will  relate  to  the  accuracy  of  the  convolution  product  and  the 
quality  of  the  system  as  a signal  processor.  Testing  of  the  gap  uni- 
formity, i.e.,  of  the  nonlinear  coupling  uniformity  is  performed  by 
convolving  a narrow  rf  pulse  with  a very  long  one;  the  output  should 
be  trapezoidal  in  shape  with  a flat  top  over  a time  corresponding  to 
the  acoustic  delay  under  the  semiconductor  less  the  narrow  pulse  length. 
Furthermore,  the  shorter  the  narrow  pulse,  the  more  accurate  the  picture 
of  the  gap  uniformity.  The  measure  of  the  amplitude  variation  in  the 
output  waveform  will  define  the  "uniformity"  of  the  gap.  It  has  been 
determined  from  experience  (code  recognition  experiments)  that  a uni- 
formity of  ± 1 dB  is  adequate  for  most  processing  applications. 

With  careful  cleaning  of  the  surfaces  of  the  semiconductor  and  of 
the  piezoelectric,  we  have  obtained  consistently  uniformities  of  ± I.5  dB 
or  less  (see  Table  5.I).  Excellent  results  (±  0.7  dB  over  11  ;isec) 
were  obtained  with  a BGO  delay  line.  This  is  attributed  to  the  softness 
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of  the  mateiial,  whose  surface  can  deform  more  easily  t m match  the 
shape  of  the  semiconductor  surface.  Oscillograms  showing  the  uni- 
formity and  the  accuracy  of  the  convolution  product  in  two  particular 
instances  are  given  in  Fig.  ^.1-2  for  the  BGO-silicon  system. 

The  airgap  height  is  subject  to  some  changes  with  time  and 

atmospheric  conditions.  One  method  to  improve  its  uniformity  and 

50 

stability  was  devised  by  Smith,  et  al.  They  ion-beam  etch  small 
posts  in  the  LiNbO^  and  pressed  a narrow  strip  of  silicon  on  them. 

This  made  the  interaction  strength  more  uniform  and  furthermore, 
solved  the  stray  capacitance  problem  we  mentioned  earlier. 

The  quantity  of  information  a system  is  able  to  process  is 
measured  by  its  Time-Bandwidth  product,  TB,  where  T is  the  maximum 
duration  of  the  input  signals  which  can  be  accurately  processed  and 
B is  the  3 dB  bandwidth  of  the  system,  as  measured  on  the  input. 

In  the  case  of  the  acoustic  surface  wave  convolvers,  T depends 
on  the  interaction  length  and  on  the  acoustic  velocity  of  the  piezo- 
electric. Technological  difficulties  limit  the  usable  length  of  the 
semiconductor  to  20  to  JO  inra;  these  difficulties  are  related  to  the 
availability  of  large  good  quality  semiconductor  wafers  and,  above 
all,  to  the  ability  to  maintain  a uniform  gap  over  long  distances. 

Let  us  point  out  that  similar  limitations  exist  with  the  integrated 
monolithic  devices.  On  the  other  hand,  this  problem  is  not  as  acute 
in  the  devices  using  only  the  piezoelectric  nonlinearity.  Low  veloc- 
ity piezoelectrics  are  suitable  to  increase  T but  they  will  make  the 
realization  of  high  frequency  interdigital  transducers  harder,  because 
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FIG.  5 . 1 “2--Bi  j^^GeO^^-airgap-Silicon  convolver. 

(a)  Autoconvolution  of  one  square  pulse; 

(b ) Autoconvolution  of  a double  pulse; 

(c)  Picture  of  the  uniformity  of  the  interaction  over  2 cm 


I 
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the  acoustic  wavelength  and,  hence,  the  finger  width  will  be  smaller. 

As  can  be  seen  from  the  derivations  of  Chapter  II,  the  nonlinear 
effect  is  a slow  function  of  frequency;  the  bandwidth  will  therefore 
essentially  be  limited  by  the  transducer  response.  Wideband  opera- 
tion can  be  achieved  with  transducers  with  small  number  of  finger 
pairs  and  adequate  tuning  circuits,  but  the  efficiency  will  be 
reduced.  In  our  devices,  we  used  finger  pair  transducers;  on  LiNbO 
they  are  tuned  by  means  of  a simple  series  inductov  and  they  yield  a 
20$  bandwidth  togethei  with  a low  Insertion  loss;  on  BGO  they  are  not 
quite  optimal  and  the  bandwidth  is  only  15^. 


Devices  built  with  k finger  pair  transducers,  and  employing  2 cm 
long  silicon  pieces,  exhibit  typical  TB  products  in  the  100  to  200 
range,  depending  on  the  piezoelectric  material  and  the  operating 
frequency  (50  to  100  MHz  on  BGO;  100  to  200  MHz  on  LiNbO  ). 

For  comparison  purposes,  we  have  gathered  in  Table  5.2,  the  char- 
acteristics of  most  of  the  acoustic  surface  wave  convolvers  reported 
on  to-date,  with  the  exception  of  those  corresponding  to  the  present 
work;  these  were  detailed  in  Table  5.1.  The  device  parameters  we  have 
been  discussing  in  this  section  are  given  when  available;  the  .^-values 
are  also  included.  The  most  efficient  structure  makes  use  of  a multi- 
strip coupler  between  the  acoustic  beam  and  the  semiconductor .^9 
Although  the  efficiency  is  very  high  (F^  up  to  -25  dB),  for  this 
reason  the  system  saturates  for  very  low  levels  (of  the  order  of  0 dBm); 
furthermore,  it  is  very  sensitive  to  reflections  at  the  strips.  This 
limits  its  applications  as  a practical  convolver.  The  semiconductor- 


*)  Type  I:  Plesoelectrlc  nonlinearity  only 


airgap-piezoelectric  systems  have  good  dynamic  ranges  and  the  TB  pro- 


duct can  be  as  high  as  50O  to  1000  when  one  uses  special  technologies 


to  maintain  a long  uniform  »ap,  and  wide  band  transducers.^^  Finally, 


very  practical  devices  usin^;  LiNbO^  and  beam-compressors  have  been 


40 

demonstrated  and  they  show  good  efficiency  and  dynamic  range  and  TB 


product  in  the  range  of  1^00.  This  structure  ir  planar,  it  uses  no 


semiconductor  and  it  looks  very  promising. 


3.2  OTHER  SIGNAL  PROCESSING  APPLICATIONS 


The  output  of  a time  invariant  linear  filter  is  the  convolution 


of  the  input  signal  with  the  filter  impulse  response.  In  the  acoustic 


parametric  convolver,  we  can  consider  one  of  the  inputs  as  a reference 


or  '-filter  impulse  response"  for  the  other.  The  system  is  thus  a linear 


filter  with  variable  impulse  response. 


The  general  purpose  of  filtering  a signal  is  to  extract  informa- 


tion of  interest  by  discriminating  against  other  information  contained 


in  the  signal,  or  against  spurious  signals.  If  the  shape  - say  f(t)  - 


or  the  spectrum  of  the  signal  is  known,  but  not  its  time  of  occurrence. 


there  exists  an  optimum  filter  - the  "matched  filter"  - which  maximizes 


the  output  signal-to-noise:  its  impulse  response  h(t)  is  the  time 


inverted  version  of  the  signal  it  is  matched  to,  i.e.,  h(t)  = f(-t) 


This  filter  is  also  sometimes  called  a correlation  filter,  since  its 


output  is  the  autocorrelation  of  the  input  signal,  i.e.. 


g(t)  = h(r)  . f(r+t)dT  . 
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(5.2-1) 


The  parametric  convolver  is  well  suited  for  this  signal  detection 
or  signal  recognition,  provided  the  signal  is  inverted  in  time  before 
it  is  used  as  the  reference.  We  shall  give  examples  of  such  opera- 
tions for  which  the  signals  f(t)  are  coded  so  that  the  autocorrela- 
tion output  is  sharply  peaked  with  low  sidelobes. 

Symmetrical  V-FM  chirp  signals  are  especially  easy  to  use  because 
of  their  even  parity.  In  the  example  shown  in  Fig.  5. 2-1,  the  input 
signal  is  gated  on  for  5*7  nsec  and  swept  in  frequency  from  10^4  MHz 
to  111  MHz  and  back  to  104  MHz.  The  output  is  a compressed  pulse 
whose  3 dB-wi..th  should  theoretic  ,j.iv  be  0.44/Af  , where  Af  is  the 
frequency  excursion. The  measured  3 dB-width  is  60  nsec,  a result 
in  good  agreement  with  the  theoretical  value  of  63  nsec.  This  experi- 
ment corresponds  to  a 90^1  compression  ratio.  The  maximum  compression 
ratio  is  approximately  2.3  times  the  Time-Bandwidth  product  of  the 
system. 

Digital  codes  have  also  been  used;  among  them,  the  biphase  Barker 
sequences.  They  have  the  property  that  their  autocorrelation  product 
has  a peak  with  amplitude  equal  to  the  number  of  bits  and  relative 
sidelobes  of  -1,  0,  and  +1  only;  all  the  lobes  are  triangular  in 
shape.  Examples  of  7”bit  and  11-bit  Barker  code  autocorrelations  are 
shown  in  Fig.  3-2-2.  The  nonuniforraity  of  the  interaction  region 
(±  1.5  dB ) is  found  to  introduce  supplementary  phase  inversions, 
which  show  up  as  extra  zeroes  in  the  output  and  sidelobes  too  low  by 
2 dB  in  the  11-bit  code  case. 
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To  operate  the  acoustic  convolver  as  a correlation  filter^  the 
reference  signal  has  to  be  time  inverted  [see  Eq,  (3,2-1)],  The  time 
inversion  of  digital  codes  is  generally  done  by  simply  generating  the 
inverted  code,  but  the  inversion  of  analog  noniviranetrical  codes  is  more 
complex.  The  parametric  convolver  can  also  perform  this  task  and  true 
correlators  can  be  built  by  cascading  two  devices. 

While  the  input  signal  is  entirely  under  the  semiconductor,  a 
5-function  or  short  pip  is  applied  on  the  semiconductor  at  twice  the 
input  frequency.  Through  the  space  charge  nonlinearity,  an  idler  is 
generated  carrying  the  modulation  of  the  input  signal.  Because  of  the 
phase  relations  between  the  pip  (2cd,  k=0)  and  the  signal  (cu,k),  the 
idler  propagates  backward  (m,  -k)  and  it  will  thus  be  an  inverted 
replica  of  the  input.  Since  the  duration  of  the  pip  T , is  small 

pip 

compared  to  the  delay  under  the  semiconductor,  T^  , the  efficiency 

of  the  device  is  reduced  from  F by  a factor  of  (T  /T  In 

pip'  A ' 

the  particular  case  shown  in  Fig,  3,2-5,  the  pip  has  a 5 dB  width  of 
250  nsec  and  a peak  power  of  22  dBm,  T^  = 5,5  ^sec  and  F^  = -45  dB  ; 
the  inverted  signal  is  found  to  be  5I  dB  below  the  input.  This  is 
within  1 dB  of  the  theoretically  predicted  50  dB  . 

Because  of  the  high  efficiency  of  the  semiconductor  convolvers, 
it  is  possible  to  cascade  two  such  devices  and  build  a true  correlator. 
This  is  illustrated  in  Fig,  5.2-3(b),  A first  device  is  operated 
simultaneously  as  a time  inverter  and  a time  delayer:  the  two  signals 
are  then  convolved  in  a second  device.  The  various  waveforms  are 
shovm  in  Fig.  3 -2-3 (a)  in  the  case  of  an  asymmetrical  triangle.  The 
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FIG.  3.2-3--Time  inversion  and  autocorrelation  of  an  asymmetrical 
triangle  with  acoustic  device. 


- 87  - 


total  insertion  loss  is  approximately  80  dB,  but,  due  to  extra  sources 
of  noise,  the  dynamic  range  is  limited  to  some  20  dB . 

Among  the  other  possible  applications  of  the  acoustic  convolvers, 
let  us  mention  the  parametric  amplifiers.  We  have  attempted  to  build 
backward  and  forward  propagating  acoustic  wave  amplifiers  of  this  type. 
I.et  us  consider  two  forward  propagating  waves  * ith  frequency 
and  propagation  constant  , k^  = , respectively, 

where  v(m)  is  the  surface  wave  velocity  which  is,  in  general,  a func- 
tion of  frequency  in  a layered  medium.  As  these  waves  propagate  under 
the  semiconductor,  they  are  parametrically  coupled  and  a third  forward 
propagating  wave  is  generated  at  the  frequency  and  with 

a propagation  constant  k^  = od^/v(ui^)  . Another  forward  wave  is  also 
generated  at  the  sum  frequency  but  we  shall  not  consider  it  here.  If 
the  three  waves  are  well  phase  matched,  i.e.,  if  = k^L  - (k^-kj^)! 
is  small,  amplification  of  the  third  wave  is  possible  with  low  input 

powers.  Unfortunately,  in  the  airgap  structure,  the  gap  introduces 

28 

a very  large  dispersion;  for  example,  on  LiNbO,  with  a gap  height 
of  1000  X,  the  phase  mismatch  is  as  high  as  lUO  Tl/cm  or  O.U5  tt  per 
wavelength  when  ~ ~ 100  MHz  and  ~ 200  MHz.  Even  with  a 

gap  of  5000  X,  the  mismatch  is  still  some  70  7T/cm  (or  0.25  7T  per  wave- 
length). The  threshold  powers  required  for  amplification  are  extremely 
high  and  practical  devices  are  not  feasible. 

To  avoid  this  difficulty,  we  have  tried  to  operate  the  amplifier 
in  the  degenerate  backward  wave  mode.  In  this  mode,  the  input  is  a 
forward  propagating  wave  at  the  frequency  and  with  propagation 
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kj  . u is  pus,pod  by  applying  a spatially  u„i£„™ 

signal  at  m = p,,. 

2 on  the  semiconductor.  Throueh  fho 

nrough  the  parametric 

coupling,  a backward  propaeatint.  • 

P opagating  wave  is  generated  at  frequency 

“P  -i  ^ and  „itb  the  propagation  constant  -k^  ; this  wave 

Stows  sinnsoidally  with  distance,  rbia  ptincipm  was  tried  oet 

expert„entally  and  a backward  growing  wave  was  detected.  However, 

system  was  unstable  because  of  the  difficulty  of  keeping  the 

-PUS  Signal  and  the  pump  in  synchronism  and  at  fined  freguencies 

powers  are  required  to  obtain  amplification. 

The  dispersion  problems  of  both  the  forward  and  backward  wave 

parametric  amplifiers  can  probably  be  eliminated  with  some  effort 

However,  as  the  convolution  and  optical  imaging  devices  appear  to  be 

Of  far  more  practical  importance  at  the  present  • • 

cne  present  time,  we  decided  to 

investigate  them  further  rather  th=.„ 

rather  than  conce.  ■ rate  on  the  problems  of 

parametric  amplification. 


^ -to  practical  application  is  the  Direct  Fourier  Transform  of 
■ • - 38 


oloctrical  signals,  demonstrated  by  Otto.38  , compression 

experiment  is  set  up  in  which  the  signal  to  be  transformed  modulates 
the  amplitude  of  one  of  the  chirps  The  f 

Ps.  The  frequency  scale  on  the  output 

IS  a linear  function  nf  me  u- 

ction  of  the  chirp  rate  and  can  thus  be  varied,  be 

shall  not  detail 

process  any  further  since  a generalised  approach 
- oroustic  transforms  will  be  described  in  the  negt  chapter. 

The  parametric  effects  in  me 

oots  in  the  tomiconductor-piesoelactrlc  struc- 

u-er.  The  expansion  in 

8.  (2.2-5)  can  be  carried  out  to  higher  order.  The  third-order 


- 89  - 


^ ■ 

4 

i I 

! 


nonlinear  term,  in  particular,  allows  for  the  coupling  of  three  waves 
We  have  demonstrated  such  a coupling  by  using  the  scii iconductor  con- 
volver as  fallows.  Two  surface  waves  are  launched  in  opposite  direc- 
tions at  the  same  freneency  ; the  corresponding  input  electrical 

“l 

powers  are  assumed  to  be  equal  to  , In  addition,  a spatially 

uniform  potential  is  generated  across  the  semiconductor  by  applying 
to  the  semiconductor  electrode,  a signal  of  frequency  cx>  and 
power  . Via  the  third-order  nonlinearity,  spatially  uniform 

potentials  are  generated  at  the  frequencies  2co^  ± m . The  corre- 
spending  detected  powers  are  P^  . Since  the  same  port  is  used 

to  feed  in  some  signals  and  to  detect  others,  a careful  filtering  is 
required  to  avoid  any  stray  mixing  outside  the  semiconductor.  In  an 
experiment  carried  out  to  detect  third-order  effects,  we  used  the 
device  whose  convolution  characteristics  are  given  in  '.he  first 
column  of  Table  5.1.  The  input  frequencies  are  mj^/27r  = 105  MHz  and 
io^/27:  = 195  MHz  . The  nonlinear  output  is  detected  at  (2m  -m  )/2n 

2mj-m  ^ ^ 

- 15  MHz  . A small  signal  analysis  predicts  P^  ^ to  be  of  the 
form 


2m  -m  CO  co  2 

"■e 


In  fact,  the  output  power  is  experimentally  gi\en  by  the  relation 


^^”“2  -11  “2  “l  0! 

P-  " = x(P,l  „) 


e nW 


where 


e nW  ' e niW 


a = I.61H  - 0.175  log  (P^-  ) . 

e iriW  ' 
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This  expression  shows  a saturation  mechanism  which  increases 
with  the  power  applied  to  the  semiconductor.  By  contrast,  acoustic 

o)i 

powers  up  to  = 25  dBm  did  not  contribute  any  further  satura- 

tion effect.  In  any  event,  the  amplitude  of  the  third-order  non- 
linearity is  fairly  high  and  it  can  be  the  basis  for  a new  type  of 
signal  processing  device.  For  example,  the  acoustic  surface  waves 
can  be  used  to  scan  a nonuniform  electric  field  distribution  at  the 
frequency  ^ present  at  the  interface  between  semiconductor  and 
piezoelectric.  If  this  field  is  set  up  by  a bulk  acoustic  wave  in 
the  piezoelectric  substrate,  an  acoustic  imaging  device  could  be 
built.  Let  us  note,  however,  that  in  the  case  of  an  airgap  device, 
a bulk  wave  impinging  normally  on  the  free  surface  of  the  piezo- 
electric material  would  not  give  rise  to  a strong  normal  component 
of  electric  field  in  the  airgap. 


- 91  - 


CHAPTER  IV 


OPTICAL  SIGNAL  PROCESSING  AND  IMAGING 


U.l  INTRODUCTION 

The  coupling  between  acoustic  wave  and  the  semiconductor  in  the 
piezoelectric-semiconductor  structure  depends  on  the  electronic  proper- 
ties of  the  semiconductor  surface;  specifically  the  conductivity  and 
the  surface  state  charge  density.  The  former  parameter  affects  the 
nonlinear  convolution  coupling  and  the  propagation  loss  directly,  the 
latter  parameter  controls  the  dc  surface  potential  and  hence,  in- 
directly, the  interaction  strength.  Because  of  the  photoconductivity 
of  the  semiconductor,  these  surface  parameters  vary  locally  according 
to  the  intensity  of  any  incident  illumination  on  the  semiconductor, 
ihe  semiconductor  can  therefore  play  the  role  of  an  optical  image 
sensor.  In  order  to  turn  the  device  into  an  imaging  system  or,  in 
other  words,  a solid  state  camera,  one  needs  a scanning  mechanism. 

One  of  the  simplest  scanning  mechanisms  is  schematically 
represented  in  Fig.  l|.l-l(a).  We  denote  by  l(z)  the  optical 
i.ntensity  pattern  illuminating  the  semiconductor.  If  f(t)  exp(jcjDt) 
and  g(t)  exp(jcrt)  are  the  two  input  electrical  signals,  the  con- 
volver output  is  given  by  the  integral 

h(t)  = p-xp(2j0Dt)^f  (t  - z/v)  • g(t  + z/v)  • B(z)dz  (U.1-1) 
where  the  nonlinear  coupling  coefficient  is  now  a function  of  space 
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FIG.  1|.  1-1-Two  techniques  used  for  scanning  optical  images. 

(a)  The  convolver  technique  with  output  at  2cd- 

(b)  The  Quate  technique  (reference  8).  ’ 
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E(^)  because  of  the  illumination.  To  first  o 
as  follows : 


rder,  we  can  write  it 


B(z) 


Bq  + b • I(z)  , 


(^.1-2) 


If  now  lot  g(t,  ge  0 constant  a.pUtndo  signal  or  a long  pulse, 
nnd  f(t)  bo  a vary  narrow  pulse  which  can  approKicato  the  Dirac 
function  6(t)  , the  narrow  pulse  will,  i„  oifect,  scan  the  semi- 

conductor  at  the  acoustic  wave  velocity  v and  t„.  output  will  b, 
given  by  the  relation 


h(t)  » exp(aj-ut).v.B(vt)  a exp(2Jcct).v.[B^.b 


(^.1-5) 


that  the  image  will  appear  as  the  modulation  of  the  output 
waveforc.  Several  other  Scanning  schea,es  ewist  and  they  have  been 
reviewed  by  Kino.^**  ihey  either  wahe  use  of  the  light  dependence  of 
the  convolution  effect  as  in  the  scheme  Just  presented,  or  of  the 
propagation  loss  Introduced  by  the  semiconductor.  The  first  device 
built  on  this  principle  was  due  to  Quate  and  his  co-workers.®  Its 
operation  is  represented  in  Fig.  h.l-l(b,.  They  sent  a long  pulse, 
called  the  reading  pulse,  along  the  uelay  line  under  the  semiconductor. 
This  pulse  was  attenuated  because  of  the  presence  of  the  smaiconductor . 

second  narrow  pulse,  called  the  scanning  pulse,  was  sent  along  the 
delay  line  in  the  opposite  direction.  When  this  pulse  was  present  at 
s particular  point,  its  interaction  with  the  semiconductor  changed  the 
surface  condition  in  this  region,  and  hence  the  attenuation  of  the 
reading  pulse.  I„  addition,  in  the  presence  of  light,  the  attenuation 
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depended  not  only  on  tl.e  ampUtnde  of  the  scanning  pulse  but  also  on 
the  carrier  density  at  the  point  0,ere  the  scanning  pulse  was  present, 
Thus  the  amplitude  of  the  reading  pulse  as  a function  of  tl„e  depended 

directly  on  the  spatial  distribution  of  carrier  density  detem.lned  by 
the  illumination  of  the  image. 

In  all  these  syste„.s,  it  will  be  noticed  that  there  is  an  output 
even  in  the  absence  of  illumination.  In  the  presence  of  illumination, 
this  output  is  modulated  and  the  modulation  should  be  of  limited  depth 
to  conserve  the  linearity  of  the  photodetection  system.  Thus  the 

IS  not  too  different  from  the  "current"  with  illumina- 
tion, and  such  systems  have  a small  dynamic  range,  m addition,  m 
ell  the  imaging  systems  where  a short  pulse  is  used  to  scan  the  image, 
the  resolution  is  limited  by  the  width  of  the  scanning  pulse.  The 
narrower  this  pulse,  the  better  the  resolution  but  the  lower  the  input 
energy  and  the  worse  the  output  signal  power.  There  will  therefore  be 
a trade-off  between  resolution  and  output  signal -to-„oise  ratio. 

Finally,  the  scan  rate  of  such  systems  is  fixed  and  it  is  equal  to  the 
acoustic  wave  velocity. 

In  order  to  solve  some  of  these  difficulties  and  thus  build  a 
useful  device.  Quote  and  his  co-workers  made  use  of  the  storage  in 
surface  states  and  the  variation  of  the  decay  time  of  these  states 
with  lllmnlnation  to  enhance  the  sensitivity  of  the  imaging  device. 

As  we  discussed  in  Chapter  II,  the  properties  of  surface  states 
depend  very  much  on  the  surface  preparation  and  they  are  difficult 
to  control,  especially  when  the  surface  of  the  semiconductor  is  left 
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unprotected  from  contamination.  For  these  reasons,  we  decided  not 
to  use  them  and,  instead,  to  take  different  approaches  to  eliminate 
difficulties  mentioned  previously. 

The  methods  adopted  are  presetted  in  the  rest  of  this  chapter, 
with  the  purpose  of  showing  the  various  principles  involved.  In 
Section  h.2,  we  describe  a simple  means  of  lowering  the  "dark  cur- 
rent" to  arbitrarily  low  levels.  A study  of  the  dark  current  level 
is  given  and  comparison  between  experiment  and  theory  made.  The 
dynamic  range  is  then  only  a function  of  the  photoconductivity  of  the 
semiconductor  and  of  its  effect  on  the  nonlinear  coupling.  We  have 
obtained  dynamic  ranges  of  30  to  35  dB  with  a threshold  illumination 
of  50  to  100  ^lW/cm  . Low  illumination  threshold  operation  (1  pW/cm^) 
was  also  obtained  by  using  storage  in  surface  states. 

A detailed  presentation  of  some  of  the  processes,  in  particular, 
the  carrier  generation  and  transport  in  the  illuminated  semiconductor, 
is  rather  involved  and  would  require  a very  extensive  experimental 
study.  The  results  presented  here  are  not  complete  but  with  the  help 
of  Chapter  II,  they  allow  for  a quantitative  understanding  of  most  of 
the  phenomena  involved  in  the  modulation  of  the  convolution  efficiency 
by  light  and  in  the  associated  signal  processing. 

In  Section  U.3  we  show  how,  in  the  semiconductor  convolver,  we 
are  able  to  obtain  directly  Fresnel  and  Fourier  spatial  transforms  of 
the  optical  image.  The  original  images  are  then  reconstructed  by 
carrying  out  the  Inverse  transform.  In  the  case  of  the  Fresnel  trans- 
form, this  is  done  by  inserting  the  output  signal,  which  is  in  the  form 
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In  the  case  of  the 


of  an  rf  chirp^  into  a dispersive  delay  line. 

Fourier  transform,  a spectrum  analyzer  or  a second  convolver  is 
used  to  obtain  the  inverse  transform.  The  employment  of  such 
transforms  gives  a major  improvement  in  signal-to-noise  ratio 
because  of  the  large  scanning  times  used  for  one  spot.  Further- 
more, these  techniques  allow  for  an  arbitrarily  narrow  output 
bandwidth  and  an  arbitrary  scan  rate.  .Such  features  are  desirable 
for  transmitting  signals  over  links  of  limited  bandwidth,  for  elimi- 
nating inessential  information  from  a picture,  and  for  differentia- 
tion against  impulse  noise.  Several  systems  have  been  tested  to 
verify  these  principles.  With  the  device  and  signal  characteristic 
used,  the  theoretical  resolution  limit  of  0.25  mm  has  been  achieved 
and  two-dimensional  pictures  are  presented. 

h.2  IMPROVED  SENSITIVITY  IMAGING  DEVICE 

The  basic  device  used  for  optical  imaging  has  been  adapted  from 
the  one  represented  in  Fig.  5. 1-1.  The  piezoelectric  is  well  polished 
on  both  sides  in  order  to  allow  the  light  to  reach  the  semiconductor 
active  surface.  In  addition,  the  semiconductor  surface  is  thermally 
oxidized  to  reduce  the  surface  recombination  centers  and  thereby 
increase  the  light  sensitivity.  Although  BGO  might  have  seemed  a 
good  choice  for  its  good  mechanical  characteristics,  it  had  to  be 
rejected  because  of  its  very  strong  optical  activity. ^5  when  a BGO 
crystal  is  illuminated,  low  mobility  charges  are  induced  within  it 
and  they  oscillate  at  a low  frequency  which  depends  on  the  light 
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intensity  and  the  presence  of  a dc  electric  field.  These  charges  set 
up  internal  fields  which  are  detected  by  the  semiconductor  and  it  is 
as  if  the  interaction  uniformity  of  the  convolver  were  highly  unstable. 

^•2.1  Dark  Signal  Reduction  — General  Principle 

Since  the  silicon  used  is  photoconductive,  illuminating  it  modu- 
lates the  carrier  density  and  hence  the  nonlinear  coupling.  The 
simplest  way  to  read-out  an  image  is  to  insert  a narrow  rf  pulse,  the 
scanning  pulse,  at  one  end  and  a long  pulse  at  the  other  end  of  the 
delay  line.  The  convolution  output  vtIII  be  modulated  and  the  output 
as  a function  of  time  corresponds  to  the  illumination  along  one  line 
of  the  image.  One  of  the  original  difficulties  of  this  system  was  its 
limited  dynamic  range  because  of  the  inherently  high  output  with  no 
light  present  (the  dark  current).  To  make  the  dark  current  very  small, 
we  passed  the  light  incident  on  the  semiconductor  through  a coarse 
optical  grating  of  period  I , as  shovm  in  Fig.  4.2-1.  This,  in 
effect,  introduces  a wavenumber  K into  the  light  beam,  such  that 
K£  = 2rr  . 

Now  the  optimum  output  is  obtained  with  two  input  signals  of 
frequencies  ^ propagation  constants  of  k^  , respec- 

tively, such  that 


k^  - k^  = ± K 

“l  ' “2  " ’ 


(4.2-la) 


(4.2-lb) 
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where  v is  the  velocity  of  the  acoustic  surface  wave.  By  using 
a nondegenerate  system  with  input  frequencies  displaced  by 

the  amount  given  in  Eq . (4.2-lb),  a strong  output  is  obtained  when 
the  device  is  illuminated,  with  a very  smiall  output  when  -.nere  is 
no  light.  Thus  the  "dark  current"  is  caused  to  be  very  small.  It 
should  be  noted  that  if  three  colored  gratings  of  different  peri- 
odicities are  used,  the  device  can  be  arranged  to  give  outputs  at 

three  different  frequencies  corresponding  to  the  three  primary 
colors . 

In  order  to  evaluate  the  "dark  current"  reduction,  we  consider 
the  input  signals  as  being  two  square  pulses  of  frequency  oj^  and 
, respectively,  and  of  duration  smaller  than  the  delay  under 

the  semiconductor.  Recalling  the  results  of  Section  2.1.1,  and  con- 
sidering again  the  center  of  the  semiconductor  as  the  origin  z = 0 , 
the  output  convolution  voltage  is  of  the  form 

V = exp  [ +a>2)t]  y^exp  [j(k^-k2)z]  dz 

wh?re  the  integral  is  evaluated  over  the  region  of  overlap  of  the  two 

acoustic  pulses.  We  find  that  V is  given  by  the  following  expres- 
sions ; 

+U)2)t]  • sin  27TZXFt  for  0 < t < T /2 

A 

(4.2-2) 

V = ^ sin  27tAF(T  -t)  for  T /2  < t < T, 

" A'  — — A 

where  AF  = (co^  - c^^)/2n  . 
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Figure  h.2-2  shows  the  envelope  of  the  output  waveform  for  various 


values  of  the  product  AF  • T . When  AF  = 0 , the  output  has  a tri 

angular  shape;  it  extends  for  a time  T and  its  maximum  amplitude  is 
V • . This  is  the  usual  result  for  degenerate  convolution.  When  a 

frequency  offset  is  introduced,  AF  / 0 , the  output  waveic  rm  is  com- 

posed of  two  symmetrical  portions  of  sine  waves.  For  very  small  off- 
sets in  frequency  satisfying  the  condition  • T < 0.5  , the  ampli- 

tude  is  maximum  at  the  center  of  the  waveform  and  the  maximum  ampli- 
tude is  vT^  sin  (ttZ^  • ) / (ttAF /T^ ) . As  AF  increases,  the  maximum 

dark  current  becomes  vT  /(ttAF  *T  ) 

The  reduction  R in  dark  current  brought  by  the  introduction 
of  the  offset  f.n  frequency  AF  is  therefore 

R = 20  log  (ttAF.T^)  (dB)  (U.2-5) 

as  soon  as  is  greater  than  0.5.  If  ZF  • T^  < 0.5  , it  will 

be  seen  that  R is  given  by  the  relation 

TTAF  . T 

R = 20  log  (dB) 

sin(7TAF  • T^) 

In  this  case,  R is  4 dB  or  less,  which  is  insufficient  in  practice. 
As  a consequence,  the  only  useful  formula  is  given  by  Eq.  (4.2-5). 

The  dark  signal  was  measured  with  two  devices  operated  in  this 
nondegenerate  mode.  The  input  signals  are  square  pulses  of  duration 
5 or  4 |is;  one  of  the  input  frequencies  is  fixed  at  105  MHz;  the  other 
allowed  to  vary  from  105  to  II5  MHz.  The  results  are  plotted  in 
Fig.  4.2-5  vs  the  frequency  offset  ZF  . Also  represented  are  the 


- lOI  - 


iiff  I 


' —111  '"iTkii  ii  - 


AF.Ta  = 0 


FIG.  4 .2-2--0utput  convolution  of  two  square  pulses  of  duration 
and  whose  center  frequencies  are  offset  by  an  amount 
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FIG.  4.2-5“-Dark  signal  reduction  du>.  to  a frequency  offset.  Two  sets 
of  experimental  da.:a  are  shown  and  compared  to  the  theo- 
retical curves  (full  lines). 


theoretical  variations  for  the  two  cases  . J os  and  = h ms 
It  will  be  seen  that  the  experimental  data  are  within  1 dB  of  the 
theoretical  predictions.  In  addition,  with  1,  ps  square  pulses  and 
« - 10  MHz  , the  dark  signal  reduction  is  1.J  dB.  The  dark  signal 
only  decreases  as  1/aF  because  of  the  relatively  large  high  fre- 
quency components  of  the  square  pulses  used.  If  a lower  dark  signal 
Is  required,  one  should  reduce  the  spectrum  of  the  inputs  so  that, 
when  the  offset  is  Introduced,  these  spectra  do  not  overlap.  This 
is  easily  achieved  by  rounding  off  the  edges  of  the  pulses  used. 
Ideally,  If  the  spectra  of  the  input  signals  are  limited  to  20/t^  , 
these  signals  will  still  very  much  resemble  square  pulses  but  the 
dark  current  will  be  zero  as  soon  as  m is  greater  than  koA^ 

If  = 4 US  ^ this  means  > lo  MHz 

The  transducers  used  in  these  experiments  have  a center  fre- 
.quency  of  lo8  MHz  and  a 3-dB  bandwidth  of  20  MHz.  They  will  there- 
fore  introduce  some  extra  loss  wi.en  they  are  operated  above  118  MHz 
or  below  98  MHz.  In  order  to  avoid  this  further  effect,  we  limited 
the  experimental  frequency  range  to  10  MHz  and  we  fixed  one  of  the 
input  frequencies  at  105  MHz.  To  accommodate  larger  frequency 

offset,  transducers  with  unequal  center  frequencies  are  used  as  we 
shall  see  in  Section 

IlRht  Sensitivity  .pj  Dyn.amic  Banco 

In  this  sub-section,  we  gi  e a more  quantitative  description  of 
the  affect  of  illumination  on  the  behavior  of  the  convolver.  The 
principal  process  is  the  carrier  generation  in  the  semiconductor. 
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Depending  on  the  surface  condition  prior  to  illumination,  different 
types  of  phenomena  can  take  place,  A complete  presentation  is  rather 
involved  because  of  the  difficulty  of  evaluating  some  of  the  optical 
parameters  like  the  quantum  efficiency  and  spectral  response,  and 
because  the  thermal  equilibrium  assumption  does  not  hold  when  the 
semiconductor  is  illuminated. 

As  a first  step,  we  study  the  cuse  of  a semiconductor  at  flat 
band  or  in  the  accumulation  regime  before  the  illumination  is  present. 
This  mode  is  relatively  simple  and  corresponds  to  the  case  of  oxi- 
dized n-type  silicon  in  the  absence  of  bias.  In  this  case,  the  light 
causes  the  carrier  concentration  to  increase  by  an  amount  AN  , The 
dependence  of  the  convolution  on  AN  is  studied  theoretically.  By 
making  simplifying  assumptions,  we  derive  the  general  form  of  AN 
Its  magnitude,  its  dependence  on  the  frequency  offset,  and  the  dif- 
fusion length  in  the  semiconductor  are  presented  and  compared  to 
experiments . 

As  a second  step,  we  study  the  response  to  light  of  the  imaging 
device  as  a function  of  the  type  of  semiconductor  used  and  the  dc  bias 
that  may  be  applied  to  it.  Experimental  convolution  vs  light  charac- 
teristics are  discussed  qualitatively  on  the  basis  of  simple  physical 
arguments . 


Li^ht sensitivity  in  flat  band  or  accumulation  case 

Recalling  the  results  of  Chapter  II,  we  find  that  the  convolution 
output  is  made  up  of  two  terms:  1)  the  nonlinear  coupling  term  (the 

U-value)  which  has  been  seen  to  be  essentially  inversely  proportional 
to  the  free  carrier  concentration,  when  the  semiconductor  surface  is  at 
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flat  band  or  accumulated,  and  2)  the  propagation  loss  term  which 
varies  exponentially  with  the  surface  resistivity. 

Let  us  now  assume  the  semiconductor  surface  is  either  at  flat 
band  or  in  accumulation,  so  that  we  can  define  an  equivalent  surface 
CO  centration  N . When  the  convolver  is  illuminated,  electron-hole 
pairs  are  created  and  the  electron  carrier  density  is  perturbed  by  an 
amount  M , and  the  hole  carrier  density  by  an  amount  AP  . For  a 
semiconductor  at  flat  band,  the  carrier  concentration  remains  uniform 
and  we  can  write  = AP  . If  the  semiconductor  surface  is  accumu- 
lated, internal  fields  exist  in  ^he  material,  the  carrier  concentra- 
tion is  not  uniform  and  the  relation  AN  = ^ is  not  strictly  true. 

In  what  follows,  we  shall  be  dealing  mainly  with  weak  accumulation 
corresponding  to  a surface  potential  of  the  order  of  ^/kx  = 1 
in  the  absence  of  light.  The  bulk  of  the  semiconductor  will  stay  at 
flat  band  under  illumination  and  M = AP  there.  By  writing  the 
expression  for  the  total  charge  to  be  neutral,  before  and  after  illumi- 
nation, we  can  evaluate  the  excess  majority  carrier  concentration  at 
the  surface  for  a given  excess  in  the  bulk.  The  details  of  this  calcu- 
lation are  not  presented  here.  It  can  be  shown  that  the  excess  carriei 
concentration  profile  in  the  semiconductor  remains  uniform  to  within 
15  to  20<^.  The  departure  from  uniformity  is  maximum  for  an  excess 
carrier  density  equal  to  5 times  the  bulk  concentration  in  the  dark. 

For  lower  and  higher  illuminations,  the  uniformity  condition  is  almost 

verified.  And  in  the  following  we  shall  therefore  consider  ^ as  being 
uniform  and  equal  to  AP 
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In  the  experiments,  the  grating  is  composed  of  dark  stripes  m cimes 
wider  than  the  white  stripes;  the  spatial  duty  cycle  is  therefore 
a = l/(nvH)  . Because  of  the  grating,  the  photon  flux  $(z)  varies 
from  zero  to  . it  can  be  expanded  in  the  following  Fourier  series 


where 


“^0  = 


and  $ = 2a  ^ ^ 

**  7Tna  M ^ n j/f  0 


The  perturbation  « can  also  be  represented  by  a aerlea  aa  follona: 


= \ ^ coa(nKz) 


(^.2-5) 


The  coefficlenta  are  linearly  related  to  through  a factor 

F(nK)  which  accounts  for  the  diffusion  effects  In  the  semiconductor. 

For  small  nK  , F la  essentially  constant  and  as  nK  increases,  the 
diffusion  becomes  stronger  and  the  corresponding  perturbation  In  car- 
rier density  Is  reduced.  An  expression  for  F(nK)  will  be  given  later. 

If  two  acoustic  waves  are  propagating  In  opposite  directions,  when 
they  reach  any  point  under  the  semiconductor,  the  product  of  their  ampli- 
tudes will  have  been  attenuated  by  an  amount  given  by  the  following 
expression: 


Q!j.(z)  dz  I A exp(-a^L) 
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where  L is  the  total  length  of  the  seraicondu''tor . The  propagation 
loss  is  therefore  not  a function  of  position  z but  depends  on  the 
average  illumination  of  the  semiconductor.  By  contrast,  the  nonlinear 
coupling  is  going  to  depend  on  z and  if  , cjd^  and  are 

the  center  frequencies  and  the  propagation  constants  of  the  two  acoustic 
waves,  respectively,  the  output  convolution  signal  appears  at  the  fre- 
quency . Its  amplitude  is  given  by  Eq.  (2.1-4)  to  be  V 

where 


L/2 


V oc  exp  (-a 


p.)  / ; 

I /o  N + 


exp[j(k2  ■ 


dz 


-L/2  ^ ^ ^'^n 


(4.2-6) 


If  the  input  signal  center  frequencies  are  identical,  - k^ 
and  the  convolver  is  operated  in  a degenerate  mode;  it  responds  to  the 
dc  component  of  the  inverse  of  the  concentration;  hence,  for  low  level 
injection  (AN^  « N)  , the  output  voltage  is  given  by  the  following 
expression 


1 / 

V,  oc  exp(-a  L)  • — I 1 - I . ^4.2-7) 

deg  r ^ ^ ^ ^ j ~ n 

In  the  high  level  injection  mode,  the  high  frequency  component 

n >0 

will  contribute  to  the  dc  component  of  l/(N  + AN)  . A general  evalua- 
tion requires  a complete  Fourier  analysis.  If  we  neglect  the  diffusion 
effects,  the  factor  F(nK)  is  independent  of  nK  . M will  there- 
fore follow  exactly  the  photon  flux;  it  will  vary  stepwise  between 
zero  and  . The  inverse  concentration  varies  stepwise  between 

1/N  and  l/(N  + Z^)  . Its  dc  component  is  straightforward  to 
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evaluate  and  leads  to  the  following  expression  for  [V(degenerate) ] : 


''deg  “ - (l  - 


N + /^^ 


(4.2-8) 


If  the  input  signal  center  frequencies  are  offset  by  Zku  = ± Kv 
the  convolver  only  responds  to  one  of  the  spatial  harmonics  of  1/(N  + ^)  ; 
and,  in  the  low  level  injection  case,  the  output  voltage  is  proportional 


to  as  follows 


non  deg 


1 Z\N 

« exp(-a  L) 4 

N 2N 


(4.2-9) 


In  the  high  level  injection  case  and  in  the  absence  of  diffusion,  we 


find  that 


_ 1 aAN  sin  Tta 

''non  deg  “ ^ 

N N + 7Ta 


(4.2-10) 


To  test  the  validity  of  the  previous  relations,  we  used  the  con- 
volver whose  characteristics  were  sumnarlted  In  the  first  column  of 
Table  2.3.  The  semiconductor  Is  an  epitaxial  silicon  with  an 

epi-layer  concentration  of  7 X 10^°  m'5  and  an  epitaxial  layer  thick- 
ness Y = 3 „ . The  epitaxial  nature  of  the  semiconductor  was  not 

Important  In  the  study  of  the  convolution  efficiency  presented  in 
Chapter  II  because  the  nonlinear  process  takes  place  In  the  first  Debye 
length  of  the  sendconductor.  On  the  contrary,  here,  the  thickness  of 
the  epl-layer  plays  an  important  role  when  light  Is  shined  at  the  send- 
conductor. The  n-^  substrate  sets  up  a potential  barrier  which  keeps 
the  holes  from  flowing  Into  the  bulk.  The  hole  current  is  therefore 
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zero  at  y = Y . In  other  words,  all  the  minority  carriers  generated 
by  the  light  are  confined  to  the  epi-layer.  They  can  diffuse  or  recom- 
bine with  electrons  in  the  layer  or  at  its  surface.  The  diffusion 
length  of  holes  in  silicon  is  always  much  larger  than  the  epi-layer 
thickness. In  addition,  we  saw  in  Chapter  II  that  because  of  the 
presence  of  surface  states,  the  surface  of  the  semiconductor  is  naturally 
accumulated.  This  accumulation  corresponds  to  a dc  surface  potential 
^ ‘^he  absence  of  light  and  the  surface  majority  carrier 
concentration  is  approximately  I9  x 10^^  m'^.  This  enhanced  surface 
concentration  reduces  the  surface  recombination  velocity^^  and  we  may 
neglect  this  effect.  It  follows  that  the  excess  hole  concentration 
will  be  almost  uniform  across  the  epi-layer.  When  steady-state  is  reached 
we  can  equate  the  current  generation  due  to  the  illumination,  to  the  sum 

of  the  diffusion  current  along  the  surface  and  the  recombination  current 
in  the  epi-layer,  thus  writing 


^n  ^ r = AP  y/t  + (nK)^  DY  ^ 

n 


(i+.2-ll) 


T Is  the  quantum  efficiency  and  r the  absorption  ratio  in  the  semi- 
conductor,  r is  the  minority  carrier  lifetime  and  D the  diffusion 
constant;  ia  the  eKcess  carrier  concentration  corresponding 

to  the  harmonic  component  of  the  flus  . Equation  (I, .2-11)  can  be 

put  in  the  final  form 


AN  = AP  = 

n n 


*n  1 r 


^nK)==DytV/r  ° 


which  provides  an  expression  for  the  factor  F(nK) 


no  - 


1 

i 


The  light  source  employed  is  a tungsten  filament  lamp  producing 
<*lte  light.  The  Illumination  power,  , was  measured  with  a 
alllcon  photodetector  and  we  estimated  the  maximum  photon  flux  to  be 


- “l  X 10^®  photons/m^/s  , (h.2-13) 
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N = 19  X 10  n , hence  the  value  of  = 36  X 10^^  ra'^ 

And  from  Eq . (U.2-l^i)  the  lifetime  necessary  to  fit  the  experimental 
data  is  found  to  be  x = 77  ps  . The  time  x was  also  measured 
directly;  the  light  is  chopped  at  a fast  rate  and  the  rise  time  of 
the  output  waveform  is  measured.  This  method  yields  x = 70  ± lO  ps 
which  agrees  well  with  the  previous  conclusion. 

Knowing  the  value  of  , we  can  now  predict  the  level  of 

the  nondegenerate  convolution  in  the  presence  of  light,  relative  to 
the  degenerate  convolution  in  dark.  This  ratio  comes  from  Eq . (I1.2-IO) 
and  we  obtain  the  result 


non  deg  w.  light 


deg  in  dark 


= 20  log  I 


N +M, 


predicted 


M 


sin  7Ta 


TTa 


= -13.5  dB  . 
(i+.2-15) 


The  directly  measured  value  of  this  same  ratio  was  -12.5  ± O.5  dB . The 
fair  agreement  with  the  predicted  value  is  a check  of  the  consistency  of 
the  theoretical  expressions  given  earlier. 

Finally,  the  ratio  defined  in  Eq.  (I1.2-I5)  was  measured  for  various 
grating  periodicities,  i.e.,  various  frequency  offsets.  In  the  experi- 
mental set-up,  the  grating  is  illuminated  with  a white  light  source  of 
given  power  and  the  image  of  the  grating  is  focused  on  the  semiconductor 
by  means  of  a lens.  By  changing  the  relative  positions  of  the  grating, 
the  lens  and  the  semiconductor,  we  vary  the  periodicity  of  the  imaged 
grating.  But,  at  the  same  time,  the  intensity  of  the  illumination 
received  by  the  semiconductor  changes.  Later,  in  conjunction  with 
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Fig.  U.2-5,  we  shall  discuss  the  variation  of  the  nondegenerate 


convolution  power  with  the  illumination  power;  we  shall  see  that 


with  high  illumination  level,  as  is  the  case  here,  the  convolution 


varies  as  * With  this  added  data,  we  can  correct  the 


measured  values  of  the  nondegenerate  convolution  power  and  hence 


normalize  them  to  a given  illumination  level.  This  is  what  has 


been  done  and  the  results  plotted  in  Fig.  1+.2-1+  are  normalized  to  the 


illumination  level  corresponding  to  the  frequency  offset  = I.75  MHz 


As  we  have  already  discussed,  for  low  frequency  offsets,  dif- 


fusion may  be  neglected.  The  nondegenerate  convolution  power  output 


is  independent  of  the  offset;  in  addition,  the  excess  carrier  concen- 


tration is  high  and  the  "high  level  injection  theory"  should  be  used. 


For  very  large  offsets,  the  excess  concentrations  are  small  because 


of  the  strong  effect  of  diffusion;  the  low  level  injection  approxi- 


mation becomes  valid  and  Eq . (h.2-^)  can  be  used  with  as  given 

by  Eq.  (U.2-12).  With  this  approximation,  we  obtain  the  theoretical 


variation  denoted  "low  level  theory"  and  shown  in  Fig.  h.2-k.  The 


values  of  the  diffusion  constant  and  life  time  of  the  minority 


carriers  (hole)  we  used  are  D = 6.5  x lo’^  m^/s  and  t = 77  ns 


respectively.  The  experimental  points  are  seen  to  be  in  good  agree- 


ment with  the  previous  conclusions,  in  the  ranges  where  each  of  the 


two  approximations  holds.  In  the  intermediary  range,  the  diffusion 


effects  are  important  but  the  low  level  injection  approximation  is 


not  valid;  several  Fourier  components  in  the  expansion  of  z:i^(z) 


should  be  kept.  They  are  affected  differently  by  diffusion  according 
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FIG.  U.2-l|--Ratio  of  nondegenerate  convolution  power  to  degenerate 
I convolution  power  in  dark  vs  frequency  offset  for  an 

n/n  silicon  sample  under  light  accumulation  conditions 
a (Plight  - 100  W/m5). 
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to  Eq.  (k.2-l2)  and,  in  a rigorous  study,  one  should  use  Eq.  (4.2-6) 
to  evaluate  the  nondegenerate  convolution  output.  Had  we  conducted 
such  a study,  we  would  have  been  able  to  smoothly  link  the  two  theo- 
retical curves  of  Fig.  4.2-4. 

Dynamic  range  and  light  sensitivity 

We  define  the  dynamic  range  of  the  imaging  device  as  the  varia- 
tion of  the  convolution  output  power  when  the  illumination  varies 
from  the  minimum  detectable  level  to  the  saturation  level.  The  mini- 
mum detectable  illumination  level  is  also  called  the  threshold  level. 

It  is  a function  of  the  dark  signal  reduction  since  the  lower  the 
dark  signal,  the  more  sensitive  the  detection  system  and  hence  the 
lower  the  threshold  level.  It  is  also  dependent  on  the  sensitivity 
of  the  convolution  phenomenon  to  any  increase  in  carrier  concentra- 
tion in  the  semiconductor.  These  two  features,  the  dark  signal  reduc- 
tion and  the  sensitivity  to  carrier  concentration  Increase,  are  func- 
tions of  the  periodicity  of  the  grating  as  we  saw  earlier.  But,  in 
addition,  the  latter  feature  depends  in  a rather  complicated  fashion, 
on  the  nature  of  the  semiconductor  and  its  surface  potential . 

To  examine  this  dependence,  we  have  measured  the  nondegenerate 
convolution  power  as  a function  of  the  illumination.  The  plots  corre- 
sponding to  the  n-  and  p-type  silicon  samples  studied  earlier  (see 
Table  2.5)  are  shown  in  Fig.  4.2-5-  The  grating  used  has  a duty  cycle 
a = 0.25  and  the  frequency  offset  is  10  MHz.  In  addition,  the  illuii- 
nation  was  measured  at  the  semiconductor,  i.e.,  through  the  grating,  and 
the  illumination  power  used  in  Fig.  4.2-5  therefore  corresponds  to 
^ our  definition  of  P 

iight 
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FIG. 


J4^.2-5--Experimental  convolution  response  to  white  illumination. 
The  parameter  q0  /kT  is  the  surface  potential.  The 
frequency  offset  is  AF  = 10  MHz  . The  grating  duty 
cycle  is  a = 0.25  and  the  illumination  is  measured 
through  the  grating,  i.e.,  P v a . 
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In  the  absence  of  any  biasing  voltage,  the  n-type  semiconductor 
has  a slightly  accumulated  surface.  We  have  studied  earlier  in  this 
section,  how  the  light  affects  the  behavior  of  the  device  for  a given 
illumination  level.  Figure  4.2-5  (top  curve  with  q^^^/kT  = 1 ) now 
shows  the  variation  of  the  nondegenerate  convolution  with  the  light 
level:  the  threshold  level  is  0.8-1.0W/m^  and  the  dynamic  range 

20-25  dB . In  addition,  for  lower  illuminations,  the  convolution  shows 
a square  law  dependence  on  the  illumination  power,  which  is  in  accordance 
with  Eq.  (4.2-9)  and  £q . (4.2-12);  as  the  light  increases,  high  level 
saturation  occurs  and  the  convolution  power  now  varies  as 

light 

For  an  even  higher  illumination,  the  excess  density  in  Eq . (4.2-10) 

will  go  on  increasing  and  we  might  expect  V , to  level  off  and  re- 

ton  deg 

main  constant. 

When  the  semiconductor  is  initially  driven  in  strong  depletion, 

rather  different  processes  take  place  which  allow  for  a much  lower 

threshold  illumination  level.  We  consider  the  same  silicon  sample  as 

previously,  and  we  apply  a bias  voltage  so  that  its  surface  potential 

becomes  q^^Q/kT  = -20  ; a corresponding  depletion  layer  is  generated 

at  the  surface  of  the  semiconductor  with  a width  given  by  f = 6.3  X = 

d 

0.94  urn  . Under  illumination,  the  minority  carriers  (holes)  created  in 
the  depletion  layer  are  pulled  almost  instantaneously  toward  the  sur- 
face, under  the  influence  of  the  large  electric  field  present  in  the 
depletion  layer.  Concurrently,  the  electrons  move  toward  the  bulk  of 
the  semiconductor.  Due  to  tuu  presence  of  surface  states,  the  holes 
vill  either  be  trapped  or  will  recombine  as  they  reach  the  surface. 

If  Tgg  is  the  life  time  for  recombination  at  the  surface,  the  mobile 
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hole  charge  at  the  surface  per  unit  area  is  found  to  be  AP  = G • t 

s ss  ' 


where  G is  the  generation  rate  of  holes  in  the  depletion  layer,  per 


unit  area.  Thus,  if  t is  large,  and  it  can  be  several  ms,  AP 

ss  ' s 


can  be  relatively  large  and  the  corresponding  light  sensitivity  very 


high. 


If  we  call  the  attenuation  constant  of  the  light  in  the 


semiconductor,  we  find  that  AP^  is  given  by  the  following  expression: 


X 10  a • Plight  • 'I  r • 1,,  (11.2-16) 


where  I is  the  depletion  layer  width.  For  lower  illumination  levels, 
APg  will  be  smaller  than  the  total  surface  state  density  and  a steady 


state  is  reached  between  light  generation  and  surface  trapping  in  which 


the  traps  are  incompletely  filled.  We  may  rewrite  Eq . (b.2-16)  in  the 


^s  = Vd^ 


(i+.2-lT) 


where 


^ Cl  . i b X 10  a • P,  . , ^ • Ti  r • T 
+ _ _sh light  * ss 


(4.2-18) 


g is,  as  in  Eq . (2.2-2’4),  a measure  of  the  "inversion"  charge  at  the 


surface.  The  presence  of  the  trapped  charge  will  partially  screen  the 


dc  bias  field  and  reduce  the  depletion  layer  width.  Since  the  effect 


of  such  a reduction  on  the  convolution  is  minor  (see  Chapter  II ),  we 
shall  neglect  this  phenomenon  here.  On  the  contrary,  by  following  the 


same  argument  as  in  Section  2.2,  the  surface  charge  AP  will  reduce 
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the  effect  of  the  acoustic  fields  on  the  semiconductor.  This  reduction 
is  characterijied  by  the  factor  ft  defined  in  Eq . (2.2-27).  It  was  seen 
to  be  dependent  on  . The  reduction  in  degenerate  convolution  and 

consequently  the  level  of  the  no.idegenerate  convolution  will  therefore 
also  depend  on  g"*"  ; furthermore  we  are  going  to  see  that,  due  to  the 

large  time  constants  associated  with  the  surface  states,  the  light 
sensitivity  in  this  process  is  increased. 

To  put  this  result  in  quantitative  terms,  we  refer  to  Fig.  k.2-5 
(top  half).  For  a given  convolution  output  (say,  2h  dB),  we  see  that 
the  illu:.lnatlon  level  required,  is  reduced  fron,  25  to 

0.3  W/m  when  the  semiconductor  is  driven  into  strong  depletion.  In 
this  case,  the  variation  in  degenerate  convolution  between  dark  and 
light  was  measured  to  be  -2.5  dB.  The  reduction  factor  should  there- 
fore  be  some  2.5  dB  corresponding  to  g"^  = 2 . By  using  Eq.  (U.2-16) 

with  i = 0.9U  = 2 X 10^  m-1  , q = 0.5  , = 7 X 10^°  m"^  , 

r = 0.7,  Xd  = 0.15  MJn  and  = 2 ms  , we  find 

= 1.5  and  = 1.5  X 10^^  ra‘^  . 

These  results  are  in  reasonable  agreement  with  the  predicted  values. 

There  is  some  inaccuracy  in  the  determination  of  a . and  as  a conse- 

ph 

quence,  an  inaccuracy  in  the  number  of  electron-hole  pairs  actually 
generated.  We  also  note  that  AP^  is  smaller  than  the  total  surface 
state  density  which  had  been  found  to  be  in  the  10^^  m"^  range. 

As  the  illumination  is  increased.  Fig.  I1.2-5  shows  that  satura- 
tion occurs.  This  is  mainly  due  to  the  complete  filling  up  of  the 
surface  states  and  to  the  depletion  layer  reduction  which  follows. 
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Similar  measurements  have  been  made  with  bulk  p-type  silicon 


(100  n^cm)  and  the  plots  of  the  convolution  output  vs  light  are  shown 
in  the  lower  half  of  Fig.  4. 2-5-  With  no  bias,  this  silicon  sample  is 
somewhat  depleted  (q^^^/kT  = -2 ) . It  exhibits  a low  threshold  level 
(O.l  W/m  ).  As  the  surface  is  driven  into  strong  depletion 
(curve  q^^^/kT  = -22),  the  threshold  level  becomes  higher  but  the 
light  sensitivity  is  improved.  By  studying  the  pulsed  bias  response 
of  the  convolver,  as  we  described  in  Chapter  II,  the  surface  state 
time  constant  has  been  estimated  to  be  in  the  200  - 500  ps  range  com- 
pared to  2 ms  in  the  case  of  the  n-type  silicon.  The  response  of  the 
system  to  light  is  a function  of  the  minority  carrier  charge  at  the 
surface,  itself  proportional  to  G ' [see  Eq . (4.2-16)].  Since 

Tgg  has  a reduced  value,  to  obtain  the  same  response,  higher  illumi- 
nation levels,  i.e.,  higher  G , are  required  with  the  p-type  silicon. 

The  presence  of  a strong  surface  depletion  together  with  long  life 

i 

surface  states  therefore  allows  for  large  threshold  level  reductions  j 

and  increase  in  sensitivity  of  the  device.  This  feature,  although  i 

interesting,  relies  on  the  properties  of  the  surface  states,  which  | 

are  difficult  to  control  and  to  reproduce.  In  addition,  when  an  MOS  I 

structure,  like  ours,  is  illuminated,  low  mobility  charges  are  created 
in  the  silicon  oxide,  i.e.,  in  slow  surface  states,  and,  under  the 

( ' 

influence  of  the  dc  field,  dipoles  are  slowly  formed.  As  a conse-  I 

quence,  after  illumination  for  some  time,  charges  corresponding  to 

the  grating  periodicity  can  be  stored  in  the  oxide,  and  the  dark  i 

current  reduction  capability  and  hence  the  dynamic  range  of  the  system  ! 


I 


arc  lost.  This  Is  because,  even  in  the  dark,  a periodic  (field) 
pattern  is  present  in  the  system.  Such  an  oxide  storage  mechanism 
took  place  in  some  of  our  experiments.  The  read-in  time  for  such 
storage  was  some  10  to  20  minutes  and  the  decay  time  of  the  storage 
was  of  the  order  of  1 or  2 days. 

All  the  arguments  presented  here  in  conjunction  with  Fig.  k,2-5 
are  rather  qualitative  and  simplistic.  Nonetheless,  they  provide  a 
fairly  good  understanding  of  the  effect  of  light  on  a convolver  using 
a semiconductor  whose  surface  is  at  flat  band  or  is  slightly  accumu- 
lated. They  also  lead  to  conclusions  which  are  in  reasonable  agree- 
ment with  experiments. 


A more  detailed  study  would  certainly  be  of  interest  in  order 
to  predict  more  fully  the  behavior  of  the  system  as  a function  of  the 
semiconductor  surface  condition.  Such  a study  would  require  very 
careful  measurements  of  the  various  characteristics  of  the  convolver 
vs  bias  and  light.  In  order  to  make  repeatable  experiments  and  to 
get  rid  of  the  difficulties  we  have  encountered,  the  configuration  of 


the  experimental  device  should  probably  be  altered.  The  major  mechani- 
cal problems  has  to  do  with  the  presence  of  the  airgap  and  its  uni- 
formity. It  is  evident  that  the  smaller  the  semiconductor  sample,  the 
easier  it  is  to  maintain  the  uniformity  of  the  gap.  Hence  in  a physics 


oriented  research  project,  smaller  semiconductor  samples  could  be  em- 
ployed with  lengths  of  the  order  of  5 mm.  One  might  then  expect  to 
obtain  uniform  interaction  regions.  The  second  main  problem  is  related 
to  the  surface  of  the  semiconductor  and  its  passivation.  Along  all  the 
steps  of  the  preparation  of  the  samples  - cleaning,  oxidation,  scribing 
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or  cutting,  mounting  — high  standards  of  cleanliness  must  be  obeyed  to 
avoid  contamination  of  the  surfaces.  Finally,  even  with  these  precau- 
tions, the  semiconductor  will  remain  in  contact  with  the  laboratory 
atmosphere.  This  will  result  in  a long-term  irreversible  drift  of  the 
characteristics  of  the  surfaces.  This  difficulty  should  be  reduced  if 
the  experimental  device  is  put  together  and  kept  under  vacuum’ 

It  should  also  be  noticed  that  when  a "dc"  bias  voltage  is  applied 
to  the  structure,  stray  effects  may  arise  due  to  the  motion  of  charges 
in  the  oxide,  as  we  stated  earlier.  To  avoid  this  phenomenon,  the  bias 
should  be  turned  on  for  small  periods  of  time  and  the  measurement  done 
rapidly.  Finally,  a study  of  the  optical  response  should  be  done  em- 
ploying monochromatic  light,  so  that  the  absorption  constants  and  the 
illumination  power  be  well  defined.  A study  as  a function  of  light 
wavelength  is  also  very  useful  to  determine  the  various  optical  transi- 
tions in  the  semiconductor. 

From  the  results  presented  here,  we  have  seen  that  a dynamic  range 
of  20-25  dB  with  a threshold  illumination  level  of  1 w/m^  are  typical 
of  the  unbiased  semiconductor  convolver.  Larger  dynamic  ranges,  froir. 

25  to  35  dB,  can  and  have  been  obtained  with  a finer  grating  and  an 
offset  of  20  MHz.  It  should  however  be  pointed  out  that  increasing 
the  frequency  offset  any  further  will  not  necessarily  increase  the 
dynamic  range.  The  dark  current  decreases  as  but  the 

response  to  kight  for  large  offsets  decreases  as  (AF)’^  . A finer 

is  nonetheless  required  if  high  resolution  imaging  devices 
are  needed  (see  Section  U.5).  To  increase  the  dynamic  range  beyond 
j5  dB , other  methods  would  have  to  be  used.  A possible  method  is  to 
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create,  at  the  semiconductor  surface,  finely  spaced  potential  wells  to 
prevent  the  carrier  diffusion  and  hence  increase  the  Modulation  Transfer 
Factor  and  the  dynamic  range.  These  potential  wells  can  be  created  by 
deposition  of  Schottky  barriers  on  the  semiconductor  or  by  using  the 
buried  layer  technology  now  currently  employed  in  Charge  Coupled  Devices. 
By  this  method,  the  light  sensitivity  would  also  be  increased  because  of 
the  possibility  of  integrating  the  photo-injected  charge  over  long  times. 
This  charge  cannot  diffuse  and  is  therefore  stored  in  the  potential  wells. 
This  again  is  a concept  used  in  the  CCD  devices. 

4.2.3  Example  of  "Direct"  Imaging 

One  of  the  simplest  optical  imaging  devices  that  can  be  built  on 
the  principles  presented  in  the  last  two  sections,  is  represented  in 
Fig.  4.2-6(a).  The  object,  in  the  form  of  a transparency,  is  projected 
through  the  piezoelectric  substrate  onto  the  semiconductor  in  the  con- 
volver. An  optical  grating  in  front  of  the  object  (a  photographic 
transparency)  introduces  a spatial  periodicity  in  the  optical  signal. 

Its  effect  on  the  interaction  is  cancelled  out  by  introducing  a rela- 
tive offset  AF  of  the  correct  amount  in  the  center  frequencies  of 
the  two  input  signals.  A narrow  scanning  pulse  at  frequency  F is 
convolved  with  a long  pulse  at  F + AF  ; the  output  at  2F  + ZliF  then 
corresponds  to  one  line  of  the  image. 

An  example  of  such  a "direct"  imaging  with  a narrow  scanning  pulse 
of  0.15  |is  width  and  a grating  periodicity  of  550  |im  (ZliF  = 10  MHz)  is 
shown  in  Fig.  4.2-6(b).  The  object  once  reduced  through  the  lens  sys- 
tem is  composed  of  4 white  stripes  1 mm  wide  and  1 mm  apart.  The  input 
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FIG.  U .2-6 (a )--Scheraatic  of  the  system  used  fcr  direct  imaging. 


aAaa^ 


SCANNING  PULSE  .I5;xs 
A-SCAN 


B-SCAN 


[~  0.5  fis 

b)--Direct  imaging  of  a 2 mm  period  pattern.  Input 
frequencies  105  and  llS  MHz. 


frequencies  are  I05  and  II5  MHz.  The  amplitude  display  of  the  output, 
called  the  A-scan  image,  shows  in  this  particular  example,  a contrast 
of  15  dB.  If  this  signal  is  used  to  modulate  the  intensity  of  the 

electron  beam  of  a CRT,  a B-scan  is  obtained.  Both  types  of  scans  are 
shown  in  Fig.  4.2-6(b). 

4.5  TRANSFORM  TECHNIQUES  APPLIED  TO  OPTICAL  IMAGING 

So  far  we  have  discussed  the  simplest  situation,  that  of  a short 
pulse  scan  of  the  Image.  In  this  scan  mode,  a high  resolution  requires 
that  the  time  taken  to  scan  one  spot  In  the  Image  be  short,  which  de- 
creases the  output  power,  and  there  Is  a trade-off  between  slgnal-to- 
nolse  ratio  and  resolution.  We  now  consider  the  possibility  of  using 
more  complicated  wavefoms  of  long-time  duration  (the  line  time)  to 
scan  one  line  of  the  Image.  In  this  case,  the  output  will  be  trans- 
form of  the  Illumination  distribution.  If  we  feed  this  output  Into  a 
correctly  designed  matched  filter,  as  In  radar  systems.  It  Is  possible 
to  obtain  a short  pulse  signal  corresponding  to  a spot  In  the  Image, 
I.e.,  to  make  an  Inverse  transform  of  the  scanning  signal.  As  In 
radar,  such  matched  filter  pulse  compression  techniques  give  rise  to 
an  Improvement  In  the  signal-to-nolse  ratio  because  any  spot  In  the 
Image  Is  now  scanned  for  the  line  time  rather  than  the  spot  time. 
Besides,  according  to  the  type  of  Input  signals  used,  generation  of 
various  transforms  with  different  properties  will  be  possible.  As  a 
goal,  the  chosen  transform  should  reduce  the  correlation  In  the  Image 
so  that  the  quantity  of  Information  to  be  transmitted  Is  lowered  and 
hence  the  required  bandwidth.  We  have  employed  the  Fourier  and  the 
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Fresnel  Transforms  because  they  are  easiest  to  implement  with  simple 

linear  FM  chirp  signals.  Other  digital  and  analog  transforms  are  also 

58 

of  interest."^ 


^•3*1  Spatial  Fresnel  Transform  and  Imaging 


Let  us  consider  the  complete  system  represented  in  Fig.  U.3-I. 

Two  linear  FM  chirp  signals  are  inserted  into  opposite  ends  of  the 

semiconductor  device,  with  frequencies  co=tD,  + ut 

a 1 ^1  ^ ^ “^2  ^2 

respectively;  they  give  rise  to  two  waves  with  linear  variation  of 
frequency  and  a square  law  variation  of  phase.  Using  the  center  of 

the  device  as  a reference,  we  find  that  the  phase  of  the  resultant 
product  signal  is 


> 


0(z,t)  - + = (cDj+co  )t  + (cc  -Ui.)  J ^ (t  ^ 

^2  2 

(4.3-1) 

With  a grating  filter  with  K = (a,^  - )/v  , the  ontput  signal 

due  to  a carrier  density  variation  n(z)  along  the  length  of  the 
devices,  caused  by  illumination  with  light  is 


F(t)  = e 


“/■ 


(z)  e 


(4-3-2) 


When 

a chirp  output  with 


, a Fresnel  transform  is  obtained,  in  this  case, 
a chirp  rate  > a chirp  duration  T and 
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FIG.  4.3-I--A  schematic  of  the  complete  imaging  system  using  Fresnel 
transforms . 
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a bandwidth  B = (n-  - |i  ) is  obtained  from  one  spot  at  z in 

the  image  [let  n(z)  = ^(^q)  integral  of  Eq . (4. 3-2)].  The 

image  can  be  reconstructed  by  inserting  the  output  signal  into  a suit- 
ably matched  filter,  matched  to  a chirp  rate 

One  advantage  of  the  use  of  the  grating  filter  is  the  elimination 
of  spurious  signals  due  to  multiple  reflections.  The  forward  propa- 
gating wave,  whose  electric  field  varies  as  exp[  j (tUj^t  - kj^z)]  , can 
be  partially  reflected  by  the  end  transducers  or  by  the  edges  of  the 
semiconductor,  because  of  the  electrical  or  mechanical  discontinuities 
they  introduce.  A backward  wave  of  amplitude  exp  [ j (cnj^t  + kj^z  )]  is 
therefore  excited.  The  combination  of  these  two  waves  will  generate 
a nonlinear  potential  which  varies  as  p exp  (2joij^t)  . When  the  con- 
volver is  operated  in  a degenerate  mode,  2(x>^  is  the  frequency  of  the 
detected  output  and  this  spurious  potential  is  detected.  This  signal 
corresponds  to  the  autoconvolution  product  of  the  input  at  o)^  .It 
is  found  to  be  -20  log  |p(  dB  below  the  true  convolution  output. 

By  contrast,  when  the  convolver  is  operated  in  a degenerate  mode,  only 
the  potential  which  varies  as  exp[  j (cOj^  + )t  ± j(k2"kj^)z]  are 
detected.  The  previously  described  spurious  signal  will  therefore 
not  interfere.  The  major  disturbance  now  comes  from  the  convolution 
product  of  the  first  reflections  Ot  the  two  input  waves.  Because  two 
reflections  are  involved  in  the  process,  this  spurious  convolution 
signal  is  -ko  log  |p  | dB  below  the  -.oain  signal.  In  practice  -4o  log  |p  | 
is  IfO  or  nore,  hence,  multiple  transit  convolutions  will  not  be  a severe 
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problem  for  the  nondeganerate  system,  consequently,  signals  of  arbitrary 
length  can,  and  have  been  used. 

In  Appendix  A,  a rigorous  mathematical  derivation  of  the  entire 
process  of  obtaining  an  optical  transform  is  given.  It  is  shown  how 
the  chirped  signals  have  to  be  chosen  to  achieve  a given  resolution. 

One  conclusion  is  that  if  we  suppose  that  the  time  length  of  the  chirps 
) bhe  minimum  resolvable  spot,  defined  as  the  distance  between 
dB  points,  equivalent  to  the  Rayleigh  criterion,  is  of  length  d 
I where 


d^  = 27Tv/(^^4-^2)Tc  = v/B^  . (1..5.5) 

Thus  B^  is  the  total  bandwidth  of  the  chirps,  or  the  bandwidth 
required  for  the  pulse  scan  of  the  minimum  resolvable  spot.  This 
does  not  take  into  account  the  resolution  limitation  implied  by  the 
optical  grating.  If  the  grating  periodicity  is  f , and  the  fre- 
quency offset  AF  = v/i  , the  resolution  limit  due  to  the  grating, 
d , is  of  the  order  of  one  period,  i.e., 

O 

dg  ~ ^ = v/AF  . 


Due  to  the  nature  of  the  input  chirped  signals  used,  the  difference 
between  the  instantaneous  frequencies  of  the  two  signals  varies  from 
^ ^ + ^q/2.  , where  AF  is  the  center  frequency  offset 

to  which  the  grating  periodicity  is  matched.  The  dark  signal  reduc- 


by  more  than  3 dB  from  R(AF)  . It  follows  that  B_  should  not 
exceed  AF  . Referring  back  to  Eqs.  (I1.5-5)  and  (k  , we  conclude 
that  the  resolution  limit  introduced  by  the  grating  is  comparable  to 
the  limit  inherent  in  the  transform  system  d ~ d . The  grating 
will,  in  addition,  have  the  effect  of  widening  the  apparent  width  of 
small  details  by  an  average  of  one  half  period  t /2  . The  practical 
minimum  resolution  is  therefore  approximately  d I.5  i 

Another  result  derived  in  Appendix  A is  that  the  effective  scan 
velocity  along  the  convolver  becomes 

Vs  = v(p^-m2)/(p^+P2)  • (^'5-5) 

Thus  an  arbitrary  scan  velocity  can  be  obtained.  The  number  of  resolv- 
able spots  is  still  N = L/d^  , where  L is  the  length  of  the  imaging 
device. 

It  is  also  shown  in  Appendix  A that  the  optimum  signal-to-noise 
ratio  improvement  in  using  this  matched  filter  chirped  scan  is  equal  to 
the  time  bandwidth  product  of  system,  or 

G = N(B^,/b^)  = N(v/v^)  = Tj,v/d^  = T^B^  . (U.5-6) 

There  is  an  improvement  of  B^/B^  because  the  bandwidth  of  the  noise 
entering  the  matched  filter  is  decreased  by  that  factor,  and  there  is 
a factor  N , the  number  of  resolvable  spots  which  arises  because  of 
the  large  increase  in  the  time  spent  obtaining  a signal  from  one  par- 
ticular spot  in  the  image. 

For  a TV  line  with  T^,  = 60  |isec  , B^  = 5 MHz  , N = 50O  , we 

would  use  an  acoustic  device  with  a delay  of  T = 6 psec  and  a 

A 
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chirp  bandwidth  = 50  MHz  ; the  signal-to-noise  ratio  improvement 
would  be  60  X 50  , i .e. , 35  dB . 

Several  experiments  have  been  conducted  with  the  same  optical 
arrangement  as  in  Fig.  h.2-6.  In  a first  stage,  we  used  grating 
filter  with  a periodicity  of  350  pm  (AF  = 10  MHz)  and  the  suboptimal 
figures  = 20  ps  , = 4 ps  , = I.5  MHz,  B^  = 7-5  HHz  leading 

to  a theoretical  30  point  resolution.  The  dispersive  filter  was  an 
acoustic  strip  line  waveguide  with  a total  dispersion  of  4o  ps  over 
3 MHz.  Half  of  the  dispersion,  i.e.,  20  ps  over  I.5  MHz,  was  actually 
used  in  these  experiments.  Examples  of  pictures  taken  with  this  device, 
by  mechanically  scanning  in  one  direction,  are  shown  in  Fig.  4.5-2.  The 
number  of  lines  is  approximately  20  and  the  experimental  resolution 
approximately  30  points,  i.e.,  d^  = 0.5  mm  . This  waij  determined  by- 
noticing  that  the  center  portion  of  the  "A"  of  "TAG"  is  well  detected. 
This  portion  measures  0. 5 mm  and  the  total  field  of  view  is  14  mm  wide. 
The  Output  signal-to-noise  ratio  is  20  dB  but  nonuniformities  in  the 
interactions  region  and  in  the  response  to  light  limit  the  grey  scale 
dynamic  range  of  the  picture. 

In  a later  stage,  we  built  a convolver  whose  input  interdigital 
transducers  have  3 dB  bandwidths  of  I08  ± 11  MHz  and  I54  ± 20  MHz, 
respectively.  A different  grating  of  periodicity  175  pm  was  used, 
allowing  for  a frequency  offset  AF  = 20  MHz  . With  this  larger 
offset,  it  is  possible  to  increase  the  chirp  bandwidths  and  hence 
the  resolution.  Examples  of  two-dimensional  pictures  taken  with 
this  new  system  are  shown  in  Fig.  4.3-3.  The  chirp  parameters  are 


FIG.  it  .3-2--Pictures  taken  with  the  system  shown  in  Fig.  I1.5-I. 

The  chirp  excursions  are  )i.5  and  5.0  MHz  for  20  ps, 
the  resolution  is  approximately  50  points,  the 
number  of  lines  is  approximately  20. 
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^•3-3-- (a)  Two-dimensional  picture  taken  with  the  Fresnel 
transform  imaging  system.  The  chirp  excursions 
are  9*8  and  T ,G  MHz  for  25  ps . The  minimum 
resolvable  spot  is  d^  ~ 0.25  (50  points). 

(^)j(c),(d)  Effect  of  bandwidth  limitations  of  the 
Transform  signal  on  the  reconstituted  image. 

(b ) no  limitation;  (c)  suppression  of  the 
lower  frequencies;  (d)  suppression  of  the  higher 
frequencies . 


as  follows:  = 25  |is  , = 17,5  MHz  , = 2.1  MHz  and  = 3-5  HS 

(corresponding  to  12. 5 mm  of  semiconductor  actually  being  used).  The 
minimum  resolvable  spot  is  then  d^  = 0.20  mm  and  the  number  of  points 
is  50.  Referring  to  Fig,  the  center  triangle  of  the  letter 

"A"  is  seen  to  be  well  resolved.  Since  its  base  is  0,65  mm  we  can  indeed 
estimate  the  minimum  resolvable  spot  at  0.2  mm.  The  resolution  limit  im- 
plied by  the  grating  is  slightly  higher,  d = I.5  x 0,175  mm  0.25 

With  this  resolution,  0,25  n™  > and  3 cm  long  silicon  samples, 

it  is  possible  to  obtain  pictures  with  a 120  point  resolution.  The  uni- 
formity of  the  interaction  under  such  long  samples  requires  a technology 
to  keep  the  airgap  uniform;^^  in  addition,  the  projected  images  have  to 
be  free  of  aberrations.  If  the  object  is  a standard  55  mm  slide,  the 
projection  ratio  is  close  to  one  so  a good  "relay"  lens  must  be  used 
and  the  illumination  of  the  object  has  to  be  uniform.  A diffuser  in 
front  of  the  light  source  will  achieve  this  uniformity  simply. 

In  the  Fresnel  transform  signal , each  spot  of  the  object  is  repre- 
sented by  a chirp  whose  center  frequency  is  a linear  function  of  th 

position  of  the  spot  and  whose  excursion  is  fixed  and  equals  B [see 

s 

Appendix  A,  Eq.  (A-1|)],  If  this  signal,  i.e.,  the  output  of  the  con- 
volver, is  fed  into  a bandwidth-limited  filter,  the  portions  of  the 
reconstituted  image  corresponding  to  the  suppressed  frequencies  will 
be  lost,  in  addition,  the  chirp  excursion  will  be  reduced  and  the  com- 
pression that  takes  place  in  the  matched  filter  will  not  be  as  effi- 
cient; therefore  the  quality  of  the  image  worsens,  the  resolution 
decreases  [see  Eq.  (A-lk)],  and  so  does  the  signal-to-noise  ratio 
[see  Eq.(A-19)].  Experimental  demonstration  of  this  phenomenon  is 
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shown  in  Fig,  U.J-J  (b,c,d).  In  Fig.  the  total  output  hand- 

width  ^ 2,1  MHz  is  noc  limited;  in  (c)  the  lower  0,5  MHz  of  the 
spectrum  is  cut-off,  and  in  fd),  the  upper  0.7  MHz  are  cut-off.  As 
predicted,  portions  of  the  image  are  lost  and  the  image  quality  is 
somewhat  deteriorated.  In  Fig.  .3-3(b ,c ,d ) , the  vertical  synchroniza- 
tion is  not  as  good  as  in  Fig.  i4^.3-3(a),  and  this  accounts  for  the 
vertical  blur. 


4.3.2  Spatial  Fourier  Transform  and  Imagine 


When  we  let  p,  - - p,  in  Eq . (4.3”2),  the  output  as  a function 

of  cime  is  the  spatial  Fourier  transform  of  the  carrier  density,  and 
hence  of  one  line  of  the  image.  A more  physical  view  of  the  process 
is  schematized  in  Fig.  4.3-4.  The  variation  of  the  propagation  con- 
stant ^ and  ^ of  the  two  acoustic  pulse  waves  are  repre- 
sented as  funct:ions  of  position  at  a given  time  t . The  difference 
between  the  propagation  constants  is  found  to  be  given  by  the  rela- 
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(4.3-7) 


It  is  seen  to  be  invariant  with  z ; so  the  output  electrode  picks 
the  Fourier  component  with  the  spatial  frequency  ^ . The  output  is 

given  by  a relation  of  the  form 

^op  =yn(z)  exp(j^z)  dz  = Jn{z)  exp(-2j  ^ z ) dz  , (4.3-8) 

where  n(z)  represents  the  nonlinear  coupling  variation  along  the 
semiconductor.  Here  n(z)  corresponds  to  the  optical  image  focused 
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or  a spectrum 


■,0 

on  the  semiconductor.  By  using  an  acoustic  convolver, 
analyzer,  and  inserting  this  signal  in  the  proper  way,  an  inverse 
Fourier  transform  of  the  input  signal  can  be  obtained,  and  hence  the 
original  image  signal  recovered.  This  device  offers  the  further  impor- 
tant advantage  that  it  tends  to  differentiate  against  impulse  noise 
because  impulse  noise  modulation  of  the  Fourier  transform  signal  gives 
rise  to  a uniform  background  on  the  image  rather  than  points  of  light. 

An  acoustic  convolver  can  be  utilized  to  take  the  Direct  Fourier 
transform  of  an  electrical  signal.  The  first  demonstration  of  this 

*2r  O 

operation  was  done  by  Otto.^  Two  linear  FM  chirps  of  opposite  slopes 
are  fed  into  the  convolver  and  the  signal  to  be  analyzed,  modulates  the 
amplitude  of  one  of  the  chirps.  If  we  call  f(t)  this  signal,  the  con- 
volution output  is  given  by  an  expression  similar  to  Eq . (A-20) 

r'^ 

h(t)  = exp(j<%t)  / f(t-z/v)  exp(-j2|j,tz/v)  dz 

with  the  change  of  variable  x = t - z/v  , and  the  notation  l/v  = T.  . 

/ A ' 

h(t)  can  be  put  in  the  following  form; 


h(t)  . cxp[J(iJj^  - / f(T)  exp(-j2ijtT)  dr  (li.J-lO) 

t<I^/2 

The  integration  process  taking  place  in  the  convolver  is  seen  to  be 
identical  to  the  one  of  a spectrum  analyzer.  It  is  to  be  noted  that 
when  the  integration  is  completed,  the  exponent  -2|it^  drops  out,  so 
h(t)  represents  the  Direct  Fourier  transform  of  f(t) 

Appendix  A gives  a derivation  of  the  Fourier  transform  imaging 
process.  It  is  shown  that  the  values  of  the  resolution,  the  scan 
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velocity  and  the  signal-to-noise  ratio  enhancement  are  still  given  by 
Eqs  . (U.5-5)>  and  (h  .^-6) , where  the  filter  chirp  rate  " 1^2 

is  replaced  by  the  spectrum  analyzer  internal  frequency  scanning  rate  a 
It  will  be  noted  that  in  this  case,  the  chirp  lengths  can,  in  principle, 
be  arbitrarily  long,  while  the  total  input  bandwidth  is  dictated  by  the 
resolution  required.  The  output  bandwidth  from  the  device,  however,  is 
reduced  by  the  ratio  of  the  delay  time  in  the  convolver  to  the 

length  of  the  chirp  T (the  scan  time).  We  have  carried  out  several 
experiments  to  demonstrate  these  principles  and  check  these  theoretical 
results . 

A first  example  is  illustrated  in  Fig.  4.5-5.  A grating  filter 
with  a 350  period  is  used  in  front  of  the  device  to  give  good  dynamic 
range,  approximately  25  dB . An  object  consisting  of  a photographic 
transparency  is  also  placed  in  the  path  of  the  incident  light,  and  a 
signal  corresponding  to  the  Fourier  transform  of  the  object  is  produced 
as  shown  in  Fig.  4.3~6(a).  In  this  case,  the  input  signals  are  10  ms 
long  chirps  of  opposite  sign  with  center  frequencies  of  II5  and  105  MHz, 
and  frequency  excursion  of  2.5  MHz.  In  the  example  shown,  we  used  a 
simple  7 period  grating  with  a periodic  spacing  of  approximately  2 mm 
as  the  illuminated  object,  and  displayed  the  output  amplitude  as  a 
function  of  time,  as  shown  in  the  figure.  The  center  peak  corresponds 
to  the  zero  spatial  frequency  component  and  the  two  main  lobes  to  the 
periodicity  of  the  grating. 

Using  the  very  slow  scan,  the  simplest  wav  to  reconstitute  the 
image  is  with  a spectrum  analyzer.  The  output  from  the  spectrum 
analyzer  is  shown  in  Fig.  4.5-6(b),  and  is  seen  to  correspond  well  to 
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FIG.  ^.3"5  Schematic  of  the  system  used  for  taking  a Fourier  transform. 


FOURIER  TRANSFORM 


2 ms 


FIG.  4.3-6(a)-“  A Fourier  transform  of  a periodic  grating  with  a period  of 
2 mm.  A grating  filter  of  550  (:m  spacing  was  used  with 
input  frequencies  of  II5  and  I05  MHz  and  chirps  with  a total 
excursion  of  2.5  MHz  and  10  ms. 


(b) 

INVERSE  TRANSFORM 
= IMAGE 


FIG.  4.3"6(b )--The  inverse  Fourier  transform  taken  with  a spectrum  analyzer. 


•r  • T jnr  ~~  1 


' 


the  configuration  of  the  original  object.  It  is  apparent  that  the 
output  bandwidth  is  determined  by  the  chirp  rate,  and  that  the  Fourier 
transform  technique  samples  any  point  on  the  image  for  the  whole  length 
of  the  chirp.  Thus  the  use  of  this  technique  provides  a considerable 
enhancement  of  the  signal-to-noise  ratio,  as  compared  to  employing  a 
direct  narrow  pulse  scan.  To  demonstrate  this  feature,  it  is  inter- 
esting to  compare  the  output  shown  in  Fig.  h.J>-6{h)  to  the  output 
obtained  with  the  direct  imaging  system  described  in  Section  h.2.'^ 
and  represented  in  Fig.  h.2-6{h).  The  objects  have  the  same  periodic- 
ity in  both  cases  and  the  convolvers  are  identical . It  is  evident  that 
the  resolution,  showing  up  in  the  sharpness  of  the  picture,  and  the 

signal-to-noise  ratio  are  much  improved  by  the  use  of  the  transform 
system. 

A measure  of  the  resolution  was  carried  out  with  a similar  set-up 
and  a grating  with  a periodicity  of  i = 175  pm  . The  frequency  offset 
to  Introduce  between  the  center  frequencies  of  the  input  chirped  signals 
is  therefore  20  MHz.  The  semiconductor  was  uniformly  illuminated  over 
5 mm  and  an  opaque  strip  superimposed  in  order  to  create  a dark  slit  in 
the  optical  signal.  This  method  was  adopted  in  preference  to  the  opposite 
method  in  which  a slit  of  light  is  used  to  illuminate  the  semiconductor 
because  the  "white”  output  signal  level  is  not  very  uniform.  With  this 
method,  the  effect  of  the  light  is  more  easily  evaluated.  Figure  4.5-7(a) 
shows  the  output  amplitude  variation  with  time  for  a uniform  white  light 
illumination  of  5 mm  of  silicon.  The  output  is  not  flat  due  to  some  spa- 
tial irregularities  in  the  airgap  coupling  and  local  variations  in  light 
sensitivity.  This  a major  inconvenience  in  making  a viable  imaging 
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INIFORM  ILLUMINATION 

(a) 


0.93mm  DARK  SLIT 

(b) 


0.66mm  DARK  SLIT  0.24mm  DARK  SLIT 

(c)  (d) 


FIG.  U .3-7--Examples  of  A-scans  obtained  with  a Fourier  transform 

imaging  system.  The  total  input  bandwidth  is  Bq  = Ih  MHz. 
(a)  Uniform  illumination  of  5 of  semiconductor; 

'b^,  (c),  (d)  Superposition  of  dark  slits  of  various 
widths . 


system.  Yet,  the  study  of  the  resolution  is  not  much  affected  by  this 
problem,  if  we  look  at  the  effect  on  the  image,  of  the  introduction  of 
a dark  slit  in  the  optical  input.  We  use  Fig.  4,5-7(a)  as  a reference, 
and  by  comparing  it  to  similar  images  taken  in  the  presence  of  the  dark 
slits,  we  obtain  the  effect  of  the  slit  alone. 

If  ®C  total  input  bandwidth,  i.e.,  twice  one  chirp  band- 

width, the  minimum  resolvable  spot  is  d^  , whose  expression  is  given 
by  Eqs.  (A-14)  or  (A-25), 

dg  = v/B^  • (4.5-11) 

As  mentioned  earlier  in  Section  4.3.1,  the  maximum  bandwidth  acceptable 
not  to  deteriorate  the  signal-to-noise  ratio,  is  equal  to  the  frequency 
offset  . Hence,  the  minimum  resolvable  spot  with  a given  system 

is 

d = • (4.5-12) 

smin 

The  grating  will,  in  addition,  have  the  effect  of  widening  the  apparent 
width  of  small  details  by  an  average  of  one  half  period  i/2  . The 

practical  minimum  resolution  is  therefore  1.5i  , i.e.,  in  this  case 
0 .25  nnn. 

In  Fig.  4.2-7,  we  show  some  of  the  A-scans  obtained  with  B = 14  MHz 

c 

and  T^  = 140  ps  , with  a uniform  illumination  (a),  and  with  dark  slits 
of  various  widths  (b),  (c),  (d).  Theoretically,  one  would  expect  a 
sinX/X  type  of  decrease  due  to  the  slit.  Experimentally,  due  to  dif- 
fusion in  the  semiconductor  (MTF),  the  very  sharp  edges  in  the  focused 
image  disappear  and  as  seen  in  Fig.  4.2-7(b),  (c),  (d),  only  the  main 
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lobe  of  this  decrease  is  visible.  Its  width  varies  as  the  slit  width 
and  the  depth  of  modulation  weakens  for  the  narrowest  slit  (0.24  mm), 
whose  width  is  comparable  to  the  resolution  limit  for  this  particular 
case  (d^  = 0.25  mm).  From  these  oscillograms  and  similar  ones  obtained 
with  of  6.4  MHZ  and  l8  MHz,  we  measured  the  apparent  widths  of  the 

dark  slits  and  compared  them  to  the  actual  values.  The  results  are 
suirmarized  in  Table  4.1.  When  the  total  input  bandwidth  is  6.4  MHz, 
dg  = 0.545  mm  and  the  wider  slits  of  width  O.93  and  0.66  mm  are  well 
detected;  their  apparent  width  is  found  to  be  within  10^  of  the  true 
value.  By  contrast,  the  narrower  slits  are  not  seen  at  all,  for  their 
width  is  much  smaller  than  the  .»*esolution  limit  of  0.545  mm.  With 
larger  bandwidths,  B^  = 14  MHz  and  B^  = 18  MHz  , the  resolution 
limit  due  to  the  transform  system,  d^  , is  reduced  to  0.250  and  O.I95  mm 
respectively,  and  the  narrowest  slit  used  (0.17  mm)  is  seen.  It  appears 
wider  because  of  the  previously  mentioned  widening  effect  of  the  grating 
when  the  spot  size  becomes  of  the  order  of  the  grating  period  I . As 
seen  in  Table  4.1,  the  apparent  width  is  comparable  to  what  we  define 
as  the  "artificial"  width 

'^artificial  “ ^real  ' 

As  a conclusion,  with  B^  = I8  MHz  , the  resolution  limit  is  approxi- 
mately 0.25  ram,  a result  similar  to  what  was  obtained  previously  in  the 
Fresnel  transform  system.  It  corresponds  to  O.55  mm,  i.e.  5 points  per 
mm  of  silicon,  as  had  been  expected. 
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TABLE  l+.l 


resolution  limit 


An  example  of  a picture  taken  with  a Fourier  transform  imaging 
system  by  mechanically  scanning  in  one  direction,  is  shown  in 
Fig.  k.2-8(a).  The  chirped  signals  have  a bandwidth  of  5«7  MHz  and 
a duration  of  hO  ps  (B^  = 11. h MHz).  The  portion  of  the  semicon- 
ductor actually  used  extends  over  11  mm,  and  corresponds  to  an 
acoustic  delay  of  T^  = 3 5 the  resulting  number  of  resolvable 

points  is  therefore  N = • T = 35 

V A 

In  an  actual  system,  the  output  of  the  imaging  device,  i.e.,  the 
Fourier  transform  signal,  is  sent  through  a transmission  channel.  The 
reconstitution  of  the  image  ~ the  decoding  by  second  Fourier  trans- 
form — is  carried  out  at  the  receiver  end  of  the  channel.  The  fre- 
quency in  the  transmitted  signal  therefore  corresponds  to  time  in  the 
decoded  image.  In  Fig.  h.2-8(b),  (c),  (d),  we  have  simulated  reduc- 
tion or  alteration  of  the  transmission  channel  bandwidth.  The  total 
bandwidth  required  for  an  accurate  reconstitution  is  800  kHz  [see  (b)]. 
Now,  if  the  lower  or  intermediary  frequencies  of  the  bandwidth  are  cut- 
off, portions  of  the  image  will  be  lost  accordingly,  but  this  will  not 
affect  the  quality  of  the  remaining  image.  This  is  a unique  feature 
of  the  Fourier  transform  system.  In  contrast,  in  a classical  direct 
imaging  system,  or  even  in  a Fresnel  Transform  system  (see  Section  k.3.1), 
a bandwidth  reduction  would  result  in  the  degradation  of  the  entire  image. 

In  some  imaging  applications,  the  quality  of  the  image  does  not  have 
to  be  uniform  over  the  entire  field  of  view;  the  center  part  should  be 
better  resolved  than  the  outer  portions.  This  is  particularly  true  in 
some  TV  tracking  systems,  and  in  such  systems,  a transform  coding  of  the 
previous  type  with  a tapered  bandwidth  would  be  suitable. 
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NO  FILTER  HI-PASS  FILTER  STOP-BAND  FILTER 

(O-BOOKHz)  (Fc  = 250KHz)  (400-750KHz) 

(b)  (c)  (d) 


FIG.  U.5-8--(a)  Two-dimensional  picture  taken  with  a Fourier  transform 
imaging  system  using  chirped  signals  of  total  bandwidth 
11. h MH^. 

(t>)j(c),(d  Demonstration  of  the  effect  on  the  reconstituted 
image  of  a frequency  filtering  of  the  Fourier  transform 
signal . 
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CHAPTER  V 
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CONCLUSIONS 

In  this  work  we  have  been  concerned  with  acoustic  surface  waves 
on  a piezoelectric  substrate  and  their  interaction  with  a neighboring 
semiconductor  separated  by  a small  airgap.  The  emphases  has  been  put 
on  the  capabilities  of  the  system  to  generate  the  convolution  product 
of  two  signals.  It  has  been  shown  that  the  semiconductor  can  produce 
a large  nonlinea*.  coupling  between  acoustic  waves  and  that  high  effi- 
ciency signal  processors  can  be  built.  In  addition,  new  types  of 
optical  imaging  systems  have  been  demonstrated  by  making  use  of  the 
fact  that  the  semiconductor  is  a photosensitive  continuous  array  which 
can  be  scanned  by  acoustic  waves. 

We  have  presented  a theory  of  the  semiconductor-airgap-piezo- 
electric  coupling.  This  theory  is  new  in  that  it  'accounts  for  any 
depletion,  accumulation  or  inversion  layer  at  the  surface  of  the 
semiconductor.  Such  layers  exist  on  real  surfaces  because  of  the 
presence  of  surface  states,  and  they  can  be  varied  artificially  by 
applying  a bias  field  across  the  structure.  Experimental  verifica- 
tion in  the  case  of  LiNbO^~Si  systems  resulted  in  a direct  measure- 
ment of  the  surface  state  distribution  and  time  constant  on  Silicon, 
and  in  the  first  demonstration  showing  that  acoustic  field  techniques 
are  potential  tools  for  the  investigation  of  semiconductor  surfaces. 

It  has  also  been  ascertained  that  the  convolver  efficiency  is  maxi- 
mum when  a depletion  layer  is  present  at  the  surface  of  the  semi- 
conductor, a result  which  has  been  left  out  in  previous  theories. 


I 
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We  have  showe  that  th^  study  of  the  convolution  output  can  result 
in  the  direct  measurement  of  semiconductor  parameters  which  are  other- 
wise difficult  to  evaluate.  This  measurement  is  also  non-destructive 
since  no  electrode  is  deposited  on  the  surface  of  the  semiconductor 
and  it  is  sinple  to  change  samples.  To  use  the  acoustic  device,  in 
its  present  form  as  an  investigation  tool  for  semiconductor  surfaces, 
one  would  have  to  employ  small  semiconductor  samples  to  ensure  a uni- 
form inceraction  between  the  piezoelectric  and  semiconductor  media. 
Particular  care  should  be  paid  to  the  preparation  and  mounting  of  the 
samples,  and  a study  under  vacuum  is  advisable. 

We  have  described,  at  length,  the  effect  of  a bias  voltage  and 
how  it  can  be  interpreted.  A study  versus  controlled  illumination 
and  temperature  would  also  be  of  interest  if  more  detailed  information 
is  needed.  Such  a study  was  beyond  the  scope  of  this  work,  but  it 
could  be  the  nucleus  of  some  future  work  in  the  field  of  semiconductor 
physics . 

The  characteristics  of  the  system  as  a signal  processor  hav/e  been 
studied.  A convolution  efficiency  of  -W  dBm,  one  of  the  highest  values 
reported,  with  time-bandwidth  product  in  the  200  range,  have  been  ob- 
served. The  results  are  in  good  agreement  with  theory.  Several  signal 
processing  applications  have  been  demonstrated.  We  have  shown  that  the 
convolver  is,  in  essence,  a programmable  linear  filter.  As  such,  it  has 
been  used  as  a compression  filter  or  a matched  filter  for  the  recogni- 
tion of  digital  and  analog  codes.  Due  to  the  parametric  nature  of  the 
wave  coupling  in  the  system,  time  inversion  of  a signal  is  possible  and 
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autocorrelation  have  been  demonstrated  for  the  first  time  by  cascading 

two  devices.  Other  parametric  applications  have  been  studied,  to  a 
lesser  extent. 

We  employed  the  photoconductivity  of  the  semiconductor  to  turn  the 
convolution  device  into  an  imaging  device.  By  introducing  a periodicity 
in  the  optical  signal  and  operating  the  system  with  offset  input  center 
frequencies,  we  have  shown  low  dark  current  and  large  dynamic  range 
operation.  These  improvements  have  been  shown  to  agree  well  with  theo- 
retical calculations.  In  addition,  a now  method  of  scanning  the  image 
has  been  demonstrated.  Long  time  signals  are  fed  into  the  acoustic  con- 
volver and  a real  time  Fourier  or  Fresnel  transform  of  the  image  is  gen- 
erated. A major  improvement  in  signal-to-noise  ratio  has  been  obtained 
and  the  scan  velocity  has  been  made  arbitrary. 

It  should  be  noted  that  these  scanning  methods  are  not  possible 
with  electron  beam  devices  like  the  vidicons,  and  they  have  not  been 
demonstrated  with  CCDs.  They  may  very  well  not  be  applicable  to  these 
devices.  With  acoustic  systems,  these  methods  are  particularly  simple 

to  implement  and  they  increase  the  versatility  of  the  system  to  a large 
extent . 

These  new  concepts  have  all  been  tested  with  the  separated  medium 
semiconductor-piezoelectric  convolver  because  of  the  experience  we  have 
had  with  this  structure  and  because  of  its  versatility.  The  airgap  uni- 
formity and  the  semiconductor  surface  treatment  are  problems  requiring 
particular  attention  in  this  device  configuration.  Reproducible  devices 
will  have  to  be  monolithic.  Such  structures  are  being  developed  and  the 
principles  presented  here  would  apply  directly  to  them. 

- 151  - 

. 


appendix  a 


IMAGING  BY  TRANSFORM  CODING 


»o  8lv.  here  a derivation  of  the  .aathematleal  reaults  associated 
With  the  convolution  and  filtering  processes  that  take  place  in  the 
complete  imaging  systems  using  transform  techniques. 

Notation: 

L — length  of  the  semiconductor 
V = acoustic  wave  velocity 

= L/v  acoustic  delay  under  the  semiconductor. 

The  input  chirped  signals  have  the  following  characteristics; 

- duration  T 

c 

- chirp  rate  and  in  rd/s^ 

- total  input  bandwidth 

- dispersive  filter  bandwidth  B^  = ^ 

Convolver  Output: 

We  take  the  center  of  the  semiconductor  as  the  origin  a . 0 . 

If  the  t»o  input  electrical  signals  are  f^(t)  and  f^(t)  and  the 

carriers  density  is  uniform,  the  output  voltage  of  the  convolver  h(t) 
is  given  by 


L/2 

- z/v)  f^Ct  + z/v)  dz  . 

-L/2 
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Now  let  and  be  the  chirped  signals  of  center  frequency 

and  and  chirp  rate  and  , respectively: 


/ "^l*^  \ 

fj^(t)  = exp  j + 1 t , 


(A-2) 


and 


exp  J (Oig  - 


(A-5) 


Due  to  the  presence  of  the  grating,  we  assume  no  dark  current  and 
hence  the  output  of  the  convolver,  from  any  point  z in  the  sample, 
proportional  to  the  light  generated  carrier  density  n(z)  and 
express  this  output  as 


IS 


we 


L/2 


h(t)  = 


= exp(jcjUjt)  n(z)  exp  j 

-L/2 


' ^2  i 2 


It  - 


1^1  + lig 


zt 


(A-4) 


The  output  corresponding  to  one  spot  of  light  at  z^  is  obtained  by 
letting  n(z)  = 5(z  - z^)  in  Eq . (A-h) 


h (t)  = exp  j 


- 


1^1  + 1^2  . >^1  - 1^2 


t X exp  j 
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A* 


h^Ct)  is  an  FM  chirp  with  a chirp  rate  and  its  center 

frequency  is  the  linear  function  of  position 


03(Zo)  = 0, z 

V 


(A-5) 


can  also  be  expressed  as  follows 


h^Ct)  = exp 


t^l  - ^2  / 2 \ 2^1  p 


where  we  note  the  time  delay  z/v  v^ere  \ is  defined  by 

5 c •' 


(A-6) 


relation 


^1  - ^^2 
^1  ^ >^2 


(A-7) 


To  reconstitute  the  original  image  n(z)  , we  must  thus  feed  h(t) 
into  a filter  matched  to  a chirp  rate  - p^  and  a chirp  length 
Tj,  . Such  a filter  must  have  an  impulse  response  proportional  to 
h(-t)  ; this  impulse  response,  denoted  f(t)^  is  therefore  given  by 


the  relation 


' \ ” ^2  P / ^ \ 

= A. exp  j 0)  t + t • IT  / — 1 

r ^ J 


(A-8) 


I5i^  - 


FT. 


where  II  is 


unity  in  the  interval  [-T^/2,  T^/2]  and  zero  out- 
side it;  A is  a normalization  factor  obtained  by  writing  that  the 
amplitude  response  of  the  filter  is  unity  over  its  bandwidth,  i.e., 


/ 


f(t)  exp(jo)t)dt 


1,  for  < 0)  < 0)^  + (^^  - 

1/2 


5Q 

One  can  find  ^ that  A = (B  /T  ) ■ j-j 

' s'  C'  * uid  not  account  for  the 

time  delay  in  the  expression  for  f(t)  . The  output  of  the  filter  is 

given  by  the  convolution  integral  of  the  input  h(t)  and  the  impulse 
response  f(t)  , hence 


8(t)  = J‘ \\{-x)‘i{t  - x)dT 


(A-9) 


By  using  E,s.  (A-1.)  and  (A-8),  g(t)  can  be  expressed  as  follous 

1/2  L/2 


s(t)  ~ ^ ^ j 


t + 


Pi  - P2 


-L/2 


-I  - t 

>vJ 


t+T^/2 

f 

X j exp 

j(Pl  - ) 

t-I^/2 

\ V / 

L S 

dx 


(A-11) 


The  integral  over  x introduces  a sin  X/X  function.  After  a few 
rearrangements,  we  obtain 
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g(t) 


(z) 


X 


sin  7TB  (t  - z/v  ) 
s s 

7TB^(t  - z/Vg) 


exp 


(A-12) 


The  sin  X/X  function  is  the  "scann:ng  function."  It  is  peaked  about 
the  point  z - v^(t)  ^ hence  is  the  actual  scan  velocity  given 

by  the  relation 


V 

s 


^1  ~ ^2 
Hi  H2 


(A-15) 


The  5 dB  width  of  the  "scanning  function"  defines  the  minimum  re- 
solvable spot  d : 


d 

s 


V 


(A-lU) 


The  number  of  resolvable  spots  is  therefore  given  by 


N = 


= Va  = "s^’c 


•(A-15) 


The  last  equality  comes  from  the  definitions  of  the  velocities  v and 
v^  ; we  can  indeed  write 
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fi 


vT.  = 


For  medium-to-high  resolution  (N  > 50),  the  scanning  function 
is  sharply  peaked  about  z = v^(t)  . In  this  case  it  can  be  shown 

that  over  the  center  lobe  of  width  2d^  , the  phase  error,  in  the 

integral  of  Eq . (A-12),  is  only  tt/N  . We  may  therefore  approximate 
the  phase  term  in  this  integral  by  its  value  at  z - v^(t)  . The 

limits  of  integration  can  also  be  removed  to  infinity  if  we  consider 
that  n(z)  is  zero  outside  the  interval  [-l‘/2,  L/2’  . Finally, 

to  obtain  the  amplitude  of  the  leconstituted  image,  we  approximate 
the  scanning  function  by  the  dirac  function 


Equation  (A-16)  follows; 


g(t) 


= 4/  — V n 


(v  t)exp  - j 


+ ^^2) 


. (A-16) 


g(t)  is  a replica  of  n,  chirped  with  a chirp  rate  approximately  equal 


to  ^ center  frequency 
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pro- 


Sinv.e  the  resolutions  have  to  be  identical^  we  can  write 
N = Ta©  ^ = TaBg  . We  then  obtain  the  S/n  improvement  G 
vided  by  the  Transform  Coding  of  the  image.  G is  given  by 


G 


C C 


(A-18) 


We  used  Eq . (A-II;).  Equivalent  forms  for  G are  readily  derived: 


(A-19) 


Design  Curves 

We  denote  by  a the  dispersion  of  the  reconstitution  filter 
(in  Hz/sec),  by  definition  of  a , ^^ra  = . Besides,  we 

call  and  B^  the  frequency  excursion  in  the  two  input  chirped 

signals  so  that  2ttB^  = and  ^rrB^  = . With  these 

additional  notations  we  give  below  a summary  of  the  results  needed 
in  designing  the  system: 


“ = »1  - »2 
“C  = »1 


h - '2 


N = B T 
s C 


C A 
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We  can  express  B , T and  B in  terms  of  the  acoustic  delay 
o u s 

time  , the  filter  dispersion  a and  the  number  of  resolvable 

spots  N . We  obtain 

= (N/a)^/^ 

Bg  = (aN)^/^ 

Figure  A-l(a)  gives  plots  of  B , T and  B vs  N in  the 

U U S 

particular  experimental  case  where  T = U psec  and  a is 

A 

3 MHz/l|0  psec.  These  plots  allow  us  to  determine  B , B„  and  T , 

the  chirp  characteristics^  to  achieve  any  resolution  N below  120. 

Figure  A-l(b)  gives  plots  of  T.  and  B vs  a for  the  three 

C s 

resolutions  N = 120,  200  and  500  . If  the  line  time  T and  the 

resolution  N are  given,  one  can  therefore  find  what  dispersive 

filter  to  use. 

The  Particular  Case  of  the  Fourier  Transform  Coding 

When  the  two  input  chirped  signals  have  equal  and  opposite  chirp 
rates,  the  output  of  the  convolver  is  given  by  Eq.  (A-U)  where 

^L/2 

h(t)  = exp  (joo^t)/  n(z)  exp  (-j2pzt/v)  dz  (A-2o) 

-L/2 
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FIG.  A-1  Design  curves  for  the  imaging  system  by  Fresnel  transforms. 

(a)  Line  time,  filter  bandwidth  and  chirp  bandwidth 

= 4 |is  and  a = 5 MHz/iiO  \is  . 

(b)  Line  time  and  filter  bandwidth  vs  filter  dispersion 
for  various  numbers  of  resolvable  spots. 
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FILTER  BANDWIDTH  BjCMHz) 
CHIRP  BANDWIDTH  B^UlOMHi) 


oi(z)  = - 2nz/v  (A-21) 

Reconstitution  of  the  image  is  done  by  feeding  h(t)  into  a spectrum 
analyzer  of  center  frequency  . If  T is  the  integration  time 

and  a the  frequency  scanning  rate  (in  Hz/s)  in  the  spectrum  analyzer, 
its  output  is  given  by 


X exp  - 27rat)t^  dt^dz 

where  is  the  normalization  factor.  The  integration  over  t^ 

is  straightforward  and  it  yields  the  following  result 


This  result  is  the  equivalent  of  Eq • (A-12)  and  we  can  draw  similar 

conclusions.  The  actual  scan  velocity  v , the  minimum  resolvable 

s 

spot  d^  and  the  number  of  resolvable  spots  N are  respectively 


Vg  = vn-a/(i 

d^  = 7T  v/^if  (A-2h) 

N = L/dg  = T^hT/tt 

To  obtain  the  amplitude  of  the  ■reconstituted  imuge^  we  approximate 
the  sin  X/X  function  of  Eq.  (A-25)  by  the  Dirac  function 

— 6(z  - v t) 

m s ^ 

aT 


Eq.  (A-25)  follows 


1 


g(t) 


V n(v  t) 
s ^ s 


(A-25) 


The  noise  in  the  transform  system  is  proportional  to  the  spectrum 
analyzer  bandwidth  aT  ; the  signal  amplitude  is  given  by  Eq.  (A. 25). 
Hence 


(S/N) 


Transform  System 


(S/N) 


Direct  System 


4kT  dT 


(v9)^ 

1+kT  0‘ 
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as  we  saw  earlier.  The  resolutions  of  the  two  systems  have  to  be 
equal,  hence  N = T^|iT/7T  = T^Q  ^ ; the  S/N  improvement  G provided 
by  the  transform  coding  of  the  image  is  therefore  given  by 

G = liT^/zr 

It  will  be  noted  that  the  bandwidth  of  the  spectrum  analyzer  aT 
is  the  equivalent  of  the  filter  bandwidth  of  the  Fresnel  Transform 

system.  Furthermore  in  practice,  the  integration  time  T is  greater 
or  equal  to  the  chirp  time  T , T should  be  used  in  Eq.(A-2U) 

(A-26)  and  all  the  results  of  the  Fourier  Transform  system  reduced  to 
those  derived  earlier  for  the  Fresnel  Transform  system. 


(A-26) 
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appendix  b 


SILICON  CLEANING  AND  OXIDATION  PROCEDURES 

We  give  here  the  Silicon  cleaning  and  oxidation  procedures  we 
have  followed  in  the  preparation  of  the  silicon  wafers  which  required 
a thermally  grown  oxide  layer.  The  cleaning  procedure  was  developed 
by  B.  T.  Khuri-Yakub  at  the  Stanford  Integrated  Circuits  Laboratory. 

Wafer  Cleaning: 

Use  a clean  quartz  beaker,  clean  tweezers  and  a clean  quartz  wafer 
holder. 

1.  Boil  in  trichlorethylene  for  5 minutes. 

2.  To  remove  any  photoresist  or  wax,  heat  in  hot  H^SO^rH^O^  (5:1) 
for  10  minutes. 

5.  Rinse  in  de-ionized  (DI)  HgO  for  at  least  1 minute. 

ll.  Submerge  in  warm  solution  of  H20:H202:NHj^0H  (5:1:1)  for  5 to  10 

minutes  (while  in  the  above  solution,  scrub  with  Q-tip  only  if  needed) 
5-  Rinse  in  DI  HgO  for  at  least  1 minute. 

6.  Submerge  in  warm  solution  of  H20:H202:HC1  (5:l:i)  for  5 to  10 

minutes . 

7.  Rinse  in  DI  HgO  for  at  least  1 minute. 

8.  Dip  in  50:1  - H20:HF  solution  for  30  seconds. 

9.  Rinse  thoroughly  in  DI  H^O  for  at  least  5 minutes;  then  rinse 
several  times  in  double  distilled  H^O  . Blow  dry  with  , and 
process.  Do  not  store. 
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Thermal  Oxidation: 


- 15  minutes  in  wet  O2  or  50  minutes  in  dry  O2  1200  C to 
obtain  approximately  1000  X of  oxide. 


- Then  5 minutes  annealing  in  at  1200  C 
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